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PREFACE 


This  report  was  prepared  by  the  Department  of  Civil  and  Environmental  Engineering, 
he  University  of  Rhode  Island,  Kingston,  Rhode  Island,  under  Contract  No.  F49620-88- 
C00523/SB588 1-0378  for  the  Research  Initiation  Program  (RIP)  sponsored  by  the  Air 
Force  Office  of  Scientific  Research  (AFOSR)  for  the  Engineering  and  Services  Laboratory, 
Headquarters,  Air  Force  Engineering  and  Services  Center  (HQ  AFESC/RDCM),  Tyndall 
Air  Force  Base,  Florida. 

This  report  is  published  by  the  Civil  Engineering  Research  Division  (CERD),  of 
Research  Development  and  Acquisition  Directorate  (RDAD),  Air  Force  Civil  Engineering 
Support  Agency  (AFCESA),  139  Barnes  Drive,  Tyndall  Air  Force  Base,  Florida  32403- 
5323. 

This  report  summarizes  the  results  of  work  to  investigate  the  dynamic  behavior  of 
compacted  unsaturated  soil  using  the  Split-Hopkinson  Pressure  (SHPB)  device  and 
represents  the  Master’s  thesis  research  of  Mr.  Blaise  J.  Fitzpatrick,  Department  of  Civil  and 
Environmental  Engineering,  University  of  Rhode  Island,  Kingston,  Rhode  Island..  A 
detailed  microstructural  analysis  was  performed  on  epoxied  specimens  of  compacted 
unsaturated  soil  to  study  the  relationship  between  microstmcture  and  macroscale  behavior. 
The  work  was  initiated  in  January,  1990  and  completed  in  December,  1991.  Dr.  George 
E.  Veyera  served  as  the  principal  investigator  at  The  University  of  Rhode  Island.  Mr. 
Charles  R.  Bailey  served  as  the  project  officer  for  HQ  AFESC/RDCM. 
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EXECUTIVE  SUMMARY 


A.  OBJECTIVE 

The  main  objective  of  this  study  was  to  determine  how  the  moisture  content  during 
compaction  affects  soil  properties,  including  compressional  wave  speed,  stress 
transmission  ratio,  quasi-static  constrained  modulus,  and  soil  microstructure.  This 
investigation  is  directed  at  developing  a  better  understanding  of  unsaturated  soil  behavior 
and  evaluates  the  influence  of  microsiructure  on  soil  behavior  under  dynamic  and  quasi¬ 
static  loading  conditions. 

B.  BACKGROUND 

The  ability  of  a  soil  to  transmit  dynamic  stresses  (energy)  is  of  particular  interest  to 
the  U.S.  Air  Force  with  respect  to  military  protection  construction  and  survivability 
designs.  The  prediction  of  ground  motions  from  conventional  weapons  detonations  and 
the  effect  on  structures  requires  information  about  the  response  of  geologic  materials  to 
intense  transient  loadings.  The  general  nature  of  stress  wave  propagation  and  energy 
transmission  in  particulate  media  such  as  soil  depends  on  a  number  of  parameters,  the 
interrelationship  between  which  is  not  fully  understood.  The  Split-Hopkinson  Pressure 
Bar  (SHPB)  device  provides  a  laboratory  experimental  capability  for  simulating  the  field 
explosive  loading  environment  and  providing  insight  into  the  dynamic  behavior  of  soils. 

Current  design  procedures  generally  assume  dry  soil  conditions  even  though  soils  in 
situ  are  actually  at  a  condition  somewhere  between  dry  and  saturated.  In  addition,  field 
compaction  specifications  for  granular  soils  only  require  a  final  compacted  density  to  be 
achieved  without  any  restriction  on  moisture.  However,  for  unsaturated  conditions,  the 
amount  of  energy  traveling  through  a  soil  and  transmitted  to  a  structure  can  be  significantly 
larger  than  for  dry  soil  conditions,  thereby  having  the  potential  for  greater  damage. 
Unsaturated  soils  are  complex,  multiphase  materials  the  behavior  of  which  is  not  well 
understood.  Previous  research  has  shown  that  compacting  moist  sands  with  varying 
amounts  of  water  prior  to  dynamic  loading  can  increase  the  stress  transmission  ratio  by  as 
much  as  a  factor  of  two.  However,  at  present,  no  theoretical,  empirical  or  numerical 
methods  are  available  for  predicting  large  amplitude  compressive  stress  wave  transmission 
and  compressive  stress  wave  propagation  velocity  in  unsaturated  soils.  A  need  exists  to 
develop  a  better  understanding  of  unsaturated  soil  behavior  at  a  fundamental  level  for  all 
types  of  loadings.  Therefore,  a  comprehensive  experimental  laboratory  study  using  the 


SHPB  for  dynamic  loading,  a  quasi-static  uniaxial  loading  device,  and  microstmctural 
analysis  techniques  was  conducted  on  specimens  of  compacted  unsaturated  soil. 

C.  SCOPE 

Specimens  of  Ottawa  20-30  sand  were  dynamically  compacted  to  a  constant  dry 
density  at  various  saturations  using  a  standard  Proctor  hammer  and  subjected  to  high  strain 
rate  testing  using  a  Split-Hopkinson  Pressure  Bar  (SHPB)  device.  Saturations  were  varied 
from  0  to  about  80  percent.  Microstmctural  analyses  of  epoxied  compacted  unsaturated 
sands  were  also  performed.  In  addition,  quasi-static  uniaxial  compression  tests  were 
conducted  to  examine  compacted  stiffnesses.  The  influence  of  soil  microstmcture  on 
compaction  energy,  stress  transmission,  and  compressional  wave  speed  for  Ottawa  20-30 
sand  was  evaluated. 

D.  METHODOLOGY 

Specimens  of  Ottawa  20-30  sand  were  dynamically  compacted  to  a  constant  dry 
density  at  different  moisture  conditions  (saturations  from  0  percent  to  about  80  percent)  in  a 
thick-walled  stainless  steel  tube  using  a  standard  Proctor  hammer.  A  series  of  SHPB  tests 
was  conducted  on  the  compacted  specimens  to  provide  information  about  stress 
transmission  and  wave  speed  characteristics.  In  addition,  a  number  of  specimens  were 
epoxied  after  compaction  and  after  SHPB  testing  to  preserve  the  microstmcture  and  obtain 
data  on  variations  in  particle  orientations  due  to  the  influence  of  moisture  present  during 
compaction  to  develop  correlations  with  dynamic  behavior.  A  series  of  quasi-static  uniaxial 
confined  compression  tests  was  conducted  to  evaluate  the  stress-strain  response  as  a 
function  of  moisture  conditions. 

E .  TEST  DESCRIPTION 

The  SHBP  at  the  Air  Force  Engineering  and  Services  Center  (AFESC/RDCM)  was 
used  to  dynamically  test  50.8  mm  diameter  specimens  of  compacted,  unsaturated  sand  in  a 
thick-walled  stainless  steel  specimen  container  at  high  strain  rates.  Specimens  were 
compacted,  epoxied  and  sectioned  to  preserve  their  microstmcture  for  analysis  using 
specialized  equipment  at  AFESC/RDCM.  Quasi-static  uniaxial  confined  compression  tests 
were  also  conducted  using  a  Materials  Testing  System  (MTS)  ate  AFESC. 
Photomicrographic  studies  and  microstmctural  analyses  were  conducted  at  the  University 
of  Rhode  Island’s  Geomechanics  Research  Laboratory  in  the  Department  of  Civil  and 
Environmental  Engineering. 
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F.  RESULTS 


Previous  research  has  shown  that  compaction  energy,  constrained  modulus  (a 
measure  of  soil  stiffness),  compressional  wave  speed,  and  stress  transmission  ratio  have 
greater  values  at  intermediate  saturations  (from  about  20  percent  to  75  percent)  than  for  dry 
soil  conditions.  The  reasons  for  this  phenomena  are  not  fully  understood.  Variations  in 
soil  microstructure  developed  during  compaction  appear  to  be  a  major  factor.  The 
experimental  investigation  presented  herein  utilized  dynamic  and  quasi- static  testing 
methods  to  evaluate  the  influence  of  soil  microstructure  on  compaction  energy,  constrained 
modulus,  wave  speed,  stress  transmission  ratio,  and  degree  of  saturation  for  Ottawa  20-30 
sand. 

In  this  study,  unsaturated  sand  specimens  were  dynamically  compacted  to  a 
constant  dry  density  at  various  saturations  using  a  standard  Proctor  hammer  and  subjected 
to  high  strain  rate  testing  using  the  split-Hopkinson  pressiffe  bar  device.  Previous  research 
has  shown  that  compacting  moist  sands  with  varying  amounts  of  water  prior  to  dynamic 
loading  can  increase  the  stress  transmission  ratio  by  as  much  as  a  factor  of  two.  The 
compaction  method  and  the  amount  of  moisture  have  also  been  shown  to  influence  static 
and  dynamic  behavior  of  sands.  The  influence  of  soil  microstructure  on  compaction 
energy,  stress  transmission,  and  compressional  wave  speed  for  Ottawa  20-30  sand  was 
evaluated.  Experimental  results  suggest  that  the  dynamic  behavior  of  the  sand  is  dependent 
on  sou  microstructure  which  is  a  function  of  the  initial  naoisture  content  during  compaction. 

A  two-dimensional  grain  orientation  analysis  was  performed  on  epoxied  specimens  of 
compacted  unsaturated  sand.  Results  from  the  analysis  suggest  that  microstructure  may  be 
a  significant  factor  affecting  the  dynamic  behavior  of  unsaturated  soUs  and  showed  vector 
magnitude  (a  measure  of  dispersion  within  a  particle  orientation  data  set)  could  be  related  to 
variations  in  saturation.  At  intermediate  saturations,  vector  magnitude  values  were  lower 
than  that  for  dry  soil,  which  indicates  a  stiffer  soil  structure.  According  to  wave 
propagation  theory,  a  stiffer  material  will  produce  greater  compressional  wave  speed  and 
stress  transmission  ratio  results  at  intermediate  saturations,  as  was  the  case  for  this 
investigation. . 

Based  on  the  results  of  this  investigation,  it  was  determined  that  soil  microstructure 
does  affect  compaction  energy,  constrained  modulus,  compressional  wave  speed,  and 
stress  transmission  ratio.  However,  a  clear  and  concise  explanation  of  the  microstructural 
development  process  is  not  ciurently  available.  Characterization  of  soil  microstructure  is  an 
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important  first  step  towards  establishing  a  better  understanding  of  stress  transmission  in 
unsaturated  soils.  Relationships  between  vector  magnitude  and  dynamic  soil  parameters 
may  be  useful  in  developing  a  model  to  predict  stress  transmission  in  unsaturated  soils 
subjected  to  dynamic  loading  from  conventional  weapons. 

G.  CONCLUSIONS 

The  results  from  the  SHPB  tests,  quasi-static  uniaxial  confined  compression  tests  and 
microstructural  analyses  indicate  that  the  amount  of  moisture  present  during  compaction  is 
probably  the  most  significant  factor  influencing  both  the  static  and  dynamic  behavior  of  the 
Ottawa  20-30  sand. 

The  compactive  energy  required  to  prepare  constant  dry  density  specimens  using  a 
standard  Proctor  hammer  was  shown  to  be  strongly  dependent  on  the  amount  of  moisture 
present  during  compaction.  The  SHPB  data  show  that  the  amount  of  energy  transmitted  is 
greatest  at  saturations  in  the  range  of  about  20  to  80  percent  compared  with  values  for  dry 
sand.  The  microstructural  analysis  results  suggest  an  inverse  relationship  between  vector 
magmtude  and  compaction  energy,  wave  speed,  stress  transmission  ration  and  constrained 
modulus  with  variations  in  the  degree  of  saturation.  Preferred  particle  orientations  were 
observed  at  saturations  of  0  to  about  20  percent  and  above  75  percent,  while  a  random 
orientation  was  observed  at  saturations  in  the  range  of  about  20  to  75  percent  In  addition, 
quasi-static  uniaxial  confined  compression  test  results  show  changes  in  specimen  stiffness 
with  changes  in  compaction  moisture  which  are  generally  consistent  with  trends  observed 
in  the  dynamic  SHPB  data.  The  variations  in  dynamic  and  static  behavior  observed  in  this 
investigation  can  be  attributed  to  the  influence  of  soil  microstructure  developed  during 
compaction  the  nature  of  which  is  predominantly  a  function  of  the  amount  of  moisture 
present. 

H.  RECOMMENDATIONS 

The  results  of  this  investigation  have  provided  important  insight  into  the  complex 
dynamic  behavior  of  unsaturated  soils  about  which  very  little  is  known.  These  findings  are 
significant  with  regards  to  military  protective  construction  design  and  require  further 
detailed  studies  in  both  the  laboratory  and  field  coupled  with  numerical  model  development 

Further  laboratory  work  using  the  SHPB  should  be  conducted  to  study  other 
parameters  including  particle  size,  particle  shape  and  variations  in  compacted  dry  density 
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under  various  confinement  conditions.  In  addition,  detailed  microstructural  analyses 
should  also  be  a  part  of  this  work. 

Scaled  laboratory  tests  should  be  performed  in  the  centrifuge  using  explosive  charges 
in  compacted,  unsaturated  soils.  This  will  allow  the  simulation  of  various  field  conditions 
in  a  laboratory  environment  and  provide  independent  verification  of  the  SHPB  results.  The 
centrifuge  is  an  important  research  tool  and  has  a  distinct  advantage  over  field  tests  in  that  a 
number  of  tests  can  be  performed  in  a  short  time  period.  The  AFESC  centrifuge  facility  is 
ideally  suited  for  this  work.  In  addition,  laboratory  and  numerical  studies  should  be 
checked  against  a  series  of  carefully  controlled,  fully  instrumented  scaled  explosive  field 
tests. 

The  laboratory  work  has  progressed  to  the  point  where  the  development  of  an  initial 
preliminary  numerical  model  in  conjunction  with  further  experimental  work  is  needed.  The 
model  should  incorporate  a  mathematical  approach  such  as  the  distinct  element  or  finite 
element  technique  based  on  the  laboratory  derived  results  for  wave  speed,  transmission 
ratio,  and  vector  magnitude,  to  predict  the  compressional  wave  speed  and  stress 
transmission  in  unsaturated  soils.  This  will  prove  to  be  a  useful  tool  to  help  explain  the 
intricacies  of  the  laboratory  observed  behavior  and  also  provide  guidance  for  the  direction 
of  continued  experimental  work.  Once  centrifuge  and  field  data  are  available,  they  can  be 
used  to  further  develop  and  improve  the  naodel. 
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SECTION  I 


INTRODUCTION 


A.  OBJECTIVES 

The  main  objective  of  this  study  was  to  determine  how  the  moisture  content  during 
compaction  affects  soil  properties,  including  compressional  wave  speed,  stress 
transmission  ratio,  quasi-static  constrained  modulus,  and  soil  microstructure.  This 
investigation  is  directed  at  developing  a  better  understanding  of  unsaturated  soil  behavior 
and  evaluates  the  influence  of  soil  microstructure  under  dynamic  and  quasi-static  loading 
conditions. 

B.  BACKGROUND 

The  ability  of  a  soil  to  transmit  dynamic  stresses  (energy)  is  of  particular  interest  to 
the  U.S.  Air  Force  with  respect  to  military  protection  construction  and  survivability 
designs.  Field  compaction  specifications  for  granular  soils  only  require  a  final  compacted 
density  to  be  achieved  without  any  restrictions  on  moisture.  In  addition,  design  procedures 
generally  assume  dry  soil  conditions.  However,  for  unsaturated  conditions,  the  amount  of 
energy  travelling  through  a  soil  and  transmitted  to  a  structure  can  be  significantly  larger 
than  for  dry  soil  conditions,  thereby  having  the  potential  for  greater  damage.  Unsaturated 
soils  are  complex,  multiphase  materials  the  behavior  of  which  is  not  well  understood. 
Ross  et  al.  (1986),  state  that  "...there  are  currently  no  theoretical,  empirical  or  numerical 
methods  available  for  predicting  large  amplitude  compressive  stress  wave  velocity  and 
stress  transmission  in  unsaturated  soils."  Therefore,  a  need  exists  to  develop  a  better 
understanding  of  unsaturated  soil  behavior. 

Previous  research  has  shown  that  compaction  energy,  constrained  modulus  (a 
measure  of  soil  stiffness),  compressional  wave  speed,  and  stress  transmission  ratio  have 
greater  values  at  intermediate  saturations  (from  about  20  percent  to  75  percent)  than  for  dry 
soil  conditions.  The  reasons  for  this  phenomena  are  not  fully  understood.  Variations  in 
soil  microstructure  developed  during  compaction  appear  to  be  a  major  factor.  The 
experimental  investigation  presented  herein  utilized  dynamic  and  quasi-static  testing 
methods  to  evaluate  the  influence  of  soil  microstmcture  on  compaction  energy,  constrained 
modulus,  wave  speed,  stress  transmission  ratio,  and  degree  of  saturation  for  Ottawa  20-30 
sand.  A  laboratory  procedure  was  developed  to  prepare  specimens  of  compacted 
unsaturated  soil  for  microstructural  analysis.  Vertical  and  horizontal  sections  were  taken 
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through  the  specimen  so  that  a  detailed  microstructural  analysis  including  axial  ratio  and 
preferred  particle  orientation  could  be  performed  on  different  planes. 

Based  on  the  results  of  this  investigation,  it  was  determined  that  soil  microstructure 
does  affect  compaction  energy,  constrained  modulus,  compressional  wave  speed,  and 
stress  transmission  ratio.  Microstructural  analyses  showed  vector  magnitude  (a  measure  of 
dispersion  within  a  particle  orientation  data  set)  could  be  related  to  variations  in  saturation. 
At  intermediate  saturations,  vector  magnitude  values  were  lower  than  that  for  dry  soil, 
which  indicates  a  stiffer  soil  structure.  According  to  wave  propagation  theory,  a  stiffer 
material  will  produce  greater  compressional  wave  speed  and  stress  transmission  ratio 
results  at  intermediate  saturations,  as  was  the  case  for  this  investigation. 

Results  from  this  investigation  may  be  useful  in  future  microstructural  studies  of 
cohesionless  soils.  Relationships  between  vector  magnitude  and  dynamic  soil  parameters 
may  be  used  in  a  developing  model  to  predict  stress  transmission  in  unsaturated  soils 
subjected  to  dynamic  loading  from  conventional  weapons. 

C.  SCOPE 

Specimens  of  Ottawa  20-30  sand  were  dynamically  compacted  to  a  constant  dry 
density  at  various  saturations  using  a  standard  Proctor  hammer  and  subjected  to  high  strain 
rate  testing  using  the  Split-Hopkinson  Pressure  Bar  device.  Saturations  were  varied  from  0 
to  about  80  percent.  Microstructural  analyses  of  epoxied  compacted  unsatmated  sands 
were  also  performed.  In  addition  quasi-static  uniaxial  compression  tests  were  conducted  to 
examine  compacted  sdtiffnesses.  The  influence  of  soil  microstructure  on  compaction 
energy,  stress  transmission,  and  compressional  wave  speed  for  Ottawa  20-30  sand  was 
evaluated. 
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SECTION  n 


LITERATURE  REVIEW 

The  literature  review  contains  five  sections  which  directly  relate  to  the  present 
investigation.  Because  the  research  reported  here  deals  with  granular  soils,  the  literature 
review  will  concentrate  on  relevant  aspects  of  these  specific  soil  types.  The  first  section 
defines  fabric  and  structure  in  soils  and  explains  the  difference  between  soil  fabric  and  soil 
structure.  The  second  section  deals  with  cunent  methods  of  microstructural  analysis.  The 
third  section  explains  sample  preparation  techniques  and  statistical  analysis  of  two- 
dimensional  particle  orientation  data.  SECTION  IV  summarizes  soil  compaction  principles 
and  investigates  compaction  energy  in  unsaturated  soils.  The  final  part  of  SECTION  n  is  a 
review  of  stress  wave  propagation  theory  and  the  effects  of  compressive  wave  propagation 
through  unsaturated  soils. 

A.  STRUCTURE  AND  FABRIC  IN  SOILS 

The  development  of  fabric  and  structure  in  soil  is  a  function  of  natural  deposition  and 
stress  history  in  the  field  and  san^le  preparation  method  in  the  laboratory.  Particle  size  and 
shape  play  an  important  role  in  the  development  of  soil  structure.  In  cohesionless  soils, 
angularity  and  sphericity  are  used  to  describe  particle  size  and  shape.  These  two 
parameters  are  very  useful  in  a  microstructural  analysis.  Shear  strength,  compressibility, 
and  mechaiucal  behavior  of  soils  are  all  affected  by  fabric/strucmre  development  and  the 
type  of  applied  stress,  i.e.,  static  or  dynamic. 

1.  Definition  ofFabric  and  Structure 

The  terms  "soil  fabric"  and  "soil  structure"  have  been  used  interchangeably 
throughout  the  geotechnical  engineering  literature.  When  fabric  was  first  introduced  during 
the  mid-1960s,  it  was  used  interchangeably  with  the  terms  "soil  structure"  and 
"microstructure"  (Casagrande,  1932;  Lambe,  1953).  However,  there  is  a  significant 
difference  between  the  two  terms.  Mitchell  (1976)  describes  soil  fabric  as  "...the 
arrangement  of  particles,  particle  groups,  and  pore  space  within  a  soil."  In  addition,  he 
also  defines  soil  structure  as  "...the  combined  effects  of  fabric,  composition,  geometrical 
and  skeletal  arrangement  of  soil  particles,  and  the  interparticle  forces  acting  on  them."  In 
sands,  soil  structure  includes  gradation,  arrangement  of  particles,  and  void  ratio. 
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The  term  "fabric”  has  different  meanings  in  the  disciplines  of  geotechnical 
engineering  and  geology.  Pettijohn  (1963)  states  that  geologists  term  fabric  as  "the 
arrangement  or  orientation  (random  or  preferred)  of  the  elements  that  compose  a  rock."  In 
igneous  or  metamorphic  petrology  the  term  fabric  includes  "the  orientation  of  grains  as 
characterized  by  their  crystal  lattice,"  Verhoggen  et  al.,  (1970).  The  difference  between  the 
geotechnical  engineering  and  geology  definitions  of  fabric  is  that  engineering  fabric  studies 
are  concerned  with  the  orientations  of  grains  and  pores,  and  geological  studies  are 
concerned  with  the  orientations  of  the  crystal  lattice  which  make  up  the  rx5ck. 

Microstructure  is  the  level  of  structure  that  requires  an  optical  microscope  for  the 
study  of  the  soil  structure.  Structure  in  sands  can  be  studied  using  a  standard  optical 
microscope.  However,  in  clays,  a  Scanning  Electron  Microscope  (SEM)  must  be  used  for 
clays.  Soil  microstructure  as  defined  by  Mitchell  (1976)  is  the  term  being  adopted  for  this 
investigation  since  this  term  refers  to  the  skeletal  arrangement  (orientation)  of  soil  particles 
and  the  forces  between  adjacent  particles.  Also,  the  soil  structure  responds  to  external 
changes  in  the  environment  such  as  applied  loads  (static  or  dynamic),  moisture  (degree  of 
saturation),  and  temperature  (freeze-thaw  cycles)  (Bowles,  1984). 

Two  useful  subconcepts  exist  concerning  structure  in  a  granular  soil:  1) 
orientation  of  an  individual  particle,  and  2)  the  position  of  a  particle  and  its  mutual 
relationship  to  other  particles.  One  of  the  most  important  structural  components  is  the 
absolute  position  of  the  particle  with  respect  to  a  fixed  coordinate  system.  The  importance 
in  this  concept  is  that  the  orientation  of  an  individual  particle  and  absolute  position  of  the 
particle  are  used  to  qualitatively  describe  soil  microstructure  (Oda,  1972a). 

2.  Particle  Size  and  Shape 

A  description  of  the  shape  or  geometric  form  of  a  particle  involves  several 
separate  but  related  geometric  concepts.  The  larger  particles  in  a  soil,  namely  silt,  sand, 
and  gravel  sizes  may  be  spherical  and/or  angular  in  shape.  Angularity  and  sphericity  are 
two  parameters  used  to  quantitatively  described  the  shape  of  soil  particles.  The  angularity 
of  a  soil  particle  is  a  measure  of  the  sharpness  of  grain  surface  features.  For  example,  an 
angularity  value  of  0. 1  indicates  that  the  soil  particle  has  sharp  edges  while  a  value  of  0.9 
would  indicate  a  well-rounded  particle.  The  sphericity  of  a  soil  particle  is  a  measure  of 
how  closely  the  particle  approaches  a  spherical  shape.  For  example,  a  sphericity  value  of 
0.1  would  indicate  that  the  soil  particle  is  irregular  in  shape  and  does  not  resemble  a  sphere, 
while  a  value  of  0.9  would  indicate  that  the  particle  approaches  the  shape  of  a  sphere. 
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Sphericity  and  angularity  are  useful  in  determining  the  origin  and  geological 
history  of  the  material.  Prior  to  weathering  and  decomposition  of  the  parent  rock,  many  of 
the  soil  particles  produced  are  initially  angular  in  shape.  With  time,  mechanical  forces 
acting  on  the  soil  particles  tend  to  reduce  the  angularity,  thereby  increasing  sphericity. 
Particles  that  have  been  weathered  mechanically  over  a  long  time  will  be  much  more 
rounded  and  have  high  sphericity  and  angularity  values  (Brewer,  1964). 

Angularity  values  are  given  by  Pettijohn  (1963)  in  Table  1,  as  well  as  a  visual 
description  of  each  class  range.  After  the  sphericity  and  angularity  values  have  been 
determined,  visual  charts  such  as  Figure  1  (Krumbien  and  Sloss,  1955)  can  be  used  as  a 
check  on  the  particle  shape. 

Particle  size  is  determined  by  performing  a  standard  mechanical  sieve  analysis, 
(ASTM  D  422-90).  Particle  shape  (angularity  and  sphericity)  can  be  determined  by 
viewing  granular  soils  under  a  microscope  and  sketching  a  two-dimensional  image  of  the 
individual  grains.  Angularity,  A,  is  determined  by  constructing  the  largest  inscribed  circle 
within  the  particle  and  by  finding  the  radius  of  individual  particle  comers,  (Krumbien  and 
Sloss,  1955).  This  is  shown  in  Figure  2a  and  determined  from  Equation  1, 

(1) 

where,  q  is  the  radius  of  individual  particle  comers,  R  is  the  largest  inscribed 
radius  within  a  particle,  and  n  is  the  number  of  individual  particles. 

Sphericity,  S,  is  determined  by  constructing  a  circumscribed  circle  around  the 
entire  particle.  The  sphericity  is  a  measure  of  the  ratio  of  projected  particle  area  over  the 
area  of  the  smallest  circumscribing  circle  (Krumbien  and  Sloss,  1955).  This  is  shown  in 
Figure  2b  and  determined  fit)m  Equation  2, 

S=^  (2) 

where,  Ap  is  the  area  of  a  particle  and  A^.  is  the  area  of  the  smallest  circumscribing  circle. 

3 .  Qiaracteristics  of  Soil  Stmctural  Elements 

In  general,  soils  consist  of  three  components,  skeleton  grains,  soil  matrix,  and 
pores  (Lafeber,  1962).  Each  of  these  components  make  up  the  structure  of  the  soil  and 
differ  in  strength  and  deformation  properties.  Skeleton  grains  are  the  coarser  particles 
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TABLE  1.  Description  of  Roundness  Classes  (Pettijohn,  1957). 


Class  Name 

Roundness  Valued 

Description 

Angular 

0-0.15 

Strongly  developed  faces  with  sharp  edges 
and  comers;  secondary  comers^  are 
numerous. 

Subangular 

0.15  -  0.25 

Strongly  developed  faces  with  somewhat 
rounded  edges  and  comers;  secondary 
comers  are  numerous. 

Subrounded 

0.25  -  0.40 

The  edges  and  comers  are  rounded  and  the 
area  of  flat  faces  is  comparatively  small; 
secondary  comCTS  are  much  rounded  and 
reduced  in  number  (5  - 10). 

Rounded 

0.40  -  0.60 

Flat  faces  are  practically  absent;  all  edges  and 
comers  are  rather  broad  curves,  and  there 
may  be  broad  re-entrant  angles;  secondary 
comers  have  disappeared. 

WeU-Rounded 

0.60  -  1.00 

There  are  no  flat  faces;  the  entire  surface 
consists  of  broad  curves. 

Note:  ^Roundness  value  is  equivalent  to  axial  ratio, 

*’S^ondary  comers  are  the  many  minor  convexities  seen  in  the  grain  profile. 
Primary  comers  are  the  principal  interfacial  edges  and  are  fewer  in  niunber  (3-5). 


which  generally  consist  of  minerals  or  rock  fragments.  The  grains  are  moderately  to  highly 
resistant  to  physical  and  chemical  weathering.  The  ratio  of  surface  area  to  volume  of  the 
grains  is  relatively  small.  Individual  grains  have  high  strength  with  respect  to  the  low  to 
moderate  stresses  encountered  in  soil  engineering  practice. 

The  soil  matrix  is  defined  as  a  material  that  encloses  other  recognizable  features  or 
materials.  The  matrix  consists  of  finer-grained  materials  (silt  and  clay  size  particles)  packed 
between  the  mineral  grains.  Brewer  (1964)  states  an  upper  size  limit  of  30  micrometers  for 
soil  matrix  materials. 

The  pores  are  the  spaces  between  the  solid  soil  particles.  Within  the  soil  matrix, 
the  pores  may  be  filled  with  water,  air  or  a  combination  of  both.  Thus,  soil  can  be  viewed 
as  a  multiphase  system.  The  pore  phase  can  affect  soil  strength  and  overall  behavior. 
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Figure  1.  Sphericity  and  Roundness  Classes  (Brewer,  1964). 


Figure  2.  Particle  Shapes  (WU,  1966). 
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4 .  Soil  Structure  and  Soil  Mechanics 


Lafeber  (1966)  indicates  that  there  are  two  major  geometrical  aspects  of  the  spatial 
arrangement  of  soil  constituents:  spatial  distribution,  and  spatial  orientation.  Spatial 
distribution  describes  relative  locations  of  various  grains  in  a  continuous,  multicomponent, 
structural  pattern.  This  pattern  describes  the  position  of  grains  with  respect  to  other  grains 
directly  surrounding  them  and  spatial  orientation  describes  the  orientation  of  the  soil 
particles  within  the  soil  matrix  to  a  given  set  of  coordinates  (Lafeber,  1963). 

Australian  soils  analyzed  by  Lafeber  (1963)  showed  that  the  frequency  of  munial 
contacts  between  skeleton  grains  is  only  a  fraction  of  the  value  which  may  have  been 
expected  if  the  grains  had  been  in  a  random  arrangement.  Continuous  subfabric  or 
substructure  was  not  evident  in  any  of  the  soils  tested  by  Lafeber  (1966),  and  an  open 
arrangement  generally  occurred  which  results  in  a  pronounced  internal  geometrical 
anisotropy  of  the  soil.  These  anisooxipies  can  be  considered  to  be  an  expression  of  the 
depositional  or  the  previous  stress- strain  history  of  the  soil,  or  a  combination  of  both. 

There  are  three  problems  concerning  the  deformation  of  randomly  packed 
granular  sand  on  the  microscopic  level  (Oda,  1972b).  They  are:  (1)  the  mechanical  role  of 
sliding  and  rolling  of  grains  during  deformation,  (2)  the  mechanism  to  form  a  shear  plane 
in  the  granular  structure  (preferred  orientation  of  shear  plane),  and  (3)  the  mechanism  to 
control  the  structure  reconstruction  (increasing  anisotropy  of  fabric). 

The  initial  structure  of  a  granular  sand  will  change  with  an  increase  of  axial  strain 
in  order  to  withstand  a  continuously  increasing  stress  ratio  (strain  hardening).  Relative 
displacement  of  grain  particles  generally  occurs  and  particles  will  come  into  contact  with 
each  other.  Newly  formed  contact  points  or  surfaces  will  develop  as  existing  contacts 
gradually  decrease  (Oda,  1972b). 

Mahmood  and  Mitchell  (1974)  analyzed  soil  composed  of  granular  materials  in 
the  silt  to  fine  sand  range  and  found  that  different  microstructural  arrangements  developed, 
depending  on  how  the  soil  was  deposited.  Mechanical  properties  also  differed  with  the 
method  of  deposition  which  was  attributed  to  variations  in  particle  orientations. 
Characterization  of  the  microstructure  was  determined  by  grain  orientation  analysis  and  by 
measurement  of  the  pore  size  distribution. 

Mahmood  and  Mitchell  (1974)  also  used  a  mercury  intrusion  porosimetry  method 
to  determine  pore  size  distribution  and  found  that  grains  developed  a  strong  preferred 
orientation  in  the  horizontal  direction  when  the  soil  was  poured,  forming  pores  between 
one  and  30  micrometers  in  diameter.  When  the  soil  was  compacted  either  by  static  (slow 
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loading  with  a  flat  piston)  or  dynamic  (tamping)  methods,  near  random  grain  arrangement 
and  pore  sizes  between  0. 1  and  10  micrometers  in  diameter  resulted. 

Juang  and  Holtz  (1986)  used  mercury  intrusion  porosimetry  techniques  on  clean 
Ottawa  sand  and  a  mixture  of  Ottawa  sand  and  kaolinitic  clay.  Specimens  were  prepared 
using  two  packing  methods,  kneading  and  pluvial  compaction.  Characterization  of  the 
effects  of  compaction  and  compaction  moisture  content  were  analyzed  by  using  a  Pore  Size 
Distribution  (PSD)  index.  Distinct  pore  size  distributions  were  observed  for  the  different 
packing  methods  investigated. 

The  shear  strength  of  crushed  basalt  was  measured  by  Mahmood  and  Mitchell 
(1974).  When  sheared  normal  to  the  preferred  orientation,  the  shear  strength  was  30 
percent  higher  than  when  shearing  was  parallel  to  the  orientation  direction.  Mahmood  and 
Mitchell  (1974)  point  out  that  this  is  to  be  expected  because  shearing  across  the 
preferentially  oriented  grains  should  involve  breakage  or  reorientation  of  many  grains. 

B .  MICROSTRUCTURAL  ANALYSIS 

A  microstructural  analysis  can  be  defined  as  a  scientific  method  used  to  quantitatively 
describe  the  arrangement  of  particles  within  a  soil  (Oda,  1972a).  The  analysis  is  important 
because  relationships  may  be  developed  between  the  statistical  results  of  the  microstructural 
analysis  and  other  relevant  soil  parameters. 

Microstructural  analyses  are  typically  performed  by  determining  the  axial  ratio 
(similar  to  sphericity)  and  prefened  orientation  of  apparent  particle  long-axes.  Oda  (1972a) 
and  Campbell  (1985)  developed  methods  for  determining  axial  ratio  and  apparent  particle 
long-axis  orientations  which  are  discussed  in  detail  in  the  following  two  subsections. 
Sample  preparation  methods  have  also  been  related  to  particle  long-axis  orientation  and  to 
static  and  dynamic  mechanical  behavior  of  soils. 

1.  Experimental  Axial  Ratio  Determination 

Microstructural  characterization  of  soils  is  closely  related  to  particle  shape, 
especially  to  the  average  flamess  or  slenderness  of  individual  particles.  Axial  ratio  is  a 
shape  factor  used  to  describe  the  two-dimensional  shape  of  soil  particles.  Oda  (1972a) 
used  the  axial  ratio  method  to  measure  constituting  grain  particles.  To  determine  axial  ratio, 
the  apparent  long-axis,  L^,  and  apparent  short-axis,  L2,  are  measured  on  a  two- 
dimensional  section  through  an  individual  particle.  The  axial  ratio  for  a  group  of  particles, 
Uavg,  is  calculated  from: 
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(3) 


where,  f  is  the  number  of  particles  in  axial  ratio  calculation.  The  value  of  Uavg  ranges  from 
zero  to  one  and  may  be  considered  as  a  coefficient  indicating  the  characteristic  of  particle 
shape.  As  navg  approaches  imity,  the  sand  grains  approach  the  shape  of  a  circle.  This 
relates  closely  to  sphericity  as  described  in  SECTION  H. 

2 .  Apparent  Long- Axis  Determination 

Apparent  long-axis  determination  is  an  important  factor  in  microstructural 
analyses.  Orientation  of  the  longest  axis  is  the  main  factor  used  to  determine  if  a  preferred 
orientation  exists  within  the  soil  structure.  Particle  orientation  is  determined  by  the 
inclination  of  apparent  long  or  short-axes  of  a  grain  with  respect  to  a  fixed  reference  axis. 
Most  orientation  problems  require  the  treatment  of  more  than  one  variable  in  a  given 
angular  measurement.  Geological  orientation  data  or  vectorial  structure  is  the  spatial 
measurement  which  indicates  directional  properties  (Pincus,  1953). 

Nonspherical  particle  orientation  is  represented  by  the  inclination  of  long  and  short-axes  of 
each  grain  with  respect  to  fixed  coordinate  axes.  In  Figure  3,  grain  Gj  is  surrounded  by 
grains  G2,  G3,  G4,  G5  with  contact  points  Pj,  P2,  P3,  P4  respectfully.  Nj  is  the  normal  to 
the  tangential  plane  P^  and  can  be  defined  by  the  angle  with  the  reference  axes  X,  Y,  and  Z 
The  mutual  relation  of  particle  G^  to  particles  G2  to  G5  is  defined  by  the  number  of  contacts 
and  directions  of  Nj  to  N4  (Oda,  1972a).  Oda  (1972a)  points  out  some  problems  in 
finding  a  apparent  long  or  short-axis  of  fine  grained  particles  in  thin  sections.  For  particles 
having  an  axial  ratio  greater  than  0.7,  the  long-axis  was  difficult  to  define  because  the 
particle  approached  the  shape  of  a  sphere.  The  state  of  three-dimensional  grain  orientations 
was  estimated  on  the  basis  of  the  degree  of  preferred  orientation  of  apparent  long-axes  by 
following  these  two  steps: 

a)  Measure  the  angles  Oj  between  apparent  long-axes  of  about  200  grains  and  the 
reference  axis  X,  and 

b)  Develop  frequency  histograms  of  0,  obtained  from  measurements. 

Figure  4  shows  typical  frequency  histograms  of  Oj  for  sand  compacted  using  the  tamping 
method.  The  V-section,  45-section,  and  H-section  represent  vertical,  45-degree,  and 
horizontal  cut  sections  through  the  specimen,  respectively.  The  frequency  distributions 
show  that  the  V-section  and  45-section  have  a  statistically  preferred  orientation  of  apparent 
long-axes,  whereas  the  H-section  does  not  have  a  preferred  orientation.  This  is  seen  in  the 
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Figure  3.  Fabric  Elements  of  Granular  Soil  (Oda,  1972a). 
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Figure  4.  Frequency  Histograms  of  in  Various  Sections  of  Sand  A  (Oda,  1972  a), 


frequency  values  for  each  9;  interval  in  which  a  greater  number  of  orientations  were 
observed  from  -30  to  +30  degrees. 

Sand  A  is  a  poorly  graded  sand  and  had  an  axial  ratio  of  0.605.  The  sand  was 
compacted  by  the  tamping  method,  and  was  characterized  by  the  arrangement  of  apparent 
long  and  short-axes  of  particles  parallel  to  the  horizontal  plane  (X  axis)  and  frequency 
histograms  are  shown  in  Figure  4.  This  sand  is  composed  of  nonspherical  particles.  Oda 
found  that  sands  composed  of  relatively  spherical  particles  have  an  orientation  which  is 
almost  isotropic.  This  means  that  particles  which  are  relatively  spherical  will  tend  to  show 
the  same  orientation  regardless  of  the  section.  The  Rayleigh  test  of  significance  can  be 
used  to  determine  whether  or  not  a  preferred  orientation  is  present  in  the  soil  structure 
(Curray,  1956).  The  test  is  widely  used  for  particle  long-axis  orientation  data. 

3 .  Particle  Orientation  and  Mechanical  Behavior  of  Soils 

Microstructural  analyses  of  cohesionless  soils  can  be  used  to  relate  particle 
orientation  to  the  mechanical  behavior.  Particle  orientation  is  characterized  by  investigating 
apparent  long-axes  of  thin  sections.  The  method  has  sufficient  accuracy  for  three- 
dimensional  mineral  grains  (Mahmood  and  Mitchell,  1974).  When  specimen  sections  are 
taken,  a  two-dimensional  shape  is  obtained,  which  is  generally  ellipsoidal;  therefore,  the 
long-axes  can  be  estimated  with  sufficient  accuracy. 

Oda  (1972a,  1972b)  observed  that  long-axis  orientations  in  sands  are  a  function 
of  particle  shape  and  method  of  compaction.  He  also  noted  that  the  initial  structure  is 
important  to  the  mechanical  properties  of  the  soil.  Oda's  (1972a)  research  consisted  of 
fixing  the  structure  of  sand  with  polyester  resin  in  order  that  the  spatial  arrangement  would 
not  be  disturbed.  From  viewing  thin  sections  under  an  optical  microscope  to  find  preferred 
orientations  and  spatial  relations,  Oda  (1972a)  observed:  (1)  the  characteristics  of 
configuration  relations  of  nonspherical  particles  are  determined  by  the  shape  of  particles 
and  by  the  method  of  compaction,  (i.e.,  sands  which  are  composed  of  flat  or  elongated 
particles  indicate  strong  anisotropic  features  of  structure  character),  and  (2)  the  initial 
structure  of  sand  has  important  influences  on  the  mechanical  properties,  such  as  mobilized 
stress  ratio  and  secant  deformation  modulus. 

Oda  (1972a)  has  shown  that  "...when  sand  grains  are  deposited  in  a  gravity  field 
irrespective  of  their  shape  character,  they  tend  to  rest  in  a  stable  position  relative  to  forces 
acting  upon  them."  An  important  finding  is  that  the  nature  and  degree  of  anisotropy  of 
sand  is  determined  by  the  way  in  which  the  sand  is  deposited  and  by  the  shape 
characteristics  of  the  sand  grains.  The  importance  of  this  concept  comes  from  the  type  of 
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soil  used  for  an  experiment.  A  soil  which  is  well  rounded  generally  will  not  show  a 
preferred  orientation.  However,  angular  elliptical  soil  grains  generally  shows  a  preferred 
orientation.  Therefore,  the  type  of  soil  used  in  an  investigation  (with  regards  to  the  shape 
characteristics)  can  greatly  affect  the  possibility  of  prefened  orientation  within  the  soil. 

Mitchell  (1976)  studied  the  relationship  between  soil  structure,  sample  preparation 
methods,  liquefaction,  and  drained  compression  behavior  under  triaxial  loading  for 
Monterey  No.  0  Sand.  Three  methods  of  compaction  were  used:  pluviation,  moist 
tamping,  and  moist  vibration.  Characterization  of  microstructure  was  based  on  particle 
long-axis  and  interparticle  orientations.  For  each  method  studied,  a  preferred  orientation  of 
long-axes  developed  in  the  horizontal  direction.  The  degree  of  preferred  orientation 
increased  in  the  order  of  moist  vibration  (nearly  random),  moist  tamping,  and  dry 
pluviation.  Two  quantitative  structural  elements  were  studied;  three-dimensional 
distributions  of  normals  to  interparticle  contact  planes,  and  the  orientations  of  apparent 
long-axes  on  a  two-dimensional  plane  from  a  thin  section.  Mitchell  (1976)  concluded  that 
moist  vibration  gave  the  most  random  orientations  of  particle  long-axes.  Dry  pluviation 
had  a  distinct  preferred  orientation  in  the  horizontal  direction,  and  moist  tamping  gave 
intermediate  values.  Higher  cyclic  triaxial  strengths  were  associated  with  the  most  random 
particle  orientations. 

Particle  orientation  diagrams  drawn  from  three  different  compacting  methods  are 
shown  in  Figure  5.  Mahmood  and  Mitchell  (1974)  found  that  when  the  grains  were  in  the 
loosest  state  (pouring),  a  strong  preferred  orientation  existed  in  the  horizontal  direction.  In 
static  compression,  the  orientatitMi  was  nearly  45  degrees  from  the  horizontal  plane.  When 
tapping  the  sample  holder,  i.e.,  dynamic  compaction,  a  near  random  grain  orientation  was 
produced.  The  shear  strength  was  observed  to  be  greater  when  samples  are  sheared  across 
the  preferred  orientation  direction  than  when  sheared  parallel  to  the  preferred  orientation. 

Ladd  (1977)  studied  the  cyclic  structural  stability  of  two  sands  using  the  cyclic 
triaxial  strength  test  The  sands  were  a  fine  sand  with  a  trace  of  silt  and  a  coarse  to  fine 
sand  with  a  trace  of  silt  Four  different  methods  of  sample  preparation  were  used:  dry 
vibration,  moist  vibration,  dry  tamping,  and  moist  tamping.  He  found  that  in  both  sands, 
samples  prepared  using  dry  methods  had  similar  cyclic  shear  strengths.  However,  similar 
but  greater  cyclic  shear  strengths  occurred  using  moist  preparation  methods.  DeGregorio 
(1990)  used  Ottawa  F-70  banding  sand  in  a  study  of  liquefaction  behavior.  Samples  were 
tested  using  a  dead-load  apparatus  and  a  hydraulic  loading  system.  DeGregorio  observed 
that  soil  strength  was  influenced  by  sample  preparation  method  and  not  by  the  type  of 
loading  system.  Soil  strength  increased  in  order  of  dry  pluviation,  moist  tamping,  and 
moist  vibration  when  sheared  with  the  dead-load  device.  The  effect  of  compaction  method 
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Figure  5.  Particle  Orientation  Diagrams  for  Crushed  Basalt: 
(a)  pouring;  (b)  dynamic  compaction;  (c)  static 
compaction  (Mahmood  and  Mitchell,  1S^4). 


Figure  6.  Technique  for  Marking  Apparent  Long-  and  Short- 
Axes  of  Sand  Particles  in  Thin  Sections  (Campbell, 
1985). 
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on  soil  behavior  was  demonstrated  in  each  study.  Grain  orientations  may  have 
significantly  influenced  the  results  of  these  two  studies,  however,  no  microstructural 
analyses  were  performed. 

The  survey  of  literature,  makes  it  apparent  that  structure  in  granular  materials 
depends  on  the  method  of  sample  preparation  and  compaction.  When  soil  samples  are 
densified  either  by  dynamic  (tamping,  plunging,  impact,  vibration)  or  static  compaction, 
internal  stresses  are  developed  and  grains  do  not  maintain  the  iiutial  structure  as  formed 
when  the  soil  is  initially  deposited.  However,  structures  developed  using  pluviation  remain 
unchanged,  since  the  soil  is  compacted  during  placement. 

C .  PARTICLE  ORIENTATION  AND  STATISTICAL  ANALYSIS 

Orientations  of  particle  long-axes  in  sands  are  determined  by  using  microstructural 
analysis  techniques.  Microstructural  analysis  is  an  important  tool  used  to  view  particles  or 
voids  in  a  soil  which  have  been  subjected  to  various  preparation  methods  and  loading 
conditions. 

1 .  Specimen  Preparation  Techniques  for  Microstructural  Analysis 

Several  sample  preparation  methods  for  microstructural  analysis  have  evolved 
over  the  past  two  decades.  For  successful  microstructural  analysis  of  granular  soils,  the 
soil  structure  must  be  preserved  so  that  it  can  withstand  sectioning  procedures.  Some  of 
the  sample  preparation  methods  used  by  various  researchers  are  given  in  the  following 
paragraphs. 

Oda  (1972a)  prepared  specimens  by  pouring  oven-dried  soil  into  a  cylindrical 
mold.  Compaction  was  performed  using  a  tapping  method  in  which  the  side  of  the  mold 
was  tapped  with  a  tamper,  and  a  plunging  method  in  which  the  hand  tamper  was  plunged 
directly  into  the  sand.  A  polyester  resin  was  used  to  fix  the  soil  structure  in  the  original 
configuration.  Thin  sections  were  prepared  and  photographed  for  microstructural  analysis. 

Mitchell  (1976)  prepared  specimens  of  Monterey  No.  0  sand  using  three  different 
compaction  methods,  pluviation,  horizontal  high  frequency  vibrations,  and  moist  tamping. 
A  polyester  resin  was  mixed  with  styrene  (20  percent  by  volume)  and  1  percent  organic 
peroxide  as  a  catalyst  to  reduce  viscosity  was  used  to  preserve  the  soil  structure.  Thin 
sections  (60  mm  thick  and  20  mm  wide)  were  prepared  for  microstructural  analysis. 

Rein  (1983)  prepared  samples  for  microstructural  analysis  by  cutting  ice-rich 
frozen  sands  with  a  wire  saw.  The  research  involved  a  study  of  the  mechanical  behavior  of 
frozen  soils.  Samples  were  prepared  by  compacting  a  sand/snow  mixture  in  plastic  tubes. 
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A  vacuum  was  applied  to  saturated  the  sample  at  0“  C  and  sections  were  taken  with  a  wire 
saw  at  -10®  C. 

Campbell  (1985)  used  a  sugar  cementation  technique  to  preserve  the  soil 
structure.  His  procedure  involved  injecting  3(X)  grams  of  granulated  sugar  dissolved  in  1 
liter  of  water  by  drawing  the  solution  through  the  specimen  under  a  low  gradient  until  full 
saturation  was  achieved.  Samples  were  dried  by  drawing  a  flow  of  dry  air  through  the 
specimen  under  a  low  gradient  Thin  sections  were  prepared  by  filling  pore  spaces  with  a 
low  viscosity  resin,  Petropoxy  154. 

Gill  et  al.  (1986)  studied  the  relationship  between  microscale  physics  and 
macroscale  response  of  bonded  soil  particles.  Prepared  specimens  were  tested  under 
uniaxial  strain  conditions  at  varying  amounts  of  macroscopic  strain.  Once  unloaded,  the 
lightly  cemented  granular  specimens  were  epoxied.  Thin  sections  then  were  taken, 
permitting  forensic  viewing  of  the  internal  microstmctural  response. 

2 .  Two-Dimensional  Analysis  of  Orientation  Data 

Analyzing  two-dimensional  orientation  data  requires  the  use  of  a  vector 
method  and  a  statistical  test  of  significance.  The  vector  method  is  mathematically  based  on 
the  same  idea  as  vectors  in  engineering  statics.  The  only  difference  is  that  the  vector 
magnitude  is  taken  as  unity  and  the  vector  is  bidirectional.  The  statistical  test  of 
significance  is  based  on  the  work  of  Rayleigh  (1894)  in  which  a  method  for  describing 
random  phases  in  sound  waves  was  developed. 

a.  Vecttjr  Method 

Curray  (1956)  notes  that  measurements  of  orientation  of  various 
geological  features  such  as  fabric  and  structural  elements  can  be  more  easily  interpreted 
through  the  use  of  statistics.  Three  relevant  statistical  measurements  can  be  made:  (1)  a 
preferred  orientation  direction,  (2)  the  degree  of  preferred  orientation,  and  (3)  the 
probability  that  the  preferred  orientation  is  real  and  not  merely  due  to  chance. 

In  the  past,  researchers  used  standard  linear  normal  distributions  and 
computed  mean  and  standard  deviation  (or  variance)  as  descriptive  statistics  (Curray, 
1956).  Jizba  (1953)  and  Chayes  (1954)  noted  difficulties  which  arise,  such  as  choosing  an 
origin  in  order  to  divide  a  circular  distribution  into  a  linear  frequency  curve.  By  changing 
the  choice  of  origin,  the  calculated  mean  and  standard  deviation  can  be  considerably 
different. 
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The  best  method  available  for  studying  particle  long-axis  orientation  is  the 
vector  method  since  the  statistics  used  are  independent  of  the  origin  chosen.  Curray  (1956) 
lists  four  requirements  for  analyzing  particle  orientation  data:  (1)  a  measme  of  the  central 
tendency  or  preferred  orientation  direction,  (2)  a  measure  of  dispersion  which  is 
independent  of  origin,  (3)  a  test  of  the  statistical  significance  of  the  results  against  a  model 
randomness,  and  (4)  a  model  distribution  such  as  the  circle  normal  or  the  linear  normal 
distributions  from  which  deviations  can  be  tested. 

Campbell  (1985)  developed  a  rapid  method  of  directly  measuring 
microstructure  in  cohesionless  soils.  A  personal  computer  and  digitizing  tablet  were  used 
to  develop  a  relationship  between  long-axis  coordination  and  contact  normal  orientation 
data.  Campbell  performed  analyses  of  thin  sections  photographs  which  involved  the 
measurement  of  particle  orientation  information  from  particle  shape,  and  orientation  of  the 
particle  with  respect  to  a  fixed  axis  system.  For  the  analyses,  approximately  250  particles 
were  selected  for  statistical  review.  Particle  apparent  long  and  short-axes  were  marked  by 
circumscribing  a  rectangle  around  each  particle.  This  technique  permitted  visualization  of 
where  the  long  and  short-axes  of  each  particle  occur  and  is  shown  schematically  in  Figure 
6. 

A  typical  histogram  of  apparent  long-axis  orientation  within  10-intervals  is 
shown  in  Figure  7.  A  rose  diagram  (angular  frequency  histogram)  developed  from  the 
same  data  is  shown  in  Figure  8.  The  rose  diagram  is  a  plot  of  the  azimuth,  0,  of  each 
particle  long-axis  on  a  360-degree  distribution  graph.  Figures  7  and  8  visually  present  the 
orientation  of  azimuths  and  the  dashed  line  indicates  a  uniform  orientation  distribution.  The 
azimuths  of  a  group  of  particles  (determined  from  the  long-axes)  can  be  detected. 
Equations  (4)  to  (10)  are  used  to  analyze  these  data.  The  technique  was  presented  by 
Curray  (1956). 

In  particle  orientation  analysis,  no  distinction  is  made  between  one  end  of 
a  particle  and  the  other,  i.e.,  the  particles  are  bidirectional.  Therefore,  all  measurements  are 
made  in  a  range  of  180-degrees.  A  360-degree  distribution  graph  may  be  generated  by 
plotting  each  grain  twice,  once  from  0-  to  180-degrees,  and  again  in  the  180-  to  360-degree 
range.  Curray  (1956)  presents  Equations  (4)  to  (10)  which  take  this  into  account  by 
doubling  the  azimuth,  9,  of  each  particle  long-axis: 

Horizontal  Component  =  Zn'  cos29  (4) 

Vertical  Component  =  Zn'  sin29  (5) 
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Figure  7.  Example  of  Histogram  Showing  Apparent  Long- 
Axis  Orientation  Data  (Campbell,  1985) 


39 

I 


Figure  8.  Example  of  Rose  Diagram  Showing  Apparent-Long 
Axis  Oreintation  Data  (Campbell,  1985) 
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where  n'  is  the  magnitude  of  each  particle  long-axis  vector  (taken  as  unity  in  this  case),  0  is 
the  azimuth  of  each  particle  long-axis,  0'  is  the  azimuth  of  resultant  vector,  r  is  the 
magnitude  of  resultant  vector,  and  L  is  the  magnitude  of  resultant  vector  in  terms  of  a 
percentage. 

Mitchell  (1976)  indicates  that  particle  orientation  can  be  represented  by  its 
apparent  long  or  short-axis  with  respect  to  a  fixed  reference  axis.  True  orientations  of  long 
and  short-axes  of  small  particles  are  difficult  to  determine.  Three-dimensional  orientations 
are  estimated  from  orientations  observed  in  sections  cut  in  different  directions.  Orientations 
were  determined  by  measuring  the  angle  %  between  apparent  long-axes  and  reference  axes 
parallel  to  the  long-axis  of  the  thin  section  shown  in  Figure  9.  Typical  rose  diagrams  are 
shown  in  Figure  10  for  data  obtained  using  an  azimuth  class  interval  of  0  =  10-degrees. 

b .  Statistical  Test  of  Significance 

A  method  developed  by  Rayleigh  (1894)  can  be  used  to  test  the 
significance  of  a  distribution  of  particle  long-axis  orientations  (Campbell,  1985).  The  term 
"significance"  implies  that  the  observed  data  do  or  do  not  show  a  preferred  orientation. 
The  Rayleigh  test  of  significance  has  been  widely  used  in  particle  orientation  analyses  and 
is  analogous  to  particle  long-axis  orientation  data.  This  can  be  explained  by  knowing  that 
each  observation  (particle  long-axis)  can  be  represented  by  a  vector  with  direction  and 
magnitude.  The  Rayleigh  test  of  significance  is  stated  as: 

p  =  e-'^/"  (10) 

where  p  is  the  probability  of  obtaining  a  greater  amplitude  by  pure  chance  combination  of 
random  phases  (for  particle  orientation  analysis,  it  is  the  probability  of  obtaining  a  greater 
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Figure  9.  Thin  Section  Locations  from  Within  Resin- 
Impregnated  Samples  (Mitchell,  1974). 
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Figure  10.  Rose  Diagrams  of  Particle  Long-Axis  Orientations  for 
Samples  of  Monterey  No.  0  Sand  Prepared  at  50% 
Relative  Density  by  Different  Methods  (Mitchell, 
1974). 
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magnitude  by  pure  chance  combinations  of  random  orientations),  r  is  the  resultant 
amplitude  (for  particle  orientation  analysis,  it  is  the  magnitude  of  the  resultant  vector),  n  is 
the  number  of  components  or  phases,  and  e  is  the  base  of  natural  logarithms. 

A  statistical  analysis  was  performed  by  Campbell  (1985)  using  the  vector 
method  and  test  of  significance  presented  by  Curray  (1956),  using  Equations  (4)  to  (10). 
For  sand  grain  orientation  studies,  a  Rayleigh  significance  level  of  0.05  is  required.  If 
p>0.05,  additional  observations  should  be  made  until  p<0.05.  For  long-axis  orientation 
data,  a  minimum  sample  size  of  about  200  observations  should  provide  an  accurate 
measurement. 

The  vector  mean  and  vector  magnitude,  must  be  determined  to  ensure  a 
suitable  estimation  of  the  mean  direction  and  variance  of  the  frequency  histograms.  Vector 
mean  is  a  measure  of  the  average  azimuth  upon  which  a  group  of  particles  lie.  Vector 
magnitude  is  a  measure  of  how  closely  the  particle  long-axis  distribution  approaches  a 
preferred  orientation.  If  the  vector  magnitude  is  zero,  then  a  completely  random 
distribution  exists.  If  the  vector  magnitude  is  unity,  then  all  observations  of  particle  long- 
axes  have  the  same  azimuth. 

D.  SOIL  COMPACnON 

Compaction  methods  are  used  in  almost  every  phase  of  engineering  construction  to 
improve  soil  conditions  on  site.  The  basis  of  compaction  theory  was  developed  by  Proctor 
in  the  late  1920s  and  a  standard  compaction  test  was  developed  in  the  early  1930s. 
Compaction  energy  depends  on  soil  type  and  amount  of  moisture  within  the  soil. 
Cohesionless  soils  differ  fiom  cohesive  soils  in  their  behavior  with  respect  to  moisture- 
density  relationships.  These  topics  are  described  in  detail  in  the  following  three  sections. 

1 .  Principles  and  Theory  of  Soil  Compaction 

Compaction  theory  first  came  into  existence  through  the  work  of  R.  R.  Proctor 
who  developed  a  controlled  specification  for  the  compaction  of  cohesive  soils  during  the 
late  1920s,  Proctor  showed  that  soil  compaction  is  a  function  of  the  following  four 
variables:  (1)  compaction  energy,  (2)  soil  type  (gradation,  particle  size  and  shape, 
cohesionless  or  cohesive  soil),  (3)  water  content,  Wc,  and  (4)  dry  unit  weight,  The 
method  of  soil  compaction  and  moisture  condition  during  compaction  strongly  influences 
both  static  and  dynamic  soil  behavior.  Soils  which  are  in  a  loose  state  tend  to  have  lower 
compression  and  shear  wave  speeds  and  shear  strength,  while  densely  compacted  soils 
have  higher  compression  and  shear  wave  speeds  and  shear  strength. 
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Holtz  and  Kovacs  (1981)  define  compaction  as  "...the  densiflcation  of  a  soil  by 
the  application  of  mechanical  energy."  Cohesionless  and  cohesive  soils  differ  in  their 
behavior  when  subjected  to  compaction  energy.  Cohesionless  soils  can  be  efficiently 
compacted  in  the  field  using  vibratory  techniques  or  by  using  dynamic  compaction  in  which 
large  weights,  10  to  40  tons,  are  dropped  from  heights  of  15  to  45  meters  (Hunt,  1986). 
Cohesive  soils  may  be  compacted  in  the  field  using  sheepsfoot  rollers  where  a  kneading 
action  occurs,  statically  by  applying  large  static  loads,  or  compacted  in  the  lab  by  a  falling 
hammer. 

The  influence  of  compaction  method  on  engineering  behavior  has  been 
demonstrated  by  a  number  of  researchers.  The  underlying  reason  for  compacting  soils  is  to 
improve  the  engineering  properties.  Holtz  and  Kovacs  (1981)  note  several  advantages  of 
soil  compaction:  (1)  detrimental  settlements  can  be  reduced  or  prevented,  (2)  soil  strength 
increases  and  slope  stability  can  be  improved,  (3)  bearing  capacity  of  pavement  structures 
can  be  improved,  and  (4)  undesirable  volume  changes,  caused  by  frost  action,  swelling, 
and  shrinkage  for  example,  may  be  controlled. 

Mulilis  et  al.  (1977)  prepared  samples  of  fine  sands  for  cyclic  liquefaction  tests 
using  1 1  different  compaction  methods.  Results  showed  that  compaction  method  strongly 
influenced  the  cyclic  liquefaction  behavior  of  the  sand.  Of  the  1 1  different  packing 
methods,  pluviated  samples  had  the  greatest  liquefaction  potential  whereas  samples 
prepared  by  vibration  and  hand  tamping  had  the  lowest  liquefaction  potential. 

Oda  (1972a)  prepared  specimens  by  pouring  oven-dried  soil  into  a  cylindrical 
mold.  Specimens  were  compacted  using  a  tapping  method  in  which  the  side  of  the  mold 
was  tapped  with  a  tamper,  and  a  plunging  method  in  which  a  hand  tamper  was  plunged 
directly  into  the  sand.  Oda  found  that  preferred  particle  orientation  is  influenced  by  the 
method  of  compaction,  observing  that  pluviated  samples  have  a  definite  particle  orientation 
and  lower  shear  strength,  whereas  taiiq)ed  samples  showed  a  more  random  orientation  with 
a  higher  shear  strength. 

A  relationship  between  compaction  and  shear  strength  was  observed  by  Olsen 
(1963)  who  studied  soils  containing  small  amounts  of  sands.  Olsen  reviewed  shear 
strength  theories  and  found  that  effective  stress  thecnies  of  soil  compaction  are  based  on  the 
assumption  that  soil  resists  the  compaction  pressure  by  the  development  of  shearing 
stresses  on  the  contact  surfaces  between  the  particles.  Prior  to  compaction,  soil  possesses 
low  shear  strength  and  an  increase  in  the  soil  density  can  be  explained  by  soil  behavior 
during  compaction  (Olsen,  1963).  As  compaction  energy  is  applied,  shearing  stresses 
developed  between  the  soil  particles  reach  the  shearing  strength  of  the  contact  surfaces,  the 
contacts  yield,  particles  slide  over  each  other,  and  soil  density  increases.  As  the  first  blow 
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is  applied,  high  pore  water  pressure  weakens  the  soil  by  reducing  the  effective  stress 
between  particles.  As  more  blows  are  applied  to  a  soil,  effective  stress  increases  due  to  an 
increase  in  residual  lateral  total  stress  and  pore  water  pressure  becomes  progressively 
smaller. 

Soil  strength  is  governed  by  effective  stress.  In  saturated  soils,  the  effective 
stress  can  be  calculated  with  reasonable  accuracy  using  Terzaghi’s  equation: 

a'  =  CT-u^  (11) 

where,  &  is  the  effective  stress,  a  is  the  total  stress,  and  u^  is  the  pore  water  pressure. 
When  a  soil  is  unsaturated,  a  third  phase,  air,  will  be  present  with  a  pressure  u^.  Olsen 
(1963)  notes  that  Bishop  (1961)  suggested  the  effective  stress  in  unsaturated  soils  could  be 
expressed  by: 

a’  =  CT-u^(%)-u^(l-X)  (12) 

where,  X  is  the  saturation  parameter  and  u,  is  the  pore  air  pressure.  The  value  of  X  varies 
depending  on  the  degree  of  saturation.  For  example,  if  the  soil  is  saturated  then  X  is  1.0 
and  if  the  soil  is  dry  then  X  is  zero.  For  unsaturated  soils  the  value  of  X  is  in  the  range 
between  zero  and  1.0.  Bishop  used  Equation  (12)  to  analyze  the  shear  strength  of  soils 
with  high  degrees  of  saturation  (usually  over  60  percent)  and  found  that  the  shear  strength 
decreased  with  increasing  saturation. 

2 .  Standard  Proctor  Test 

The  standard  Proctor  test  (ASTM  D  698-90)  was  developed  by  Proctor  in  the  late 
1920's  as  a  control  specification  for  the  compaction  of  cohesive  soils.  The  research  was 
prompted  by  dam  construction  for  the  Los  Angeles  Water  District.  The  present  standard 
Proctor  test,  designated  as  ASTM  D  698-90,  involves  dropping  a  25.4  N  (5.5  pound) 
weight  from  a  height  of  0.305  m  (12.0  inches)  onto  the  soil  contained  in  a  compaction 
mold  of  volume  0.944  m^  (1/30  ft^).  Three  layers  are  used  with  25  blows  applied  to  each 
layer.  The  test  is  repeated  several  times  using  increasing  water  contents.  A  relationship 
between  dry  density.  Yd,  and  water  content,  Wg,  can  be  developed  from  the  data.  The  total 
compactive  energy,  CE,  applied  to  a  soil  using  the  standard  Proctor  test  is  592.7  kJ/m^ 
(12,375  ft-lbffft^).  This  method  is  considered  a  dynamic  impact  compaction  test. 
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3 .  Compaction  Energy  in  Unsaturated  Soils 

Relationships  between  dry  density  and  moisture  content  can  be  developed  for 
dynamically  compacted  cohesionless  soils.  The  curve  in  Figure  1 1  shows  that  as  the  water 
content  increases  the  dry  density  initially  decreases.  With  a  further  increase  in  water 
content  the  density  increases  until  a  higher  saturation  level  is  reached.  Lambe  and  Whitman 
(1969)  citing  Foster  (1962),  state  that,  "...this  type  of  compaction  curve  occurs  when  the 
material  is  permeable  enough  to  impede  the  development  of  positive  pore  pressures  during 
compaction."  Lambe  and  Whitman  (1969)  attribute  low  densities  at  lower  water  contents  to 
bulking,  a  phenomena  which  occurs  when  the  capillary  forces  of  the  soil  resist  particle 
rearrangement  At  lower  water  contents  (lower  degrees  of  saturation),  capillary  tension  in 
the  pores  inhibits  the  tendency  of  the  cohesionless  soil  particles  to  be  densely  compacted. 
As  the  water  content  increases  the  soil  particles  are  lubricated,  and  thus  frictional  resistance 
between  particles  decreases  allowing  the  soil  to  be  ccnnpressed  into  a  denser  packing. 

In  cohesive  soils,  the  moisture-density  curve  starts  at  a  low  value  of  water  content  and  dry 
density.  The  curve  then  rises  to  some  peak  value  (optimum  moisture  content)  after  which 
the  dry  density  decreases  with  an  increase  in  water  content  as  illustrated  in  Figure  12.  A 
change  in  soil  structure  occurs  as  the  moisture  content  increases.  The  soil  structure  below 
optimum  moisture  content  tends  to  be  flocculated,  as  the  moisture  content  increases  the 
structure  of  the  cohesive  soil  becomes  dispersed-  The  structure  of  cohesive  soil  changes  as 
the  compactive  effort  increases  and  maximum  dry  density  increases  as  optimum  moisture 
content  decreases.  The  difference  between  moisture-density  curves  in  sand  and  clays  can 
be  attributed  to  the  grain  size  distribution,  shape  of  the  soil  grains,  specific  gravity  of  soil 
solids,  and  the  amount  and  type  of  clay  minerals  present  (Das,  1985). 

As  the  number  of  blows  applied  to  a  soil  increases,  i.e.,  an  increase  in 
compaction  energy,  the  dry  density  and  shearing  strength  of  the  soil  also  increase.  At  some 
point,  the  soil  may  become  sufficiently  strong  such  that  no  interparticle  contacts  yield  when 
a  blow  is  applied.  In  such  a  case,  the  soil  is  said  to  have  reached  its  maximum  density  for 
the  particular  compaction  method  used  (Olsen,  1963). 

Olsen  (1963)  noted  that  prior  to  the  second  blow  the  soil  is  in  a  much  denser 
condition  than  it  was  at  the  beginning  of  compaction.  The  denser  soil  produces  more  lateral 
confinement  Shearing  deformations  break  down  some  of  the  interparticle  contacts,  expel 
pore  air,  and  increase  the  dry  density.  The  compaction  effort  is  less  for  the  second  blow 
than  it  was  for  the  first  and  the  increase  in  dry  density  is  smaller,  this  is  caused  by  the 
greater  lateral  confinement  At  low  water  contents  the  dry  density  increases  quickly  for  the 
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Water  Cootent 


Figure  11.  Compaction  Curve  for  Cohesionless  Soils 
(Lambe  and  Whitman,  1969). 


Figure  12.  Compaction  Curve  for  Cohesive  Soils  (Das,  1985). 
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first  few  blows.  However,  further  density  increases  with  additional  compactive  effort  are 
small  as  shown  graphically  in  Figure  13. 

Ross  (1986)  and  Veyera  and  Fitzpatrick  (1990a)  compacted  unsaturated  Ottawa 
20-30  sand  to  a  constant  dry  density  using  a  standard  Proctor  hammer.  Their  findings 
show  that  the  amount  of  energy  required  to  compact  the  soil  to  a  constant  dry  density 
changes  with  saturation.  At  saturations  between  20  percent  and  60  percent  the  amount  of 
energy  required  to  compact  the  soU  to  a  constant  dry  density  was  much  higher  than  at  zero 
and  80  percent  saturation.  In  addition,  Ross  (1989)  also  performed  similar  tests  using 
Eglin  sand.  Figure  14  shows  the  relationship  between  saturation  and  energy  required  to 
compact  the  Eglin  soil  to  a  constant  dry  density.  Ross  found  that  a  22  percent  increase  in 
the  amount  of  energy  was  required  to  compact  the  soil  up  to  40  percent  saturation.  From 
this  point  on  the  amount  of  energy  required  to  compact  the  soil  decreased.  Ross  (1986) 
states  "It  appears  that  at  intermediate  saturations,  the  capillary  stress  between  the  sand 
particles  increases  the  stiffness  of  the  sand." 

E.  WAVE  PROPAGATION  THEORY 

Compression  and  shear  waves  are  generated  when  a  disturbance  is  introduced  into  a 
medium.  Earthquakes  produce  motions  in  the  ground  due  to  stress  waves  (compression 
and  shear)  created  by  the  rupture  of  the  stressed  earth  mass.  Compressional  waves  can  be 
generated  by  dropping  a  rock  being  into  a  still  body  of  water.  The  disturbance  from  the 
rock  will  create  compressional  waves  which  travel  outward  in  all  directions  from  the  source 
(shear  waves  are  not  developed  in  water).  Waves  also  may  be  generated  by  artificial  means 
such  as  blasting,  bombardment  during  war,  pile  driving,  and  aircraft  landing  (Prakash, 
1981).  The  following  two  sections  deal  with  wave  propagation  theory  and  specifically 
compressive  wave  velocity  in  unsaturated  soils. 

1.  Wave  Propagation  Principles 

When  a  dynamic  load  is  applied  to  a  body,  the  entire  body  is  not  disturbed  at  the 
time  of  initial  loading.  Regions  closest  to  the  disturbance  are  affected  first,  and 
deformations  produced  by  the  disturbance  spread  throughout  the  body  in  the  form  of  stress 
waves.  Wave  propagation  through  a  single  material  (such  as  a  rod)  can  be  determined  from 
Newton’s  Second  Law  of  motion  and  elasticity  theory  if  the  applied  stresses  are  within  the 
elastic  limits  of  the  material.  Equation  2.13  was  presented  by  Kolsky  (1963)  which 
determines  the  compressive  wave  propagation  velocity,  ,  in  a  rod  for  an  elastic  material: 
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Number  of  Blows  per  Layer  (N) 


Figure  13.  Effect  of  the  Number  of  Stress  Applications  oi 
Compacted  Dry  Density  of  a  Silty  Clay  with  an  Optimum 
Moisture  Content  of  12%  (Olsen,  1963). 
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where,  E  is  the  elastic  Young  s  modulus  and  p  is  the  mass  density.  Equation  13  gives  the 
compressional  wave  speed  in  a  bounded  elastic  material. 

Wave  propagation  through  several  materials  can  also  be  determined  from  wave 
propagation  theory.  When  propagating  compressional  waves  reach  a  boundary  between 
different  materials,  part  of  the  wave  will  be  reflected  back  into  the  first  medium,  and  a  part 
of  the  wave  will  be  transmitted  into  the  second  medium.  If  the  boundary  is  normal  to  the 
propagation  direction,  then  only  compressional  and  tensile  waves  will  be  reflected  and 
transmitted.  If  the  boundary  is  not  normal  to  the  propagation  direction,  then 
compressional,  tensile  and  shear  waves  will  be  reflected  and  transmitted.  Only  normal 
boundaries  will  be  considered  here.  Figure  15  shows  the  case  where  three  materials  are  in 
contact  with  one  another.  The  stress  at  each  interface  can  be  determined  using  Equations 
(14)  to  (17).  These  four  equations  were  given  by  Rinehart  (1975)  to  determine  reflected 
and  transmitted  stresses  between  different  materials. 

From  Figure  15,  assuming  normal  boundaries,  the  amount  of  stress  reflected 
back  into  the  first  material,  ,  can  be  determined  by: 


(  P2^C2  '  Pl^C  1  \ 
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and  the  amount  of  stress  being  transmitted,  Cj  ,  into  the  second  medium  is  determined  by: 
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where  V  is  the  velocity,  and  p  is  the  mass  density.  The  subscripts  T,  R,  and  I  refer  to 
transmitted,  reflected  and  incident  stresses,  respectively,  and  subscripts  1  and  2  refer  to  the 
first  and  second  mediums,  respectively.  The  product,  pV,  is  termed  the  "acoustic 
impedance"  of  the  material. 

The  magnitude  of  stress  transmitted  into  a  third  medium,  ,  can  be  determined 
by  replacing  the  incident  stress  with  the  stress  transmitted  into  medium  2  and  by  changing 
the  values  of  the  acoustic  impedances  in  Equation  (15): 
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(16) 


cs'  =(  P3^C3 

^  P2^C2  P3^C3 


If  the  first  and  third  tnediums  are  the  same  (subscripts  1  and  3),  but  differ  from 
the  second  medium,  a  relationship  between  the  incident  pulse  and  the  pulse  transmitted  into 
the  third  medium  can  be  determined.  This  particular  case  directly  applies  to  the  split- 
Hopkinson  pressure  bar  where  a  soil  sample  (medium  2)  is  placed  between  two  steel  bars 
(mediums  1  and  3).  By  combining  Equations  (15)  and  (16)  for  this  special  condition,  the 
following  equation  is  obtained; 


(  P2^C2  Pi^ci  \ 
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Equation  17  is  referred  to  as  the  "transmission  ratio."  By  comparing  Equations 
(15)  and  (17),  the  magnitude  of  the  stress  transmitted  into  the  third  medium  is 
approximately  twice  that  in  the  second  medium  if  the  acoustic  impedance  of  the  second 
medium  is  considerably  smaller  than  the  first  and  third. 

2 .  Compressive  Wave  Velocity  in  Unsaturated  Soils 

In  two  and  three  phase  soil  media  there  are  several  factors  which  affect  the 
velocity  of  wave  propagation.  A  two  phase  soil  media  consists  of  soil  particles  and  water, 
whereas  a  three  phase  soil  media  consists  of  soil  particles,  air,  and  water.  Allen  et  al. 
(1980)  developed  a  testing  program  to  measure  the  influence  of  void  ratio  and  degree  of 
saturation  on  compression  wave  propagation.  Soils  were  tested  at  different  void  ratios,  at 
maximum  or  near  maximum  relative  density,  and  in  the  saturated  and  unsaturated 
condition.  Experimental  results  showed  that  void  ratio  was  the  most  influential  soil 
parameter  on  wave  propagation  velocity  in  a  saturated  medium.  For  a  given  sand  and  range 
of  void  ratio,  Allen  et  al.  (1980)  found  that  the  wave  velocity  generally  decreased  as  the 
void  ratio  increased.  Frequency,  hydrostatic  pressure,  temperature,  specific  gravity  of  the 
solid  particles,  and  permeability  had  less  effect  on  the  wave  propagation  in  saturated  sands. 
However,  in  the  nearly  saturated  state,  the  samration  level  was  the  most  dominant  factor. 

Ishihara  (1967)  showed  that  compression  waves  transmitted  through  the  fluid 
depended  primarily  on  the  bulk  compressibility  of  the  soil-fluid  system.  Small  amounts  of 
air  can  have  an  appreciable  influence  on  the  bulk  compressibility  of  water.  Based  on 
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mixture  theory  it  can  be  shown  that  0.10  percent  of  air  in  the  mix  will  reduce  the 
compression  wave  velocity  in  water  to  about  one  fourth  of  the  value  when  saturated 
(Richart  et.  al.,  1970). 

Hughes  and  Kelly  (1952)  investigated  the  effects  of  water  content  on 
compressional  wave  velocity  in  sandstone.  They  found  a  distinct  trend  within  the  results  in 
that  small  increases  in  saturation  beyond  about  95  percent  caused  a  sharp  drop  in  velocity 
as  the  saturation  approached  100  percent  They  also  demonstrated  that  as  confirung 
pressure  increases,  the  compressional  wave  velocity  in  the  dry  sandstone  samples  increase 
at  a  much  faster  rate  than  velocities  in  the  saturated  samples.  As  the  saturation  increases  at 
a  given  pressme,  the  velocity  at  that  pressure  decreases.  This  can  be  seen  more  clearly  in 
Figure  16.  In  the  saturation  range  of  15  percent  to  about  90  percent,  saturation  has  little 
effect  on  the  velocity  at  different  pressures.  However,  at  saturations  greater  than  90 
percent  there  is  a  considerable  decrease  in  the  velocity.  Also  at  higher  pressiues  the 
decrease  in  the  velocity  is  much  more  dramatic. 

High  strain  rate  stress  wave  propagation  tests  using  the  split-Hopkinson  pressure  bar 
(SHPB)  were  performed  by  Pierce  (1989).  Relationships  between  compressional  wave 
velocity  and  stress  transmission  with  variations  in  saturation  and  confining  stress  were 
studied.  Pierce  compacted  sands  to  a  constant  dry  density  using  four  2.54  cm  (1  inch)  lifts 
by  striking  a  steel  rod  with  a  rubber  mallet  Samples  were  fully  saturated,  and  then 
desaturated  to  a  particular  level  by  use  of  a  pressure  plate  apparatus.  Pierce  found  that 
saturations  between  zero  and  80  percent  do  not  significantly  influence  the  values  of 
compressive  wave  velocity,  transmission  ratio,  or  quasi-static  stiffness  for  Ottawa  20-30 
and  Eglin  sands  compacted  dry.  However,  retention  curves  for  Ottawa  20-30  and  Eglin 
sands  show  that  the  influence  of  capillary  forces  on  effective  stress  is  negligible.  Soil 
mictostructure  developed  during  compaction  may  be  influenced  by  the  capillary  stresses 
which  affect  soil  behavior  for  static  and  dynamic  cases. 

Ross  (1986,  1989),  Veyera  (1989),  and  Veyera  and  Fitzpatrick  (1990a) 
performed  SHPB  tests  on  dry  and  unsaturated  compacted  sands  including  Ottawa  20-30 
sand,  Ottawa  F-58  sand,  Tyndall  sand,  and  Eglin  sand.  Figures  17  and  18  show 
compressional  wave  speed  and  transmission  ratio  as  a  function  of  saturation  in  which  the 
data  have  been  normalized  to  the  average  value  at  zero  saturation  (dry  soil).  These  figures 
show  a  typical  behavior  found  in  all  soils  tested,  which  is  that  the  wave  speed  and 
transmission  ratio  increase  up  to  about  20  percent  saturation.  At  saturations  of  about  20 
percent  to  60  percent  the  wave  speed  and  transmission  ratio  remain  relatively  constant.  As 
the  saturation  increases  past  60  percent,  the  wave  speed  and  transmission  ratio  decrease. 
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Figure  16.  Dilatational  Wave  Propagation  Velocity  as  a  Function 
of  Saturation  and  Confining  Pressure  (Hughes  and  Kelly,  1952). 
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OTTAWA  20-30  SAND 


Figure  17.  Normalized  Transmission  Ratio  as  a  Function  of 
Saturation  for  Ottawa  20>30  Sand  (Veyera,  1989). 


NORMALIZED  WAVE  SPEED 


OTTAWA  20-30  SAND 


Figure  18.  Normalized  Wavespeed  as  a  Function  of  Saturation 
for  Ottawa  20-30  Sand  (Veyera,  1989). 
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The  stiffness  of  soil  is  dependent  on  the  moisture  present  during  compaction. 
Pierce  (1989)  and  Ross  (1986)  obtained  different  results  of  wave  propagation  velocity  and 
transmission  ratio  for  the  same  soils  using  different  procedures  for  moisture  addition.  The 
reason  for  the  difference  is  probably  due  to  the  sample  preparation  methods.  Pierce 
compacted  sands  dry  and  then  saturated  and  desaturated  the  sands  using  the  pressure  plate 
apparatus,  thus  the  effect  of  moisture  during  compaction  was  not  present.  Ross  compacted 
sands  for  both  dry  and  unsaturated  conditions,  thus  the  effect  of  moisture  during 
compaction  was  present  Pierce  concluded  that  soil  microstructure  may  be  influenced  by 
the  saturation  level  during  compaction  and  that  this  difference  may  change  the  state  of  stress 
in  the  sand. 

In  reviewing  previous  research  it  is  quite  evident  that  the  compressional  wave 
speed  and  stress  transmission  ratio  are  both  greatly  affected  by  the  amount  of  moisture 
present  during  compaction.  The  sample  preparation  method,  void  ratio,  density  and 
stiffness  of  the  soil  also  produce  variations  in  compressional  wave  speed  and  stress 
transmission  ratio. 
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SECTION  m 


EXPERIMENTAL  PROCEDURES 

i 

SECTION  ni  discusses  two  important  topics;  equipment  and  instrumentation 
(including  standard  and  specially  designed  items),  and  procedures  used  for  this 
investigation.  The  equipment  required  for  high  strain  rate  testing  and  sample  preparation  of 
soils  is  located  at  the  Air  Force  Engineering  and  Service  Center,  Tyndall  AFB,  Florida. 

•  ■  The  split-Hopkinson  pressure  bar  (SHPB)  and  thin  sectioning  devices  are  specially 

designed  equipment  used  to  test  and  prepare  specimens  for  microstructural  analysis. 

A.  LABORATORY  EQUIPMENT  AND  INSTRUMENTATION 

1.  Split-Hopkinson  Pressure  Bar 

The  split-Hopkinson  pressure  bar  (SHPB)  was  used  in  this  investigation  to 
measure  stress  transmission  and  compressional  wave  propagation  velocity  in  unsaturated 
Ottawa  20-30  sand.  The  strain  rates  produced  by  the  SHPB  simulate  blast  loadings  such  as 
those  generated  from  conventional  weapons  and  can  be  used  to  test  materials  in 
compression  at  strain  rates  on  the  order  of  10/sec  to  10''^/sec.  Strain  rates  for  this 
investigation  were  on  the  order  of  13/sec  to  19/sec.  Several  researchers  have  used  the 
SHPB  to  study  dynamic  stress  transmission  in  soils  (Veyera  and  Fitzpatrick,  1990b; 
Veyera,  1989;  Ross,  1989;  Felice  et  al.,  1987;  Ross  et.  al.,  1986).  A  schematic  of  a  soil 
specimen  in  the  SHPB  device  ready  for  testing  is  shown  in  Figure  19. 

The  SHPB  device  produces  a  compressive  stress  wave  in  a  steel  bar  to  load  the 
specimen  by  the  impact  of  a  striker  bar  on  an  incident  bar.  The  specimen  is  set  between 
two  steel  bars,  the  incident  and  transmitter  bar.  The  compressive  stress  wave  travels 
through  the  incident  bar  and  into  the  specimen.  A  portion  of  the  compressive  stress  wave 
is  transmitted  through  the  specimen  to  the  transmitter  bar  and  a  portion  is  reflected  back  into 
the  incident  bar.  As  the  compressive  stress  wave  travels  through  the  incident  bar, 
specimen,  and  transmitter  bar  the  SHPB  instrumentation  (strain  gages,  signal  amplifiers, 
and  digital  storage  oscilloscope)  records  the  stress  pulse  thus  allowing  the  determination  of 
the  stress  transmission  ratio  and  compressional  wave  speed.  A  detailed  discussion  of  the 
SHPB  theory  and  principles  is  given  by  Ross  (1989). 
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2.  MTS  Quasi-Static  Loading  Device 


A  Materials  Testing  System  (MTS)  device  was  used  in  this  investigation  to  quasi- 
statically  load  the  dry  and  unsaturated  soil  specimens.  The  constrained  modulus,  which  is 
a  measure  of  soil  stiffness,  was  obtained  by  using  this  system  which  can  be  used  to  test 
materials  in  the  quasi-static  strain  rate  range,  i.e.,  strain  rates  of  10^/sec  to  lO'^sec.  For 
this  investigation  the  maximum  load  was  1237  Newtons  (55(X)  pounds)  and  the  minimum 
load  value  was  set  at  zero  Newtons.  Several  cycles  of  loading  and  unloading  can  be 
achieved  if  so  desired.  For  this  investigation  loading  rates  for  loading  cycles  were  on  the 
order  of  4.45  kN/min  (1000  Ib/min)  and  unloading  cycles  were  on  the  order  of  13.34 
kN/min  (3000  Ib/min). 

3.  Specimen  Container 

Three  different  specimen  containers  were  designed  specifically  for  this  research. 
Each  container  was  15.24  cm  (6  inches)  long  with  an  inside  diameter  of  5.08  cm  (2.00 
inches).  The  specimen  containers  were  designed  such  that  they  would  be  compatible  with 
the  SHPB.  The  difference  between  the  three  specimen  containers  was  in  the  container  wall 
thicknesses  which  were  2.54  cm  (1.00  inches),  1.27  cm  (0.50  inches),  and  0.1588  cm 
(0.0625  inches).  For  this  investigation,  the  2.54  cm  (1.00  inches)  wall  thickness  specimen 
container  was  used.  The  other  stainless  steel  specimen  containers  were  used  in  a 
continuation  of  this  work  to  simulate  different  degrees  of  boundary  stiffness  that  might 
typically  be  encountered  in  the  field,  the  results  of  which  are  not  part  of  this  research. 

Pedestals,  wafer  seals,  and  compaction  pistons  were  designed  to  fit  with  the 
specimen  containers  during  the  compaction  procedure  using  a  standard  Proctor  hammer 
ASTM  D  698-90  (ASTM  1990).  A  pedestal  was  used  to  hold  the  wafer  seal  in  place  such 
that  the  top  of  the  wafer  seal  was  5.08  cm  (2.00  inches)  up  from  the  bottom  of  the 
specimen  container.  The  wafer  seals  with  Viton  O-rings  were  used  to  minimize  moisture 
leakage  fiom  the  base  of  the  specimen  containers  during  compaction  and  at  both  ends  of  the 
specimen  for  containment  during  SHPB  testing.  Four  compaction  pistons  were  designed 
such  that  four  individual  layers  each  of  2.54  cm  (1.00  inches)  could  be  compacted.  Figure 
20  shows  an  assembly  drawing  of  the  specimen  container. 

4.  Sectioning  Apparatus 

Two  sectioning  devices  were  used  in  developing  specimen  sections  for 
microstructural  analysis  which  are:  a  water-cooled  Buehler  Metaseu  Powermet  Automatic 
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Abrasive  Cutter  (AAC)  used  for  initial  rough  cutting,  and  a  Buehler  Datamet 
Micro-Processor  Grinding  and  Polishing  unit  used  for  grinding  and  polishing.  A  section  is 
a  cut  made  through  the  specimen  which  exposes  the  preserved  soil  structure.  Cutting  plane 
locations  are  discussed  in  detail  in  later  SECTION  III.  The  AAC  is  stocked  with 
interchangeable  cutting  wheels  and  is  used  to  make  an  initial  rough  cut  through  the 
specimen.  An  aluminum  alignment  device  was  fabricated  such  that  uniform  cuts  through 
the  specimen  could  be  obtained.  The  alignment  device  had  dimensions  of  7.62  x  7.62  x 
0.32  X  12.70  cm  long  (3  x  3  x  0.125  x  5  inches).  Upon  proper  alignment,  the  AAC 
automatically  makes  a  complete  cut  through  the  specimen. 

The  grinding  and  polishing  unit  is  used  as  a  final  step  to  prepare  specimens  prior 
to  the  microstructural  analysis  and  is  used  in  conjunction  with  the  following  items:  (a)  a 
grooved  cast  iron  polishing  wheel  with  #600  Polishing  Grit,  (b)  a  smooth  cast  iron 
polishing  wheel,  (c)  a  polishing  cloth,  (d)  Metadi  Fluid,  (e)  diamond  paste,  and  (f) 
Micropolish  "C".  The  grinding  and  polishing  procedure  involves  three  separate  steps:  (1) 
initial  grinding,  (2)  rough  polishing,  and  (3)  fine  polishing.  These  three  steps  are  briefly 
outlined  in  SECTION  HI  of  this  thesis.  A  more  detailed  description  of  the  AAC,  grinding 
and  polishing  unit,  and  procedures  is  given  by  Veyera  and  Fitzpatrick  (1990a). 

5 .  Hardware  and  Software  for  Microstructural  Analysis 

The  mathematical  analysis  of  microstructural  data  is  performed  using  Macintosh 
nsi  and  Macintosh  SE/30  personal  computers.  A  Summagraphics  (1991)  digitizing  tablet 
was  interfaced  with  the  Macintosh  to  perform  particle  orientation  analyses.  Two  con^uter 
programs  developed  by  Rock  Ware,  Inc.  (1991),  "Digitize"  and  "Rosy,"  were  used  to 
determine  the  orientation  of  hundreds  of  particles  based  on  particle  long-axis  orientation. 
The  "Digitize"  program  allows  the  user  to  digitize  line  segments  which  are  then  transformed 
to  azimuth  data  and  stored  for  statistical  analysis.  The  program  "Rosy"  generates  angular 
frequency  histograms  or  rose  plots,  and  performs  a  complete  statistical  analysis  of  the  data 
including  the  number  of  particles  measured,  vector  mean,  confidence  angle,  and  vector 
magnitude.  Rosy  also  uses  the  Rayleigh  test  of  significance  to  determine  if  a  preferred 
particle  orientation  exists  within  the  data.  The  output  from  the  program  is  used  to 
determine  if  preferred  particle  orientation  exists  in  the  samples  as  a  function  of  saturation, 
compaction  energy,  confinement  conditions,  and  load  intensity. 
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B .  LABORATORY  PROCEDURES 
1 .  Specimen  Preparation 
a.  Compaction 

Specimens  of  Ottawa  20-30  sand  were  dynamically  compacted  to  a  constant 
dry  density  of  1.715  g/cm^  (107.0  pcf)  at  varying  degrees  of  saturation  in  the  2.54  cm 
thick  walled  stainless  steel  specimen  container  using  the  standard  Proctor  hammer.  An 
assembly  drawing  of  a  typical  specimen  container  is  shown  in  Figure  20. 

A  standard  Proctor  hammer,  ASTM  D  698-90  (ASTM  1990),  was  used  to 
apply  a  controlled  amount  of  compactive  effort  per  impact  to  each  soil  specimen  (7.5  Joules 
or  5.5  ft-lbs  per  impact).  Prior  to  compaction,  each  specimen  container  and  bottom  steel 
wafer  insert  were  sprayed  with  a  dry  Teflon®  lubricating  agent  (TFL).  All  test  specimens 
were  formed  using  four  individually  compacted  layers  of  equal  weight  such  that  a  final 
specimen  length  of  10.16  cm  (4.00  inches)  would  be  obtained.  Four  individual  lifts,  each 
having  a  final  thickness  of  2.54  cm  (1.00  inches),  were  used.  Saturations  were  varied 
from  zero  (dry)  to  near  80  percent.  At  saturations  of  zero,  each  layer  of  dry  soil  was 
poured  directly  into  the  tube  and  compacted. 

In  preparing  moist  specimens,  the  required  amount  of  water  for  a  given 
degree  of  saturation  (at  final  density  when  compacted)  was  added  to  the  originally  dry  soil, 
thoroughly  mixed  in  and  then  allowed  to  equilibrate  before  compacting.  Since  the 
specimens  ranged  in  degree  of  sattu^tion  from  zero  to  about  80  percent,  the  tests  were 
conducted  on  unsaturated  specimens.  With  respect  to  this  study,  the  term  "unsaturated" 
implies  that  both  continuous  air  and  water  phases  exist  in  the  soil  (e.g.,  there  are  no 
isolated  air  or  water  pockets  in  the  grain  structure).  For  most  soils,  this  generally  occurs  at 
saturations  less  than  about  85  percent  (Corey,  1977). 

After  impact  compaction  using  the  standard  Proaor  hammer  was  completed, 
some  specimens  were  tested  in  the  SHPB  device.  The  same  saturation  range  was  used  for 
both  NO  SHPB  (compacted  only)  and  SHPB  (compacted  and  tested  in  SHPB  device) 
specimens  in  moist/moist  tests  (i.e.,  zero  to  about  80  percent).  Upon  completion  of  SHPB 
tests,  specimens  were  dried  in  an  oven  for  a  minimum  of  12  hours  at  1 10®  C  to  remove  all 
moisture.  This  procedure  applied  to  specimens  dynamically  compacted  with  the  standard 
Proctor  hammer  only  (NO  SHPB),  and  to  those  specimens  dynamically  compacted  with  the 
standard  Proctor  hammer  and  tested  using  high  strain  rates  (SHPB). 
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Other  specimens  were  dried  in  an  oven  after  impact  compaction  for  a 
minimum  of  12  hours  at  1 10“  C  to  remove  all  moisture  and  then  tested  in  the  SHPB  device 
(moist/dry  tests).  The  was  done  to  determine  the  influence  of  moisture  during  testing. 

b .  Split-Hopkinson  Pressure  Bar  Tests 

For  this  investigation,  the  SHPB  device  was  used  to  study  dynamic  stress 
transmission  characteristics  and  compressional  wave  velocity  as  a  function  of  saturation 
and  microstracture  in  Ottawa  20-30  sand.  The  SHPB  has  been  used  by  several  researchers 
to  study  dynamic  stress  transmission  in  soils  (e.g.,  Veyera  and  Fitzpatrick,  1990b;  Veyera, 
1989;  Pierce,  1989;  Felice  et.  al.,  1987;  Ross  et.  al.,  1986).  Dry  and  unsaturated 
specimens  were  compacted  using  the  standard  Proctor  hammer  and  tested  at  high  strain 
rates  using  the  SHPB.  Transmission  ratio  and  wave  speed  at  several  saturations  were 
measured. 

The  conditions  of  the  testing  program  using  the  SHPB  were  as  follows:  (a)  a 
20.32  cm  (8.00  inches)  long,  5.08  cm  (2.00  inches)  diameter  steel  striker  bar  was  used; 

(b)  the  average  velocity  of  the  striker  bar  was  approximately  1305  cm/sec  (514  inches/sec); 

(c)  the  diameters  of  the  incident  bar,  specimen,  and  transmitter  bar  were  all  5.08  cm  (2.00 
inches);  and  (d)  incident  stress  levels  were  approximately  290,000  kPa  (42,059  psi)  in  the 
incident  bar. 


c.  (Juasi-Static  Uniaxial  Q>mpression  Tests 

Specimens  were  compacted  using  four  layers  in  the  2.54  cm  (1.00  inches) 
walled  tube  using  the  standard  Proctor  hammer  as  shown  in  Figure  20.  The  top  flange 
plate  was  removed  and  the  compaction  piston  used  for  compacting  Layer  1  was  placed  on 
the  soil.  The  specimen  was  placed  in  the  MTS  device  and  two  complete  load-unload  cycles 
were  applied.  Maximum  and  minimum  loads  were  recorded  on  a  data  sheets  and  analog 
stress-strain  plots  were  generated  by  the  MTS  device. 

For  comparison  purposes,  constrained  modulus  values  in  uniaxial 
compression  were  determined  from  the  analog  stress-strain  plots  by  taking  a  secant  line  at  a 
predetermined  stress  level.  The  constrained  modulus  was  determined  for  each  saturation 
by  constructing  a  secant  line  between  stress  levels  of  zero  and  10,974  kPa  (1,592  psi). 
The  reason  for  choosing  a  quasi-static  stress  level  of  10,974  kPa  was  that  SHPB  stress 
levels  in  the  soil  specimens  were  on  the  order  of  1 1,000  kPa.  This  analysis  was  performed 
to  obtain  information  about  variation  in  soil  stiffness  with  saturation.  The  data  are  useful 
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because  the  stiffness  of  the  soil  directly  affects  the  wave  speed  and  transmission  ratio. 
Equation  18  defines  the  constrained  modulus  as: 


Lo 


(18) 


where  A<J  is  the  change  in  compressive  stress,  Ae  is  the  change  in  axial  strain,  AP  is  the 
change  in  applied  static  load,  A  is  the  cross  sectional  area,  Lq  is  the  initial  compacted 
specimen  length,  and  AL  is  the  change  in  specimen  length. 

d.  Epoxying  of  Conq)acted  Specimens 

A  number  of  commercially  available  epoxies  and  bonding  agents  are  routinely 
used  in  standard  petrographic  analyses  of  geologic  materials.  In  this  study,  several  criteria 
were  set  forth  in  selecting  a  suitable  bonding  agent,  including  the  following  requirements: 

(1)  strong  bonding  with  soil  grains  and  the  ability  to  remain  intact  during 
cutting  and  polishing  procedures; 

(2)  relatively  low  viscosity  material; 

(3)  straightforward  and  simple  preparation  procedure; 

(4)  no  shrinkage  during  curing; 

(5)  relatively  rapid  curing  time  under  controlled  conditions;  and 

(6)  bonding  agent  can  be  easily  colored  or  dyed  (as  necessary). 

A  commercially  available  bonding  agent,  "EPOTEK  301®,"  an  epoxy 
produced  by  Epoxy  Technologies,  Inc.,  was  used  to  preserve  the  structure  of  compacted 
specimens  for  the  microstmctural  analysis.  EPOTEK  301  is  a  high  strength,  spectrally 
transparent,  low  viscosity  (100  cps),  two  component  epoxy.  Curing  time  takes  about  1 
hour  in  a  dry,  temperature  controlled  oven  at  65  ±  2°C.  Specimens  can  also  be  cured  at 
room  temperature  overnight  When  fully  cured,  the  epoxy  is  strongly  bonded  to  the 
specimen  grains  and  noaintains  the  soil  structure  during  cutting  and  polishing  operations. 

Since  the  epoxy  is  clear,  it  was  necessary  to  introduce  color  to  provide 
adequate  contrast  between  the  pores  and  grains.  Several  coloring  agents  and  dyes  were 
investigated  and  the  most  acceptable  results  were  obtained  using  commercially  available 
biological  stains  in  powdered  form  (Sudan  IV-dark  red;  Sudan  Black-dark  blue  to  black). 

The  epoxy  used  does  not  mix  or  cure  if  water  is  present,  therefore  specimens 
initially  compacted  moist  at  varying  degrees  of  saturation  in  the  specimen  containers  were 


45 


dried  in  a  temperature  controlled  oven  at  110°C  for  at  least  12  hours.  Specimens  were 
then  allowed  to  cool  to  room  temperature  before  epoxying.  Careful  handling  of  the 
specimen  container  was  essential  during  all  phases  of  preparation  to  insure  that  the 
compacted  soil  structure  was  not  disturbed  (i.e.,  by  jarring,  vibration,  tipping,  etc.). 

Material  quantities  of  EPOTEK  301  epoxy  were  determined  based  on  a 
calculation  of  the  available  pore  space  for  a  medium  granular  soil  (Ottawa  20-30  sand) 
compacted  to  a  void  ratio  of  0.545.  Allowances  were  made  for  additional  material  required 
to  cover  particles,  minor  losses  in  handling,  etc.,  with  a  minimum  of  waste.  The  mixing 
ratio  for  the  epoxy  component  parts  is  4:1  (resin  to  hardener)  as  specified  by  the 
manufacturer.  Figure  21  shows  a  typical  specimen  being  epoxied.  Details  of  the 
laboratory  procedure  for  epoxying  compacted  specimens  are  given  by  Veyera  and 
Htzpatrick  (1990a).  The  general  procedure  is  briefly  outlined  as  follows: 

(1)  Specimens  compacted  moist  are  dried  overnight  in  a  temperature 
controlled  oven.  Specimens  compacted  dry  are  epoxied  immediately  (Step  2). 

(2)  EPOTEK  301  epoxy  is  carefully  introduced  into  the  specimen  under 
atmospheric  pressure  to  saturate  the  void  spaces.  The  specimen  is  then  cured  in  a 
temperature-controlled  oven  at  65  ±  2‘’C  for  about  1  hour. 

(3)  After  curing,  the  specimen  is  extruded  finom  the  specimen  container. 

2.  Specimen  Sectioning 

a.  Extrusion 

After  the  epoxy  has  fully  cured,  the  next  step  is  to  extrude  each  specimen 
from  the  specimen  container.  For  oven-cured  specimens,  the  containers  must  be  cooled 
before  extrusion.  The  procedure  uses  a  modified  Marshall  Testing  machine  for  extruding 
specimens,  and  requires  a  loading  piston,  loading  support  blocks  and  a  flange  plate. 

b.  Initial  Rough  Cutting 

A  water-cooled  Buehler  Metaseu  Powermet  ACC  was  used  for  the  initial 
rough  cutting  of  epoxied  specimens.  A  simple  aluminum  specimen  alignment  device  for 
making  uniform  cuts  was  fabricated.  Once  a  specimen  is  properly  aligned,  the  ACC 
automatically  makes  a  complete  cut  through  the  specimen.  Each  cut  section  is  clearly 
marked  to  record  the  original  orientation.  Figure  22  shows  the  location  of  the  various 
cutting  planes. 
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Cutting  the  specimens  for  use  in  microstructural  analysis  involves  several 
steps  which  are  listed  below; 

(1)  first  cut  removes  about  0.3125  cm  (0.125  inches)  from  each  end  of  the 
specimen; 

(2)  second  cut  splits  the  specimen  in  half  along  its  long-axis  (Figure  23); 

(3)  third  cut  splits  one  of  the  specimen  halves  along  its  diameter  (Figtire  24). 

(4)  a  vertical  cut  about  0.3125  cm  (0.125  inches)  in  thickness  is  made  along 
the  long-axis  of  one  of  the  halves  from  the  initial  cuts.  Step  2  (Figure 
25); 

(5)  a  horizontal  cut  about  0.3125  cm  (0.125  inches)  in  thickness  is  made 
along  the  diameter  of  one  of  the  halves  from  the  initial  cuts.  Step  3 
(Figure  26). 

c .  Grinding  and  Polishing 

The  final  preparation  procedure  for  the  specimens  prior  to  microstructural 
analysis  involves  several  steps  and  uses  a  Buehler  Datamet  Micro- Efrocessor  Grinding  and 
Polishing  Unit.  Details  of  the  laboratory  procedure  for  the  grinding  and  polishing  of 
contacted  specimens  are  given  by  Veyera  and  Fitzpatrick  (1990a).  The  general  procedure 
is  briefly  outlined  as  follows: 

(1)  Initial  Grinding  -  This  step  is  used  to  true  the  planar  surface  obtained 
during  the  rough  cutting  process  and  to  prepare  the  specimen  surface  for 
polishing. 

(2)  Rough  Polishing  -  This  step  is  used  to  rough  polish  the  specimen  section 
surface  prior  to  detailed  finish  polishing. 

(3)  Fine  Polishing  -  This  step  is  used  to  develop  a  smooth,  highly  polished 
surface  finish  that  provides  adequate  contrast  between  the  specimen 
grains  and  the  expoy-filled  pores. 

Once  the  final  polishing  has  been  completed,  specimens  are  examined  under  a 
microscope  to  check  if  the  surface  has  been  properly  prepared  and  if  adequate  contrast 
between  pore  spaces  and  individual  grains  has  been  developed.  Properly  prepared 
specimens  can  then  be  used  in  microstructural  analysis  studies. 

3.  Photomicrographs  of  Sectioned  Specimens 

Photographs  of  sectioned  specimens  are  used  to  determine  axial  ratio  and  particle 
long-axis  orientation  based  on  the  procedure  developed  by  Oda  (1972a)  and  Campbell 
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Figure  23.  Initial  Specimen  Cut  Along  the  Long  Axis  ( Section  A- A). 


Figure  24.  Initial  Specimen  Cut  Along  theDiametral  Axis  (Section  B  -  B). 
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Figure  25.  Final  Cut  of  Specimen  Section  Along  the  Long  Axis. 


Figure  26.  Final  Cut  of  Specimen  Cut  Along  the  Diametral  Axis . 


(1985).  A  35  mm  camera  with  special  microscope  attachments  was  mounted  on  a  reflected 
light  microscope.  Direct  lighting  was  used  to  illuminate  the  cut  specimen  section  from 
above.  The  microscope  was  set  up  such  that  sections  of  each  specimen  could  be 
photographed  using  a  magnification  of  lOx  on  all  specimen  sections.  Standard  35  mm 
color  slide  film  (ASA  100)  was  used  to  photograph  specimens.  The  optimum  shutter  speed 
was  1/60  second,  which  was  obtained  through  a  careful  study  of  photographs  using 
various  shutter  speeds  and  provided  the  optimum  visual  contrast  between  sand  grains  and 
epoxy. 

A  series  of  photographs  was  taken  in  the  vertical  and  horizontal  directions.  The 
grid  layout  and  the  identification  schemes  for  Sections  A-A  and  B-B  are  shown  in  Figure 
27 .  A  record  of  frame  number  and  corresponding  grid  number  was  maintained  for  proper 
identification  of  each  slide.  A  total  of  19  grids  were  photographed  for  each  specimen. 
Grids  a-g  represent  the  horizontal  direction  on  Section  A-A,  grids  h-o  represent  the  vertical 
direction  on  Section  A-A,  and  grids  p-s  represent  the  semicircular  section.  The  properly 
identified  slides  are  then  used  to  analyze  the  microstructure  of  the  compacted  specimens. 

4.  Microstructural  Analysis 

To  perform  particle  long-axis  orientation  and  axial  ratio  analyses,  a  two- 
dimensional  image  of  the  grain  and  void  arrangements  is  needed.  The  two-dimensional 
image  is  obtained  by  projecting  the  slides  through  a  standard  slide  projector  onto  a  21.59 
cm  by  27.94  cm  (8.5  in  x  1 1  in)  sheet  of  paper  attached  to  the  digitizing  tablet  (see  Figure 
28).  Individual  grains  are  digitized  by  locating  the  longest  axis  of  the  particle.  Digitized 
data  files  are  clearly  identified  according  to  degree  of  samration,  grid  number,  and  whether 
or  not  the  specimen  was  tested  in  the  SHPB  device. 

a.  Axial  Ratio 

To  determine  the  axial  ratio  of  an  individual  particle  and  a  group  of  particles, 
the  procedure  suggested  by  Oda  (1972a)  was  foUowed: 

(1)  Determine  the  length  of  the  particle's  shortest  axis,  L2,  by  finding  the 
largest  possible  circular  that  can  be  inscribed  on  the  particle  (see  Figure 
29a). 

(2)  Determine  the  length  of  the  particles  longest  axis,  Li,  by  finding  the 
smallest  possible  circular  that  can  be  circumscribed  on  the  particle  (see 
Figure  29b). 
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a.  Short  Axis  (1^  )  Determination. 


Individual  Particle 


b.  Long  Axis  (L^ )  Determination. 


Figure  29.  Determination  of  Axial  Ratio  Parameters  for  an  Individual  Particle. 
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The  axial  ratio  for  a  group  of  particles,  Uavg ,  is  given  by  Equation  (5).  The 
axial  ratio  for  an  individual  particle,  Uj ,  is  given  by  Equation  (19). 


Hs 


(19) 


Axial  ratio  data  are  useful  in  classifying  particle  shape  and  also  provide  an 
indication  of  the  potential  particle  arrangements  that  can  be  expected  in  a  compacted 
specimen.  An  example  would  be  that  as  particles  approach  a  spherical  shape,  particle 
arrangement  become  more  random  Oda  (1972a)  indicated  that  an  axial  ratio  of  0.7  or  less 
was  required  for  preferred  particle  long-axis  orientation.  Table  1  provides  classes  and 
descriptions  of  particle  shapes  based  on  axial  ratio  data  (roundness  ratio).  The  Ottawa  20- 
30  sand  used  in  this  research  would  be  classified  as  well  rounded  according  to  Table  1. 

b.  Particle  Long-Axis  Orientation 

Particle  long-axis  orientation  is  a  standard  method  used  in  the  microstructural 
analysis  of  granular  soils.  The  procedure  used  to  determine  long-axis  orientation  of 
individual  particles  was  based  on  the  method  developed  by  Campbell  (1985)  which 
involves  determining  the  longest  axis,  Li,  and  the  shortest  axis,  L2,  of  a  particle.  Figure 
30  shows  the  determination  of  particle  long-axis,  Li,  and  orientation  angle,  0,  with  respect 
to  a  set  of  Cartesian  coordinate  system  reference  axes. 

The  two  endpoints  of  the  longest  axis  are  digitized  and  read  by  the  "Digitize" 
program  as  line  segment  data.  The  orientation  angle,  9,  of  an  individual  particle  is 
determine  by  converting  line  segment  data  to  azimuth  data.  The  azimuth  data  is  read  into 
the  "Rosy"  program  which  then  calculates  the  Vector  Mean,  Vector  Magnitude,  Rayleigh 
test  of  significance,  and  also  generates  angular  frequency  histograms  (rose  plots)  for  a 
group  of  particles.  For  each  rose  plot  developed,  orientation  measurements  (class 
intervals)  were  made  at  lO-degree  intervals  using  a  95  percent  confidence  interval.  The 
output  firom  Rosy  allows  the  user  to  determine  if  a  preferred  particle  orientation  exists  in  the 
specimen. 

Particle  long-axis  orientation  analyses  were  performed  in  the  vertical  and 
horizontal  directions  on  Section  A- A  and  on  the  semicircular  Section  B-B  (see  Figure  27). 
The  vertical  analysis  was  further  extended  to  examine  long-axis  orientations  in  each 
compacted  layer.  Vertical  grids  (h-o),  were  separated  into  four  layers:  (a)  Layer  #1  -  grids 
o,  n,  and  m;  (b)  Layer  #2  -  grids  m,  1,  and  d:  (c)  Layer  #3  -  grids  d,  k,  and  j;  (d)Layer  #4  - 
grids  j,  i,  and  h.  Results  from  the  long-axis  orientation  analyses  provide  information  about 
preferred  particle  orientations  that  may  occur  either  locally  or  globally  within  the  specimen. 
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SECTION  IV 


EXPERIMENTAL  RESULTS 

This  chapter  contains  five  sections  which  show  the  experimental  results  obtained  in 
the  present  investigation.  The  first  section  gives  information  pertaining  to  the  physical  soil 
properties.  The  second  section  deals  with  specimen  compaction.  The  third  section 
presents  the  Split-Hopkinson  pressure  bar  results.  Section  four  gives  the  results  of  quasi¬ 
static  uniaxial  compression  tests,  and  the  final  section  is  devoted  to  results  obtained  in  the 
microstructural  analysis. 

Most  of  the  figures  presented  in  this  chapter  show  saturations  ranging  from  zero  to 
about  75  percent  Several  other  figures  include  data  for  saturations  of  about  95  percent. 
Specimens  prepared  for  the  microstructural  analysis  had  saturations  from  zero  to  about  75 
percent  Blow  count  compressional  wave  speed,  and  stress  transmission  ratio  data  were 
obtained  for  saturations  from  zero  to  about  95  percent 

A.  PHYSICAL  PROPERTIES  OF  SOIL 

A  commercially  available  granular  soil  designated  as  "Ottawa  20-30  sand"  was  used 
in  this  investigation.  Approximately  227  kg  (500  pounds)  of  the  sand  were  obtained  from 
the  Ottawa  Silica  Company  and  random  samples  were  taken  from  the  bulk  quantity.  The 
material  is  a  uniformly  graded,  subrounded  to  rounded,  medium  sand  with  no  fines  and  is 
classified  as  SP  (poorly  graded  sand)  according  to  the  Unified  Soil  Qassification  System 
(USCS).  Physical  index  properties  for  the  Ottawa  20-30  sand  are  summarized  in  Table  2 
and  a  grain  size  distribution  cmve  is  shown  in  Figure  31. 

B .  SPECIMEN  COMPACTION 

Each  specimen  of  Ottawa  20-30  sand  was  dynamically  compacted  to  a  constant  dry 
density  of  1.715  g/cm^  (107.0  pcf)  using  a  standard  Proctor  hammer.  The  amount  of 
compactive  effort  required  to  obtain  to  this  constant  dry  density  varied  with  the  amount  of 
moisture  present.  Records  of  the  number  of  blows  required  to  achieve  this  dry  density  for 
each  layer  and  saturation  were  maintained.  Table  3  summarizes  the  blow  count  data  and 
Figure  32  shows  the  average  number  of  blows  as  a  function  of  saturation  required  to 
compact  the  Ottawa  20-30  sand  to  a  constant  dry  density  of  1.715  g/cm^.  Results  are 
based  on  data  collected  fix)m  56  specimens  compacted  in  the  2.54  cm  (1.00  in)  thick  walled 
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TABLE  2.  Physical  Properties  of  Ottawa  20-30  Sand. 

uses  Classification  SP 

Specific  Gravity  2.65 

D50  particle  size  0.70  mm 

aCu  1-40 

bCc  1.03 

^Percent  passing  #100  sieve  <1  % 

<lMaximum  dry  density  1,763  kg/m3 

<lMinimum  dry  density  1 ,587  kg/m^ 

Maximum  void  ratio  0.669 

Minimum  void  ratio  0.504 


Note:  a<;k)efficient  of  Uniformity 

^Coefficient  of  Curvature 


CU.S.  Standard  Sieve 
<lData  from  Ottawa  Silica  Sand 
Company 


specimen  container.  Figure  32  indicates  that  the  total  number  of  blows  sharply  increases 
from  zero  to  about  20  percent  saturation,  then  slightiy  increases  from  about  20  percent  to 
about  40  percent  saturation  and  decreases  thereafter  with  increasing  saturation. 

Figure  33  shows  the  total  number  of  blows  required  to  compact  the  Ottawa  20-30 
sand  to  a  constant  dry  density  of  1.715  g/cm^  for  each  individual  2.54  cm  (1.00  in)  thick 
layer.  The  data  show  that  the  greatest  compaction  effort  is  required  in  the  first  layer  at  all 
saturations.  The  first  layer  shows  a  similar  trend  as  that  found  in  Figure  32.  Layers  2  and 
3  show  a  similar  trend  with  respect  to  each  other  in  that  the  total  number  of  blows  increases 
from  zero  to  about  20  percent  saturation.  As  the  saturation  increases  beyond  20  percent, 
the  total  number  of  blows  steadily  decreases.  In  Layer  4,  the  total  number  of  blows  from 
about  20  percent  to  about  75  percent  saturation  is  relatively  unchanged. 

C.  SPLIT-HOPKINSON  PRESSURE  BAR  TESTS 

High  strain  rate  compressive  tests  were  conducted  on  Ottawa  20-30  sand  using  the 
SHPB  device.  Compressional  wave  speed  and  stress  transmission  ratio  were  studied  as  a 
function  of  saturation  and  microstructure.  Some  specimens  were  compacted  moist  and 
tested  moist  in  the  SHPB  (moist-moist  tests)  and  others  were  compacted  moist  and  tested  in 
the  SHPB  after  drying  (moist-dry  tests). 

Tables  4  and  5  give  the  raw  data  obtained  from  SHPB  tests  for  compressional  wave 
speed  and  stress  transmission  ratio  for  specimens  which  were  compacted  moist  and  tested 
moist  (moist-moist)  and  specimens  which  were  compacted  moist  and  tested  dry  (moist- 
dry),  respectively.  Figures  34  through  37  show  the  variation  within  the  data  about  the 
average  compressional  wave  speed  and  stress  transmission  ratio  as  a  function  of  saturation. 
Figures  34  and  35  show  compressional  wave  speed  and  stress  transmission  ratio  as  a 
function  of  saturation  from  data  obtained  in  27  SHPB  tests  for  specimens  compacted  moist 
and  tested  moist  (moist-moist).  Figure  34  shows  that  the  average  compressional  wave 
speed  increases  from  zero  to  about  20  percent  saturation.  The  compressional  wave  speed  is 
relatively  constant  from  20  percent  to  about  75  percent  saturation,  and  slightly  decreases 
thereafter.  Figure  35  shows  that  the  stress  transmission  ratio  increases  from  zero  to  about 
20  percent  saturation  and  has  a  sharp  increase  from  about  20  percent  to  about  40  percent 
saturation.  The  stress  transmission  ratio  then  steadily  decreases  from  about  40  percent  to 
about  95  percent  saturation. 

Figures  36  and  37  show  compressional  wave  speed  and  stress  transmission  ratio  as  a 
function  of  saturation  from  data  obtained  in  25  SHPB  tests  for  specimens  compacted  moist 
and  tested  dry  (moist-dry).  Figure  36  shows  that  the  compressional  wave  speed  increases 
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TABLE  3.  BLOW  COUNT  TEST  RESULTS  FOR  OTTAWA  20-30  SAND 

COMPACTED  USING  THE  STANDARD  PROCTOR  HAMMER. 


Container  Wall  Thickness  =  2.54  cm  (1.00  in) 

Dry  Density  =  1.715  g/cc  (107  pcf) 

Void  Ratio  =  0.545 

Sample  Lwigth  =  10.16  cm  (4.00  in) 

Compacted  four  individual  2.54  cm  (1 .00  in)  layers 

Testing  performed  by  Dr.  George  E.  Veyera  and  Blaise  J.  Fitzpatrick  of  the  University  of  Rhode  Island. 
Research  paformed  at  AFESC/RDCM,  Tyndall  AFB,  Florida. 

Date:  Summer  1990 


LAYER  &  NUMBER  OF  BLOWS 

TOTAL  NUMBER 

1 

2 

2 

OF  BLOWS 

0.0 

11 

6 

10 

4 

31 

0.0 

9 

8 

8 

3 

28 

0.0 

8 

7 

10 

5 

30 

0.0 

11 

9 

6 

4 

30 

0.0 

10 

4 

8 

4 

26 

0.0 

18 

10 

12 

5 

45 

0.0 

8 

7 

8 

4 

27 

0.0 

16 

7 

8 

5 

36 

0.0 

8 

5 

5 

6 

24 

0.0 

18 

10 

12 

6 

46 

0.0 

10 

5 

12 

4 

31 

0.0 

15 

14 

18 

11 

58 

0.0 

13 

12 

24 

5 

54 

0.0 

19 

9 

12 

9 

49 

AVERAGE 

12 

8 

11 

5 

37 

STD  DEV 

3.89 

2.71 

4.79 

2.09 

10.89 

LAYER  &  NUMBER  OF  BLOWS _ 

TOTAL  NUMBER 

S% 

1 

2 

2 

OF  BLOWS 

18.7 

16 

21 

23 

10 

70 

18.7 

25 

22 

34 

16 

97 

18.7 

29 

25 

24 

12 

90 

18.7 

34 

29 

20 

19 

102 

18.7 

35 

29 

27 

11 

102 

18.7 

40 

29 

37 

14 

120 

18.7 

56 

27 

40 

13 

136 

18.7 

55 

39 

24 

16 

134 

18,7 

55 

40 

50 

13 

158 

18.7 

40 

32 

38 

15 

125 

18.7 

53 

34 

33 

13 

133 

AVERAGE 

40 

26 

32 

14 

115 

STD  DEV 

13.00 

3.39 

8.70 

2.44 

24.14 
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TABLES,  (continued) 


s% 

1 

2 

2 

i 

OF  BLOWS 

37.4 

43 

28 

26 

12 

109 

37.4 

42 

30 

22 

15 

109 

37.4 

37 

27 

37 

15 

116 

37.4 

58 

25 

30 

10 

123 

37.4 

53 

35 

24 

10 

122 

37.4 

60 

30 

29 

15 

134 

37.4 

72 

31 

31 

15 

149 

37.4 

70 

30 

32 

20 

152 

37.4 

45 

27 

30 

15 

117 

37.4 

56 

31 

26 

17 

130 

AVERAGE 

54 

29 

29 

14 

126 

STD  DEV 

1127 

3.13 

4.12 

2.91 

14.38 

S% 

1 

2 

2 

OF  BLOWS 

56.2 

35 

26 

25 

8 

94 

56.2 

31 

22 

25 

9 

87 

56.2 

41 

20 

22 

8 

91 

56.2 

65 

26 

30 

10 

131 

56.2 

36 

23 

27 

8 

94 

56.2 

57 

27 

28 

13 

125 

56.2 

50 

22 

34 

15 

121 

56.2 

56 

25 

23 

9 

113 

56.2 

41 

24 

27 

18 

no 

AVERAGE 

46 

24 

27 

11 

107 

STD  DEV 

11.02 

2.52 

3.46 

3.41 

15.40 

S% 

1 

2 

2 

i 

OF  BLOWS 

74.9 

25 

21 

25 

8 

79 

74.9 

54 

IS 

21 

9 

99 

74.9 

47 

16 

17 

12 

92 

74.9 

70 

28 

10 

10 

118 

74.9 

70 

17 

15 

12 

114 

74.9 

62 

25 

26 

11 

124 

74.9 

52 

16 

24 

8 

100 

74.9 

60 

18 

15 

10 

103 

74.9 

75 

19 

20 

11 

125 

74.9 

55 

21 

18 

11 

105 

74.9 

50 

20 

11 

10 

91 

74.9 

45 

22 

27 

11 

105 

AVERAGE 

55 

20 

19 

10 

105 

STD  DEV 

12.98 

3.72 

5.48 

1.30 

13.28 
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TABLES,  (continued) 


LAYER  &  NUMBER  OF  BLOWS _  TOTAL  NUMBER 


1 

1 

2 

i 

OF  BLOWS 

93.6 

18 

11 

12 

10 

51 

93.6 

12 

16 

15 

11 

54 

93.6 

23 

16 

13 

12 

64 

AVERAGE 

18 

14 

13 

11 

56 

STD  DEV 

4.50 

2.36 

125 

0.82 

5.56 

AVERAGE 

LAYER  &  AVERAGE  NUMBER  OF  BLOWS _  TOTAL  NUMBER 


S% 

1 

2 

2 

i 

OF  BLOWS 

0.0 

12 

8 

11 

5 

37 

18.7 

40 

30 

32 

14 

115 

37.4 

54 

29 

29 

14 

126 

56.2 

46 

24 

27 

11 

107 

74.9 

55 

20 

19 

10 

105 

93.6 

18 

14 

13 

11 

56 

NORMALIZED 

LAYER  A  NORMALIZED  NUMBER  OF  BLOWS _ 

TOTAL  NUMBER 

1 

2 

2 

OF BLOWS 

0.0 

1.0 

1.0 

1.0 

1.0 

1.0 

18.7 

33 

3.8 

2.9 

2.8 

3.1 

37.4 

43 

3.6 

2.6 

2.8 

3.4 

56.2 

3.8 

3.0 

23 

2.2 

2.9 

74.9 

4.6 

2.5 

1.7 

2.0 

2.8 

93.6 

1.5 

1.8 

12 

2.2 

13 
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TABLE  4.  SHPB  TEST  RESULTS  FOR  OTTAWA  20-30  SAND  COMPACTED  MOIST 
AND  TESTED  MOIST. 


Container  Wall  Thickness  =  2.54  cm  (1.00  in) 

Dry  Density  =  1.715  g/cc  (107  pcO 

Void  Ratio  =  0.545 

Sample  Length  =  10.16  cm  (4.00  in) 

Testing  performed  by  Dr.  George  E.  Veyera  and  Blaise  J.  Fitzpatrick  of  the  University  of  Rhode  Island. 
Research  performed  at  AFESC/RDCM,  Tyndall  AFB,  Florida. 

Date:  Summer  1990 


total 


UNIT  WEIGHT 

WAVE  SPEED 

WAVE  SPEED 

TBANSMTSSION 

SSl 

g/gc 

(tUsec) 

rWsecl 

RATIO 

0.0 

0.00 

1.715 

1193 

363.6 

0.0168 

0.0 

0.00 

1.715 

1572 

479.1 

0.0298 

0.0 

0.00 

1.715 

1282 

390.7 

0.0199 

0.0 

0.00 

1.715 

1812 

552.3 

0.0511 

0.0 

0.00 

1.715 

1759 

536.1 

0.0517 

0.0 

0.00 

1.715 

1622 

494.4 

0.0330 

0.0 

0.00 

1.715 

1630 

496.8 

0.0270 

0.0 

0.00 

1.715 

1967 

599.5 

0.0434 

0.0 

0.00 

1.715 

2083 

634.9 

0.0670 

0.0 

0.00 

1.715 

1873 

570.9 

0.0457 

0.0 

0.00 

1.715 

2026 

617.5 

0.0749 

18.7 

3.85 

1.781 

1978 

602,9 

0.0350 

18.7 

3.85 

1.781 

1610 

490.7 

0.0429 

18.7 

3.85 

1.781 

2008 

612.0 

0.0635 

18.7 

3.85 

1.781 

1852 

564.5 

0.0534 

37.4 

7.70 

1.847 

1905 

580.6 

0.0653 

37.4 

7.70 

1247 

1967 

599.5 

0.0726 

37.4 

7.70 

1.847 

1817 

553.8 

0.0421 

56.2 

11J5 

1.913 

1754 

534.6 

0.0394 

562 

11.55 

1.913 

1921 

585.5 

0.0455 

562 

11.55 

1.913 

1852 

564.5 

0.0657 

74.9 

15.40 

1.979 

1822 

555.3 

0.0356 

74.9 

15.40 

1.979 

1873 

570.9 

0.0521 

93.6 

1925 

2.045 

1873 

570.9 

0.0313 

93.6 

1925 

2.045 

1727 

526.4 

0.0415 

93.6 

1925 

2.045 

1792 

546.2 

0.0307 

AVERAGE 

AVERAGE 

AVERAGE 

NORMALIZED 

NORMALIZED 

WAVE  SPEED 

WAVE  SPEED 

TRANSMISSION 

WAVE  SPEED 

TRANSMISSION 

&Sl 

fft/sec^ 

(m/sec) 

RATIO 

RATIO 

0.0 

1711 

521.5 

0.0418 

1.000 

1.000 

18.7 

1862 

567J 

0.0487 

1.088 

1.165 

37.4 

1896 

577.9 

0.0600 

1.108 

1.435 

56.2 

1842 

561.4 

0.0502 

1.077 

1.201 

74.9 

1848 

563.1 

0.0439 

1.080 

1.049 

93.6 

1797 

547.7 

0.0345 

1.050 

0.825 
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TABLE  5.  SHPB  TEST  RESULTS  FOR  OTTAWA  20-30  SAND  COMPACTED  MOIST 
AND  TESTED  DRY. 


Container  Wall  Thickness  =  2.54  cm  (1.00  in) 

Diy  E)ensity  =  1.715  g/cc  (107  pcQ 

Void  Ratio  =  0.545 

Sample  l.ength  =  10.16  cm  (4.00  in) 

Testing  performed  by  Dr.  George  R  Veyera  and  Blaise  J.  Fitzpatrick  of  the  University  of  Rhode  Island. 
Research  performed  at  AFESC/RDCM,  Tyndall  AFB,  Florida. 

Date:  Summer  1990 


total 


UNIT  WEIGHT 

WAVE  SPEED 

WAVE  SPEED 

TRANSMISSION 

Wc% 

g/cc 

(ft/sec> 

(m/sec) 

RATIO 

0.0 

0.00 

1.715 

1193 

363.6 

0.0168 

0.0 

0.00 

1.715 

1572 

479.1 

0.0298 

0.0 

0.00 

1.715 

1282 

390.7 

0.0199 

0.0 

0.00 

1.715 

1812 

552.3 

0.0511 

0.0 

0.00 

1.715 

1759 

536.1 

0.0517 

0.0 

0.00 

1.715 

1622 

494.4 

0.0330 

0.0 

0.00 

1.715 

1630 

496.8 

0.0270 

0.0 

0.00 

1.715 

1967 

599.5 

0.0434 

0.0 

0.00 

1.715 

2083 

634.9 

0.0670 

0.0 

0.00 

1.715 

1873 

570.9 

0.0457 

0.0 

0.00 

1.715 

2026 

617.5 

0.0749 

18,7 

3.85 

1.781 

2002 

610.2 

0.0577 

18.7 

3.85 

1.781 

2110 

643.1 

0.0576 

18.7 

3.85 

1.781 

2090 

637,0 

0.0743 

37.4 

7.70 

1.847 

1705 

519.7 

0.0536 

37.4 

7.70 

1.847 

2144 

653.5 

0.0780 

37.4 

7.70 

1.847 

1899 

578.8 

0.0519 

56.2 

11J5 

1.913 

1569 

478.2 

0.0482 

56.2 

1U5 

1.913 

1684 

5133 

0.0491 

56.2 

1155 

1.913 

2110 

643.1 

0.0605 

56^ 

1155 

1.913 

1990 

606.5 

0.0761 

74.9 

15.40 

1.979 

1495 

455.7 

0.0198 

74.9 

15.40 

1.979 

2008 

612.0 

0.0726 

74.9 

15.40 

1.979 

1889 

575.7 

0.0612 

74.9 

15.40 

1.979 

1910 

582.1 

0.0678 

AVERAGE 

AVERAGE 

AVERAGE 

NORMALIZED 

NORMALIZED 

WAVE  SPEED 

WAVE  SPEED 

TRANSMISSION 

WAVE  SPEED 

TRANSMISSION 

(fUsec\ 

(m/secl 

RATIO 

RATIO 

0.0 

1711 

521.5 

0.0418 

1.000 

1.000 

18.7 

2067 

630.0 

0.0632 

1.208 

1.512 

37.4 

1916 

584.0 

0.0612 

1.120 

1.464 

56.2 

1838 

560.2 

0.0585 

1.074 

1.400 

74.9 

1825 

5563 

0.0553 

1.067 

1.323 
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Figure  34.  Wave  Speed  as  a  Function  of  Saturation  for  Ottawa  20-30  Sand  Compacted 
Moist  and  Tested  Moist. 
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Figure  37.  Transmission  Ratio  as  a  Function  of  Saturation  for  Ottawa  20-30  Sand 
Compacted  Moist  and  Tested  Dry. 


from  zero  to  about  20  percent  saturation  and  gradually  decreases  thereafter.  Figure  37 
shows  the  stress  transmission  ratio  increasing  from  zero  to  about  20  percent  saturation  and 
gradually  decreasing  thereafter. 

D.  QUASI-STATIC  UNIAXIAL  COMPRESSION  TESTS 

A  total  of  13  quasi-static  uniaxial  compression  tests  were  performed  on  the  Ottawa 
20-30  sand  using  the  MTS  loading  device  at  several  saturations.  The  quasi-static  tests  were 
performed  to  obtain  information  about  the  stress-strain  response  and  stiffness  of  the 
compacted  unsaturated  soils  for  comparison  with  SHPB  trends.  Two  sets  of  tests  were 
performed:  the  first  used  four  compacted  layers  (as  for  SHPB  tests),  and  the  second  used 
one  compacted  layer.  Figures  38  through  49  show  the  stress-strain  curves  generated  from 
each  set  for  one  loading  and  unloading  cycle. 

1.  Four  Compacted  Layers 

Figures  38  through  43  show  stress- strain  curves  generated  at  different  saturations 
using  four  layers  compacted  to  a  constant  dry  density  of  1.715  g/cm^.  Each  specimen  was 
moist  during  testing.  In  all  cases  the  unloading  portion  of  the  stress-strain  curve  is  much 
steeper  than  the  loading  portion.  This  indicates  that  the  soil  is  stiffer  at  the  completion  of 
the  first  load-unload  cycle. 

2.  One  Compacted  LayCT 

Figures  44  through  49  show  stress-strain  curves  generated  at  different  saturations 
using  one  layer  compacted  to  a  constant  dry  density  of  1.715  g/cm^.  Each  specimen  was 
moist  during  testing.  In  all  cases  the  unloading  portion  of  the  stress-strain  curve  is  much 
Steeper  than  the  loading  portion.  This  indicates  that  the  soil  is  stiffer  at  the  completion  of 
the  first  load-unload  cycle. 

3 .  Constrained  Modulus 

Figure  50  shows  the  constrained  modulus  as  a  function  saturation  determined 
from  the  stress-strain  curves  generated  by  the  MTS  loading  device  (see  Figtues  38  through 
43).  The  constrained  modulus  was  determined  using  maximum  and  minimum  stress  levels 
for  the  first  loading  cycle  of  10,974  kPa  (1,592  psi)  and  zero  kPa  and  Equation  (18). 

Figure  50  shows  that  the  constrained  modulus  sharply  increases  from  zero  to  about  20 
percent  saturation  and  has  a  slightly  increases  from  about  20  percent  to  about  40 
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Figure  50.  Constrained  Modulus  as  a  Function  of  Saturation  for  Quasi-Static 
Testing  of  Ottawa  20-30  Sand  Compacted  Moist  and  Tested  Moist. 


percent  saturation.  The  constrained  modulus  then  decreases  from  about  40  percent  to  about 
60  percent  saturation  and  slightly  increases  from  about  60  percent  to  about  75  percent 
saturation,  sharply  decreasing  thereafter. 

E.  MICROSTRUCTURAL  ANALYSIS 

1 .  Axial  Ratio 

Figure  5 1  shows  the  results  of  the  axial  ratio  particle  analysis  and  is  a  composite 
plot  of  more  than  2500  individual  particles  which  were  manually  analyzed  using  the  method 
developed  by  Campbell  (1985).  The  results  for  the  Ottawa  20-30  sand  indicate  that  the 
average  axial  ratio,  navg ,  is  0.674  and  the  standard  deviation  is  0.106.  This  indicates  that 
the  particles  are  approaching  a  spherical  shape.  Axial  ratio  values  ranged  from  a  minimum 
of  0.253  to  a  maximum  of  0.968.  Based  on  the  axial  ratio  data  obtained,  the  Ottawa  20-30 
sand  would  be  classifted  as  consisting  of  well  rounded  particles  according  to  the  criteria 
shown  in  Table  1. 

2.  Rose  Diagrams 

Particle  long-axis  orientations  are  represented  by  rose  diagrams  (angular 
frequency  histograms).  A  total  of  70  rose  diagrams  were  generated  from  azimuth  data 
obtained  by  digitizing  a  number  of  two-dimensional  particle  images  using  two  main 
categories  (see  Appendix  A  Figures  A.l  through  A.70).  The  first  set  of  rose  diagrams 
includes  specimens  compacted  using  the  standard  Proctor  hammer  only  (NO  SHPB)  and 
the  second  set  of  rose  diagrams  include  specimens  compacted  using  the  standard  Proctor 
hammer  gni  then  tested  in  the  SHPB  device  (SHPB). 

The  two  sets  of  rose  diagrams  (NO  SHPB  and  SHPB)  were  further  subdivided 
into;  (a)  horizontal,  vertical,  and  semicircular  rose  diagrams  (a  total  of  30  rose  diagrams), 
and  (b)  Layers  1,  2,  3,  and  4  in  the  vertical  plane  (a  total  of  40  rose  diagrams).  Rose 
diagrams  for  each  of  the  four  layers  were  developed  using  grids  h,  i,  j,  k,  d,  1,  m,  n,  and  o 
from  Section  A-A  in  the  vertical  plane  (see  Figure  27).  Each  layer  consisted  of  two  and 
one-quarter  grids,  (eg.  Layer  1  used  grids  o,  n,  and  one-fourth  of  grid  m).  The  complete 
set  of  rose  diagrams  can  be  found  in  Appendix  A. 

Statistics  provided  within  the  "ROSY"  software  provide  the  following  data:  the 
number  of  particles  measured,  the  maximum  percentage  of  particles  that  fall  within  a  given 
class  interval,  the  mean  and  standard  deviation  of  this  percentage,  vector  mean,  confidence 
angle,  vector  magnitude,  and  Rayleigh  test  of  significance.  A  summary  of  the  statistical 
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data  generated  by  ROSY  for  each  rose  diagram  has  been  compiled  in  Tables  6  and  7.  A 
description  of  each  heading  used  in  rose  diagrams  and  in  Tables  6  and  7  is  given  in  Table 
8.  A  typical  rose  diagram  is  presented  in  Figure  52.  This  figure  shows  how  to  interpret 
the  rose  diagram  statistics  output  by  Rosy.  Comparisons  of  the  results  from  the  statistical 
analysis  of  all  70  rose  diagrams  will  be  discussed  in  SECTION  V. 

Table  6  shows  the  final  results  of  the  number  of  particles  measured  for  each  rose 
diagram,  N,  maximum  percentage,  vector  mean,  contidence  angle,  vector  magnitude,  and 
Rayleigh  test  of  significance  as  a  function  of  saturation  for  the  horizontal,  vertical,  and 
semicircular  sections,  and  Layers  1,  2,  3,  and  4  in  the  vertical  plane  on  Section  A- A  (see 
Figure  27),  for  NO  SHPB  specimens.  Table  7  shows  the  final  results  as  a  function  of 
saturation  for  the  horizontal,  vertical,  and  semicircular  sections,  and  Layers  1, 2,  3,  and  4 
in  the  vertical  plane  on  Section  A-A  (see  Figure  27),  for  SHPB  specimens.  All  data  were 
collected  from  rose  diagrams  generated  by  the  ROSY  program. 

The  Rosy  data  given  in  Tables  6  and  7  show  that  the  semicircular  plane  has  the 
highest  Rayleigh  values  and  lowest  vector  magnitude  values.  This  indicates  that  there  is  no 
preferred  particle  orientation  in  any  semicircular  plane.  The  number  of  particles  measured 
in  the  semicircular  plane  (150  to  3 10)  are  much  lower  than  those  measured  in  the  horizontal 
or  vertical  planes  (280  to  760).  Only  4  grids  (p  through  s)  were  available  for  measurement 
in  the  semicircular  plane,  whereas  the  horizontal  plane  had  7  grids  (a  through  g),  the 
vertical  plane  had  and  9  grids  (h  through  o,  including  d).  Oda  (1972a)  and  Campbell 
(1985)  stated  that  a  minimum  of  200  to  300  particles  need  to  be  measured  for  a  particle 
orientation  analysis.  Therefore,  due  to  the  limited  amount  of  data  for  the  semicircular 
plane,  those  results  will  not  be  analyzed  in  SECTION  V.  However,  data  for  the  vertical 
and  horizontal  planes  show  Rayleigh  values  less  than  0.05  and  high  vector  magnitude 
values  which  indicates  there  is  a  preferred  particle  orientation  in  these  planes.  Therefore, 
the  microstructural  analysis  presented  in  SECTION  V  will  focus  on  results  for  the  vertical 
and  horizontal  planes. 

3.  Particle  Long- Axis  Orientation 

Figure  53  shows  vector  magnitude  as  a  function  of  saturation  in  the  horizontal 
(upper  portion  of  Layer  2  and  lower  portion  of  Layer  3),  vertical  (Layers  1  through  4),  and 
semicircular  planes  for  NO  SHPB  specimens  (see  Figure  25).  In  the  horizontal  and  vertical 
planes  the  vector  magnitude  tends  to  be  lower  at  intermediate  saturations  and  higher  at  zero 
p>ercent  and  about  75  percent  saturation.  In  the  semicircular  plane,  the  vector  magnitude 
tends  to  be  higher  at  intermediate  saturations  and  lower  at  0  percent  and  80  percent 
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TABLE  6.  STATISTICAL  RESULTS  FOR  MICROSTRUCTURAL  ANALYSIS 
OF  NO  SHPB  SPECIMENS, 


sa- 

N 

Mat.  % 

V.M. 

Conf  Angle 

R.  May. 

Baylcleli 

Section 

0.0 

615 

8.6 

277.8 

1734 

0.182 

0.0000 

horizontal 

18.7 

487 

9.7 

280.8 

23.68 

0.152 

0.0000 

horizontal 

37.4 

618 

95 

274.4 

17.73 

0.177 

0.0000 

horizontal 

56.2 

636 

9.6 

883 

22.62 

0.137 

0.0000 

horizontal 

74.9 

672 

9J 

270.0 

12.69 

0.236 

0.0000 

horizontal 

S% 

Max.  % 

VJVf. 

Conf  Angte 

R.  Mag. 

Rayleigh 

Section 

0.0 

607 

9.1 

275.0 

1332 

0340 

0.0000 

vertical 

18.7 

287 

10.1 

85.6 

5236 

0,089 

0.1028 

vertical 

37.4 

626 

10.1 

276.8 

1739 

0,182 

0.0000 

vertical 

56.2 

729 

11.7 

87.9 

1331 

0.220 

0.0000 

vertical 

74.9 

764 

9.8 

89.1 

12.44 

0.226 

0.0000 

vertical 

S% 

a 

Max.  % 

VJVf. 

Conf  Anghi 

R-  Mag- 

Rayleigh 

Section 

0.0 

245 

9.4 

79.4 

227.13 

0.025 

0.8587 

semi-circular 

18.7 

250 

9.2 

272.9 

6339 

0.078 

0,2146 

semi’circular 

37.4 

310 

8.7 

253 

87.98 

0.052 

0.4256 

scmi-ciicular 

56.2 

265 

9.8 

291.9 

53.65 

0.091 

0,1089 

semi-circular 

74.9 

237 

8.4 

383 

85.43 

0.061 

0.4183 

semi-circular 

SSl 

Max.  % 

Conf  Angle 

R-May- 

Rayleigh 

Laver  # 

0.0 

69 

21.7 

290.0 

2532 

0363 

0.0001 

1 

18.7 

147 

13.6 

2793 

1832 

0345 

0.0000 

1 

37.4 

160 

11.2 

783 

2134 

0385 

0.0000 

1 

56.2 

164 

12.8 

82.8 

22,66 

0.267 

0.0000 

1 

74.9 

111 

10.7 

86.4 

31.15 

0.190 

0.0017 

1 

S% 

Max.  % 

VJVf. 

Cggf  Angle 

R.  Mag. 

Rayleigh 

Layer# 

0.0 

197 

11.2 

2703 

2934 

0.190 

0.0008 

2 

18.7 

151 

9.9 

733 

39.64 

0.164 

0.0173 

2 

37.4 

179 

123 

2763 

27.80 

0312 

0.0003 

2 

56.2 

181 

13.8 

883 

19.97 

0388 

0.0000 

2 

74.9 

189 

10.1 

81.7 

2434 

0.234 

0.0000 

2 

S% 

Max-  % 

VJVf- 

Conf  Angle 

R-  Ma^- 

Rayleigh 

Laver# 

0.0 

213 

8.9 

89.1 

33.28 

0.165 

0.0030 

3 

18.7 

90 

8.9 

2993 

60.07 

0.137 

0.1833 

3 

37.4 

186 

11.8 

296.4 

31.96 

0.181 

0.0022 

3 

56.2 

202 

9.9 

89.7 

40.17 

0.137 

0.0228 

3 

74.9 

218 

12.4 

275.0 

1934 

0.270 

0.0000 

3 

S% 

£[ 

Max.  % 

VJVf. 

Cggf  Angle 

R.  Mag. 

Rayleigh 

Layer  ff 

0.0 

128 

123 

2753 

16.26 

0.414 

0.0000 

4 

18.7 

15 

20.0 

336.9 

84.43 

0.238 

0.4261 

4 

37.4 

101 

12.9 

2903 

53.17 

0.146 

0.1176 

4 

56.2 

181 

12.7 

271.4 

28.07 

0306 

0.0005 

4 

74.9 

179 

11.2 

270.2 

27.85 

0.211 

0.0003 

4 
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TABLE  7.  STATISTICAL  RESULTS  FOR  MICROSTRUCTURAL  ANALYSIS 
OF  SHPB  SPECIMENS. 


s% 

a 

Max,  % 

VJVf. 

Conf  Angle 

R.  May. 

Raxleisb 

Section 

0.0 

420 

143 

279.4 

9.05 

0.407 

0.0000 

horizontal 

18.7 

742 

12.4 

274.0 

1238 

0334 

0.0000 

horizontal 

37.4 

556 

10.4 

271.6 

13.79 

0341 

0.0000 

horizontal 

56.2 

554 

103 

2733 

16.62 

0302 

0.0000 

horizontal 

74.9 

510 

10.4 

275.8 

2039 

0.174 

0.0000 

horizontal 

li 

Max.  % 

V.M, 

Conf  Angle 

R.  Mag. 

Rayleigh 

Section 

0.0 

544 

13.6 

279,6 

9.43 

0348 

0.0000 

vertical 

18.7 

738 

10.8 

88.8 

10.70 

0366 

0.0000 

vertical 

37.4 

455 

143 

278.6 

1135 

0324 

0.0000 

vertical 

56.2 

632 

10.4 

272.1 

13.67 

0326 

0.0000 

vertical 

74.9 

607 

11.9 

2813 

12.07 

0365 

0.0000 

vertical 

S% 

Max,  % 

VJVf. 

Conf  Angle 

R.  May. 

Ravlclyh 

Section 

0.0 

159 

10.7 

312.0 

15838 

0.039 

0.7813 

semi-circular 

18.7 

320 

9.7 

270.7 

40.76 

0.107 

0.0253 

semi-circular 

37.4 

218 

113 

278.9 

13539 

0.039 

0.7141 

semi-circular 

56.2 

229 

92 

2913 

52.73 

0.099 

0.1083 

semi-circular 

74.9 

191 

8.4 

293.7 

7238 

0.079 

0.3060 

semi-circular 

S% 

a 

Max,  % 

VJVf. 

Conf  Angle 

R,  Mag. 

Ravlclgh 

Laver# 

0.0 

113 

143 

275.4 

2131 

0340 

0.0000 

1 

18.7 

182 

11.0 

87.1 

21.49 

0367 

0.0000 

1 

37.4 

87 

13.8 

2763 

1833 

0.444 

0.0000 

1 

56.2 

239 

11.7 

833 

1733 

0390 

0.0000 

1 

74.9 

145 

12.4 

883 

29.73 

0318 

0.0010 

1 

S% 

a 

Max.  % 

VJVf. 

RkMggft 

Ravleigh 

Laver  # 

0.0 

129 

133 

274.9 

19.93 

0348 

0.0000 

2 

18.7 

205 

133 

271.4 

16.96 

0317 

0.0000 

2 

37.4 

114 

16.7 

209.4 

1739 

0.406 

0.0000 

2 

56.2 

139 

13.7 

84.8 

3738 

0.177 

0.0126 

2 

74.9 

170 

12.9 

2753 

15.06 

0386 

0.0000 

2 

S% 

a 

Max,  % 

VJVf. 

Conf  Angle 

R.  May. 

Rayleigh 

Laver  # 

0.0 

174 

14.9 

285.1 

1632 

0357 

0.0000 

3 

18.7 

179 

11.7 

274.9 

2337 

0350 

0.0000 

3 

37.4 

137 

10.9 

88.4 

29.12 

0330 

0.0007 

3 

56.2 

111 

103 

284.1 

22.69 

0357 

0.0000 

3 

74.9 

145 

153 

281.8 

1833 

0346 

0,0000 

3 

S% 

M 

Max.  % 

VJVf. 

Conf  Angle 

R.  May. 

Ravleigh 

Laver# 

0.0 

128 

123 

2803 

18.74 

0363 

0.0000 

4 

18.7 

172 

13.4 

80.1 

24.81 

0341 

0.0000 

4 

37.4 

117 

18.8 

275.4 

22.42 

0317 

0.0000 

4 

56.2 

161 

113 

277.7 

26.92 

0329 

0.0000 

4 

74.9 

147 

93 

307.6 

31.18 

0305 

0.0020 

4 

90 


saturation.  Figure  54  shows  vector  magnitude  as  a  function  of  saturation  in  the  horizontal, 
vertical,  and  semicircular  planes  for  SHPB  specimens.  In  the  horizontal  and  vertical  planes 
the  vector  magnitude  tends  to  be  lower  at  all  saturations  and  higher  at  zero  percent 
saturation.  In  the  semicircular  plane  the  vector  magnitude  tends  to  be  higher  at  intermediate 
saturations  and  lower  at  zero  percent  and  about  75  percent  saturation.  SHPB  specimens 
should  have  higher  vector  magnitudes  since  they  were  tested  in  the  split-Hopkinson 
pressure  bar  device  after  compaction. 

TABLE  8.  DESCRIPTION  OF  STATISTICS  USED  IN  ROSY  PROGRAM. 


-  The  number  of  particles  measured  for  each  rose  diagram. 

Class  Interval  -  The  width  of  each  petal  on  the  rose  diagram  (10-degree  intervals  for 
this  investigation). 

Max.  %  -  The  maximum  percentage  of  particles  that  occur  on  a  given  cell  (or  petal  of 
the  rose  diagram). 

Mean  %  -  Mean  percentage  of  the  average  percent  of  all  data. 

Standard  Deviation  -  Standard  deviation  of  the  mean  percentage. 

Vector  Mean  -  the  average  azimuth  of  all  particles  measured. 

Conf.  Angle  -  Confidence  angle.  For  a  confidence  interval  of  95  percent,  this  shows 
that  95  percent  of  the  particles  measured  lie  within  some  angle  ±  of  the  vector  mean. 

R.  Magnitude.  -  Vector  magnitude,  values  range  from  0.0  to  1.0,  data  sets  which 
have  large  dispersion  will  have  R.  Mag.  values  near  0.0,  data  sets  with  very  small 
dispersions  will  have  R.  Mag.  values  near  1.0. 

Rayleigh  -  Rayleigh  test  of  significance,  if  a  value  of  0.05  or  greater  occurs,  then  the 
data  is  said  to  be  no  good  and  no  preferred  particle  orientation  exists.  If  a  value  <  0.05 
is  encountered,  then  a  preferred  particle  orientation  exists  within  the  sample. 


Figures  55  through  58  compare  vector  magnitude  as  a  function  of  saturation  for 
each  individual  compacted  layer  (NO  SHPB  or  SHPB).  In  nearly  all  cases,  the  vector 
magnitude  is  higher  for  SHPB  specimens  as  compared  to  NO  SHPB  specimens.  For 
SHPB  specimens  the  vector  magnitude  tends  to  be  higher  at  zero  percent  saturation  and 
lower  at  all  other  saturations.  For  NO  SHPB  specimens  a  general  trend  occurs  throughout 
most  of  the  four  layers  in  that  the  vector  magnitude  tends  to  be  lower  at  intermediate 
saturations  and  higher  at  zero  percent  and  about  75  percent  saturation. 
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Magnitude  as  a  Function  of  Saturation  for  Layer  #1 
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Figure  56.  Vector  Magnitude  as  a  Function  of  Saturation  for  Layer  #2 


Magnitude  as  a  Function  of  Saturation  for  Layer  #3. 


SECTION  V 


ANALYSIS  OF  RESULTS 
A.  COMPACTION  ENERGY 

A  summary  of  the  compaction  energy  results  using  the  standard  Proctor  hammer  is 
presented  in  Figure  59,  which  shows  normalized  total  compactive  energy  as  a  function  of 
saturation.  As  can  be  seen,  the  amount  of  energy  required  to  compact  the  Ottawa  20-30 
sand  to  a  constant  dry  density  of  1.715  g/cm^  varied  with  saturation.  This  trend  has  been 
observed  by  several  other  researchers  for  Ottawa  20-30  sand  and  other  sands  (Ross  et 
al.,1986;  Pierce,  1989;  Veyera,  1989).  The  results  from  Figure  59  show  that  the  amount 
of  required  compactive  energy  sharply  increases  from  zero  to  about  20  percent  saturation, 
then  slightly  increases  from  about  20  percent  to  about  40  percent  saturation  and  decreases 
thereafter. 

There  appears  to  be  an  increase  in  soil  stiffness  at  intermediate  saturations  which  is 
most  likely  due  to  capillary  tension  forces  present  during  compaction.  The  particles  are 
most  resistant  to  rearrangement  at  saturations  between  about  20  percent  and  40  percent. 
The  higher  compaction  energy  required  at  these  saturations  may  be  due  to  the  particles 
resisting  rearrangement,  thus  requiring  more  energy  for  compaction.  At  higher  saturations 
(greater  than  about  60  percent),  the  moisture  may  somewhat  "lubricate"  the  soil  particles  as 
surface  tension  forces  are  reduced.  This  would  allow  the  particles  to  slide  past  one  another 
during  compaction  with  less  resistance,  therefore  decreasing  the  amount  of  compactive 
energy  necessary. 

Figure  60  shows  compactive  energy  results  for  each  individual  2.54  cm  (1.00  inches) 
thick  layer.  The  data  have  been  normalized  to  the  average  compactive  energy  required  for 
the  first  layer  at  zero  saturation.  As  shown,  more  compactive  energy  was  required  to 
compact  the  first  layer  at  all  saturations  as  compared  with  Layers  2,  3,  and  4.  This  may  be 
explained  as  follows:  prior  to  the  first  blow  for  Layer  1,  the  soil  is  in  its  loosest  state  from 
initial  placement  After  the  first  blow  is  applied,  the  density  of  the  soil  increases  because 
during  the  first  blow  the  particles  are  forced  to  rearrange  themselves  into  a  much  denser 
packing.  A  compressional  stress  wave  generated  by  the  standard  Proctor  hammer  is 
transmitted  through  the  soil  and  into  the  steel  pedestal  below.  As  the  energy  from  this 
wave  propagates  through  the  soil,  it  compacts  and  increases  the  density  of  the  soil.  When 
the  compressional  stress  wave  encounters  the  steel  pedestal  it  is  reflected  back  as  a  tensile 
wave  which  may  loosen  the  layer  being  compacted. 
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Each  of  the  four  layers  show  the  same  general  trend  as  that  found  in  Figure  59,  i.e., 
the  amount  of  energy  required  to  compact  the  soil  to  a  constant  dry  density  sharply 
increases  from  zero  to  about  20  percent  saturation  and  decreases  thereafter.  The  trend  is 
similar  yet  less  pronounced  in  Layers  2,  3,  and  4  than  it  is  in  Layer  1  (Figure  60). 
Different  stiffnesses  and  densities  exist  between  the  compacted  soil  and  steel  pedestal. 
Differences  in  compaction  energy  between  the  four  layers  may  be  due  the  magnitude  of  the 
reflected  tensile  wave  being  reduced  as  the  specimen  is  compacted  to  final  height,  creating 
less  disturbance  in  Layers  2, 3,  and  4.  Also,  the  amount  of  energy  needed  to  compact  the 
soil  to  the  required  dry  density  decreases  with  each  layer. 

Figure  61  shows  cumulative  compactive  energy  results  (total  blows  applied  to  each 
layer  including  those  of  succeeding  layers)  for  each  individual  2.54  cm  (1.00  inches)  thick 
layer.  The  cumulative  energy  applied  to  Layer  1  represents  the  total  amount  of  compaction 
energy  (see  Figure  32)  applied  to  the  sample.  Each  of  the  four  curves  in  Figure  61  show 
the  same  general  trend,  i.e.,  less  energy  was  required  to  compact  each  layer  to  a  constant 
dry  density  at  zero  percent  saturation  and  more  energy  was  required  to  compact  each  layer 
to  a  constant  dry  density  at  intermediate  saturations  (20  percent  to  60  percent).  Similar 
behavior  is  expected  between  the  four  layers  because  the  same  soil  was  compacted  to  the 
same  constant  dry  density.  Thus,  it  seems  that  the  effect  of  moisture  (capillary  tension 
forces)  during  compaction  on  each  individual  layer  may  be  significant 

B .  SPLIT-HOPKINSON  PRESSURE  BAR  RESULTS 

1 .  Compressional  Wave  Speed  and  Transmission  Ratio 

The  split-Hopkinson  pressure  bar  was  used  to  measure  compressional  wave 
speed  and  stress  transmission  ratio  as  a  function  of  saturation  for  Ottawa  20-30  sand. 
Figures  62  and  63  show  results  for  specimens  compacted  moist  and  tested  moist  (moist- 
moist),  and  Figures  64  and  65  show  results  for  specimens  compacted  moist  and  tested  dry 
(moist-dry).  Specimens  were  tested  under  moist-moist  and  moist-dry  conditions  to 
determine  if  the  presence  of  water  during  SHPB  testing  significantly  affected  the 
compressional  wave  speed  and  stress  transmission  ratio,  or  if  the  observed  behavior  was 
more  strongly  influenced  by  microstructure  alone. 

a.  Moist-Moist  Results 

Results  obtained  from  moist-moist  tests  in  this  investigation  are  shown  in 
Figures  62  and  63  and  indicate  trends  similar  to  those  observed  in  previous  research 
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Figure  61.  Cummulative  Blow  Count  as  a  Function  of  Saturation  for  Ottawa  20-30 
Sand  for  Layers  1  to  4. 
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Figure  62.  Normalized  Average  Wave  Speed  as  a  Function  of  Saturation  for 
Ottawa  20-30  Sand  Compacted  Moist  and  Tested  Moist. 
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Figure  63.  Normalized  Average  Transmission  Ratio  as  a  Function  of 
Saturation  for  Ottawa  20-30  Sand  Compacted  Moist  and  Tested  Moist. 
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Figure  64.  Normalized  AverageWave  Speed  as  a  Function  of  Saturation  for  Ottawa  20-30 
Sand  Compacted  Moist  and  Tested  Dry. 
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Figure  65.  Normalized  Average  Transmission  Ratio  as  a  Function  of  Saturation  for 
Ottawa  20-30  Sand  Compacted  Moist  and  Tested  Dry. 


(Ross,  1989)  using  the  same  SHPB.  The  data  have  been  normalized  to  the  compressional 
wave  speed  and  stress  transmission  ratio  at  zero  saturation  (dry  soil).  Figures  62  and  63 
show  that  at  intermediate  saturations  (20  percent  to  75  percent),  compressional  wave  speed 
and  stress  transmission  ratio  are  approximately  1.11  and  1.44  times  greater  in  magnitude, 
respectively,  than  for  dry  soil.  As  saturations  increase  to  about  75  percent,  both  the 
compressional  wave  speed  and  stress  transmission  ratio  are  roughly  1.08  and  1.05  times 
greater  in  magnitude,  respectively,  than  for  dry  soil.  Beyond  75  percent  saturation, 
compressional  wave  speed  and  stress  transmission  ratio  decrease,  however  they  will 
increase  as  the  soil  approaches  a  saturated  condition  (i.e.,  S  =  100  percent).  Saturated  soil 
will  have  a  slightly  greater  value  of  compressional  wave  speed  than  the  theoretical 
compressional  wave  speed  in  water  due  to  a  higher  stiffiiess  and  total  density  of  the  soil- 
water  mixture.  The  theoretical  compressional  wave  speed  in  water  is  1433  meters/sec 
(4700  feet/sec)  at  atmospheric  pressure  and  at  a  temperature  of  20°  C  (68°  F). 

The  results  for  compressional  wave  speed  and  stress  transmission  ratio  as  a 
function  of  saturation  obtained  in  this  study  are  similar  to  those  obtained  by  other 
researchers  who  also  performed  SHPB  tests  on  the  Ottawa  20-30  sand  compacted  moist 
and  tested  moist  (Veyera,  1989;  Ross,  1989;  Felice  et  al.,  1987;  and  Ross  et  al.,  1986). 
At  intermediate  saturations  (20  percent  to  75  percent)  these  values  were  greater  than  those  at 
zero  percent  saturation,  and  may  be  due  to  an  increase  in  soil  stiffness  from  capillary  forces 
(see  SECTION  V).  This  applies  to  SHPB  specimens  which  were  compacted  moist  and 
tested  moist  and  to  those  compacted  moist  and  tested  dry. 

Results  of  compressional  wave  speed  and  stress  transmission  ratio  as  a 
function  of  saturation  in  this  investigation  differed  firom  results  presented  by  Pierce  (1989) 
who  compacted  specimens  dix  saturated  them  to  the  desired  level  by  means  of  a 
pressure  plate  apparatus,  thus  the  influence  of  moisture  during  compaction  was  not  present 
Pierce  did  not  observe  any  changes  in  transmission  ratio  or  wave  speed  for  specimens 
prepared  as  stated.  His  results  were  not  replicated  in  this  study,  and  this  can  be  attributed 
to  different  sample  preparation  techniques.  In  this  investigation  specimens  were  compacted 
moist,  thus  the  influence  of  moisture  during  compaction  is  present  Differences  between 
results  obtained  by  Pierce  (1989),  Ross  (1986),  Veyera  (1989),  and  the  present 
investigation  show  that  when  moisture  is  present  during  compaction,  the  compressional 
wave  speed,  stress  transmission  ratio,  and  constrained  modulus  are  greatly  influenced. 
Therefore,  it  is  evident  that  the  specimen  preparation  method  can  significantly  affect  test 
results. 


108 


b.  Moist-Dry  Results 

Results  obtained  from  moist-dry  tests  in  Figures  64  and  65  show  that  at  a 
saturation  of  about  20  percent,  compressional  wave  speed  and  stress  transmission  ratio  are 
approximately  1.21  and  1.52  times  greater  in  magnitude,  respectively,  than  for  dry  soil. 
The  data  have  been  normalized  to  the  compressional  wave  speed  and  stress  transmission 
ratio  at  zero  saturation.  As  the  saturation  increases  beyond  about  20  percent  to  about  75 
percent  in  both  Figures  64  and  65,  the  compressional  wave  speed  and  stress  transmission 
ratio  steadily  decrease  to  approximately  1.07  and  1.34  times  the  magnitude  of  that  for  dry 
soil,  respectively. 

2 .  Comparison  Between  Moist-Moist  and  Moist-Dry  Results 

A  comparison  between  the  results  obtained  from  moist-moist  and  moist-dry  tests 
is  shown  in  Figures  66  and  67.  Both  test  series  produced  similar  results  in  that  at 
intermediate  saturations  (20  percent  to  60  percent),  compressional  wave  speed  and  stress 
transmission  ratio  had  greater  values  than  at  zero  saturation.  This  further  reinforces  the 
hypothesis  that  soil  microstructure  developed  during  compaction  may  be  the  primary 
phenomenon  affecting  the  ccxnpressional  wave  speed  and  stress  transmission  ratio. 

The  results  obtained  in  moist-dry  testing  conditions  showed  higher  values  of 
compressional  wave  speed  and  stress  transmission  ratio  at  all  compacted  saturations  over 
moist-moist  results  (Figures  66  and  67).  This  suggests  that  moisture  present  during  SHPB 
testing  also  affects  the  compressional  wave  speed  and  stress  transmission  ratio.  Moist- 
moist  tests  were  conducted  on  a  three  phase  material  exists,  i.e.,  soil,  water,  and  air, 
whereas  moist-dry  tests  were  conducted  on  a  two  phase  material,  i.e.,  soil  and  air. 
Therefore,  if  water  is  present,  frictional  resistance  between  particles  may  be  lower  and  the 
compressional  wave  speed  transmission  ratio  will  have  lower  values  during  SHPB  testing. 

Capillary  tension  forces  are  present  during  moist-moist  tests,  however  their  effect 
on  compressional  wave  speed  and  stress  transmission  ratio  is  negligible.  Pierce  (1989) 
studied  capillary  stresses  in  unsaturated  Ottawa  20-30  sand  and  found  that  they  are 
extremely  small  (on  the  order  of  7  kPa  (1  psi))  as  compared  to  stress  levels  encountered  in 
split-Hopkinson  pressure  bar  testing  (on  the  order  of  15,000  kPa  (2,175  psi)).  Therefore, 
they  have  little  effect  on  dynamic  testing.  However,  the  influence  of  these  capillary  tension 
forces  in  the  development  of  soil  microstructure  during  compaction  appears  to  be  a 
significant  factor  which  affects  both  the  static  and  dynamic  soil  behavior. 
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C .  QUASI-STATIC  COMPRESSION  TESTS  AND  CONSTRAINED  MODULUS 


Figures  68  and  69  show  comparison  stress-strain  curves  generated  from  the  MTS 
device  using  four  compacted  layers  and  one  compacted  layer,  respectively,  at  each 
saturation.  The  stress-strain  curve  at  zero  percent  saturation  is  compared  to  stress-strain 
curves  at  samrations  of  18.7  percent,  37.4  percent,  56.2  percent,  74.9  percent,  and  93.6 
percent  Dry  soil  values  (S  =  0  percent)  are  used  as  a  basis  for  comparison  in  normalizing 
the  results. 

Figure  68  shows  that  the  slope  of  the  stress-strain  curve  for  zero  percent  saturation 
was  steeper  than  93.6  percent  whereas  the  slopes  of  the  stress-strain  curves  for  saturations 
of  18.7  percent,  37.4  percent,  56.2  percent  and  74.9  percent  were  generally  steeper  than 
zero  percent  saturation,  for  specimens  with  four  compacted  layers.  Figure  69  shows  that 
the  slope  of  the  stress-strain  curve  for  zero  percent  was  steeper  than  37.4  percent  and  93.6 
percent,  whereas  the  slope  of  the  stress-strain  curve  for  saturations  of  18.7  percent,  56.2 
percent  and  74.9  percent  was  steeper  than  zero  percent  saturation,  for  specimens  with  one 
compacted  layer.  The  general  trend  in  Figures  68  and  69  is  an  increase  in  overall  stiffness 
with  increasing  saturation. 

In  Figures  68  and  69,  the  stress-strain  curves  which  have  the  steepest  slope  indicate  a 
stiffer  soil.  In  comparing  stress- strain  results  for  four  layers  and  one  layer  it  can  be  seen 
that  specimens  with  one  layer  were  much  stiffer  than  those  with  four  layers.  This  may  be 
due  to  the  initial  length  of  the  specimen.  For  one  layer,  the  initial  specimen  length  was 
2.54  cm  (1.0  inches).  Thus  the  amount  of  strain  (deformation)  produced  is  smaller  and  the 
specimen  behaves  as  a  stiffer  material.  Both  specimens  were  subjected  to  the  same  uniaxial 
compressive  load.  In  general,  as  saturations  are  increased  to  about  75  percent,  the  soil 
becomes  stiffer  as  compared  with  zero  percent  saturation.  As  saturations  increase  past  75 
percent,  the  soil  becomes  less  stiff.  This  compares  well  with  SHPB  test  results  for 
compressional  wave  speed  and  stress  transmission  ratio  in  that  a  stiffer  material  is  present 
to  saturations  of  about  75  percent 

Figure  70  shows  that  the  constrained  modulus  (normalized  to  results  at  S  =  0  percent) 
was  found  to  be  approximately  1.08  times  greater  at  saturations  to  about  40  percent  and 
was  approximately  1.03  times  greater  at  saturations  between  about  60  percent  and  75 
percent  A  further  increase  in  saturation  showed  that  the  constrained  modulus  was 
approximately  0.97  times  the  constrained  modulus  at  zero  percent  saturation.  This  ttend 
(higher  values  at  intermediate  saturations)  was  also  observed  in  the  compaction  data  (Figure 
59),  compressional  wave  speed  data  (Figure  62),  and  stress  transmission  ratio  data 
(Figure63)  as  a  function  of  saturation.  The  behavior  may  be  due  to  the  resistance  of 
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Figure  70.  Normalized  Constained  Modulus  as  a  Function  of  Saturationfor 
Ottawa  20-30  Sand  Compacted  Moist  and  Tested  Moist. 


particles  to  relative  movements  at  intermediate  saturations.  By  applying  a  quasi-static 
loading,  the  soil  is  compressed  into  a  packing  greater  than  that  from  the  initial  compaction 
procedure.  As  in  the  impact  compaction  process  using  the  standard  Proctor  hammer, 
capillary  forces  probably  tend  to  resist  particle  rearrangements  at  lower  saturations  in  these 
tests.  Capillary  forces  decrease  with  increasing  saturation,  thereby  reducing  the  frictional 
resistance  between  soil  particles  which  may  allow  particles  to  move  relative  to  one  another 
with  greater  ease. 

Figure  71  shows  the  normalized  constrained  modulus  and  normalized  total  blow 
count  (normalized  to  results  at  S  =  0  percent)  as  a  function  of  saturation.  The  same  trend  is 
noticeable  for  constrained  modulus  and  total  blow  count,  i.e.,  high  values  are  obtained  at 
intermediate  saturations.  The  constrained  modulus  and  total  blow  coimt  can  be  related  to 
soil  stiffness.  A  stiffer  material  will  produce  a  higher  constrained  modulus  and  more 
energy  will  be  required  to  compact  a  stiffer  material  to  a  constant  dry  density.  Therefore, 
the  relationships  shown  in  Figure  71  are  funher  evidence  that  moisture  present  during 
compaction  affects  soil  stiffness  and,  in  turn,  the  constrained  modulus  and  required 
compaction  energy. 

D.  AXIAL  RATIO 

More  than  2,500  particles  were  measured  to  determine  the  axial  ratio  of  the  Ottawa 
20-30  sand  used  in  this  investigation.  The  average  axial  ratio,  navg  ,  was  found  to  be 
0.674  with  a  standard  deviation  of  0.106.  Axial  ratio  values  ranged  from  0.253  to  0.968. 
Ottawa  20-30  sand  is  classified  as  a  well  rounded  sand  according  to  Table  1. 

Oda  (1972a)  proposed  that  the  minimum  axial  ratio  required  for  preferred  particle 
orientation  is  0.7.  An  important  note  is  that  axial  ratio  is  a  "shape  parameter"  and  is  not  a 
function  of  saturation  within  a  soil.  However,  particle  shape  does  affect  the  probability  of 
obtaining  a  preferred  particle  orientation.  With  an  average  axial  ratio  of  0.674  the  particles 
are  well-rounded  and  there  is  a  low  resistance  to  movements  from  applied  loads.  In  other 
words,  well  rounded  particles  will  have  less  resistance  to  movement  due  to  either  static  or 
dynamic  loading  than  would  angular  particles.  From  the  two-dimensional  images  of 
sections  through  the  specimen,  most  of  the  particles  were  elliptical  in  shape. 

E.  COMPARISONS  OF  PARTICLE  ORIENTATION  ANALYSIS  RESULTS 

Rose  plots  generated  from  the  two-dimensional  particle  orientation  analyses  are  given 
in  Appendix  A.  A  detailed  description  of  the  statistics  used  in  the  Rosy  program  has  been 
given  in  SECTION  IV.  From  an  in-depth  review  of  all  statistical  parameters  generated  by 
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the  Rosy  program,  it  was  determined  that  the  vector  magnitude  was  the  only  parameter  that 
related  compressional  wave  speed,  transmission  ratio,  microstructure,  and  saturation 
together. 

Vector  magnitude  is  a  measure  of  the  dispersion  of  azimuths  in  a  data  set  and  values 
range  from  zero  to  one.  Vector  magnitude  values  near  zero  indicate  that  there  is  no 
preferred  particle  orientation.  The  data  obtained  in  this  study  implies  that  low  values  of 
vector  magnitude  may  be  related  to  higher  values  of  soil  stiffness,  compaction  energy, 
compressional  wave  speed,  and  transmission  ratio,  i.e.,  an  inverse  relationship  may  exist 
between  the  soil  parameters.  As  previously  explained  in  SECTION  V,  at  intermediate 
saturations  (20  percent  to  60  percent)  the  soil  is  stiffer  and  resists  rearrangement  of  soil 
particles.  Thus,  a  lower  degree  of  preferred  particle  orientation  would  be  expected  and 
vector  magnitude  values  are  expected  to  be  lower  at  intermediate  saturations  which  may  be 
attributable  to  the  resistance  of  particle  movements  at  intermediate  saturations. 

1.  Vector  Magnitude  Results 

Two  separate  categories  were  analyzed  from  the  two-dimensional  orientation  data: 
NO  SHPB  results,  and  SHPB  results,  respectively.  The  reason  for  performing  the 
different  analyses  is  to  observe  the  differences  in  soil  microstructure  in  compacted 
specimens  before  and  after  testing  in  the  split-Hopldnson  pressure  bar. 

a.  Horizontal  and  Vertical  Planes 

Figures  53  and  54  show  the  vector  magnitude  as  a  fimction  of  saturation  in 
the  horizontal  and  vertical  planes  for  NO  SHPB  and  SHPB  specimens,  respectively.  As 
expected,  SHPB  specimens  had  higher  values  of  vector  magnitude  at  all  saturations  for 
each  plane.  Since  more  energy  is  applied  to  the  soil,  the  particles  are  forced  to  develop  a 
higher  degree  of  preferred  orientation. 

The  Ottawa  20-30  sand  used  in  this  investigation  has  been  classified  as  a 
poorly-graded,  well-rounded  sand.  The  mean  axial  ratio  of  the  sand  is  0.647  and  most  of 
the  particles  are  somewhat  elliptical  in  shape,  having  smooth  surfaces.  For  this  shape, 
most  of  the  soil  particles  would  tend  to  have  a  vector  mean  in  the  range  of  80-  to  110- 
degrees  or  260-  to  280-degrees  from  initial  deposition  for  drv  soil  conditions  only  (note 
that  vector  mean  values  are  bidirectional).  The  principle  behind  this  idea  is  that  an  elliptical 
shaped  particle  falling  through  the  air  will  tend  to  rest  on  its  longest  axis  parallel  to  the 
plane  of  deposition.  This  was  confirmed  by  Mitchell  (1976)  in  which  dry  pluviated 
samples  showed  a  strong  degree  of  preferred  particle  orientation  parallel  to  the  plane  of 


118 


deposition.  As  energy  is  applied,  the  sand  is  forced  to  rearrange  itself  into  a  denser 
packing  and  develop  high  vector  magnitude  values. 

It  was  expected  that  the  vector  magnitude  as  a  function  of  saturation  would 
exhibit  an  inverse  relationship  as  compared  to  compaction  energy,  compressional  wave 
speed  and  stress  transmission  ratio  as  a  function  of  saturation.  In  the  cases  of  compaction 
energy,  compressional  wave  speed  and  stress  transmission  ratio  as  a  function  of  saturation, 
each  of  these  exhibited  a  typical  trend  in  that  at  intermediate  saturations  (20  percent  to  60 
percent)  higher  values  were  observed.  Data  for  vector  magnitude  as  a  function  of 
saturation  imply  that  lower  values  are  encountered  at  intermediate  saturations  (20  percent  to 
60  percent).  In  general,  this  trend  was  observed  in  both  the  horizontal  and  vertical  planes, 
as  shown  in  Figures  53  and  54.  SHPB  specimens  exhibited  greater  vector  magnitude 
values  for  both  the  horizontal  and  vertical  planes. 

b.  Layers  1  Through  4 

Figures  72  and  73  show  the  vector  magnitude  as  a  function  of  saturation  and 
layer  number  (developed  from  vertical  plane  data)  for  NO  SHPB  and  SHPB  specimens, 
respectively.  The  data  have  been  normalized  to  the  vector  magnitude  at  zero  saturation  for 

Vector  magmtude  values  for  each  layer  were  assumed  to  decrease  with  successive 
layers,  i.e..  Layer  1  should  have  the  highest  vector  magnitude  and  Layer  4  should  have  the 
lowest  vector  magmtude.  Since  Layer  1  has  had  much  more  energy  applied  to  it  from  the 
iOlal  number  of  blows,  the  particles  probably  have  a  higher  degree  of  preferred  orientation, 
i.e.,  a  higher  vector  magnitude. 

Figures  55  through  58  show  the  vector  magnitude  as  a  function  of  saturation 
for  each  layer  and  compare  vector  magnitude  values  for  NO  SHPB  and  SHPB  specimens. 
As  the  amount  of  energy  applied  to  a  soil  increases  the  degree  of  preferred  particle 
orientation  should  increase.  This  is  seen  in  each  of  the  four  layers  and  saturations  in  that 
the  SHPB  specimens  generally  had  higher  vector  magnitude  values. 

2 .  Qjrrelation  of  Vector  Magnitude  with  Average  Total  Blow  Count 

The  next  step  in  the  analysis  was  to  take  the  product  of  vector  magnitude  and 
compaction  energy  as  a  function  of  saturation  and  to  note  any  distinct  trends  that  may  be 
occurring  in  the  data.  These  products  as  a  function  of  saturation  may  constitute  a  first  order 
or  second  order  function  on  a  best  fit  curve  which  may  be  of  use  in  future  studies.  This 
comes  from  the  inverse  relationship  between  vector  magnitude  and  compaction  energy, 
compressional  wave  speed,  and  transmission  ratio  as  a  function  of  saturation. 
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Figure  72.  Vector  Magnitude  Normalized  to  0%  at  Layer  #1  as  a  Function  of 
Saturation  for  Ottawa  20-30  Sand  (NO  SHPB). 


Figure  73.  Vector  Magnitude  Normalized  to  0%  at  Layer  #1  as  a  Function  of 
Saturation  for  Ottawa  20-30  Sand  (SHPB). 


a.  Horizontal  and  Vertical  Planes 


Figures  74  and  75  show  the  product  of  vector  magnitude  and  average  total 
blow  count  as  a  function  of  saturation  for  the  horizontal  and  vertical  planes  for  NO  SHPB 
and  SHPB  specimens,  respectively.  Both  the  horizontal  and  vertical  plane  results  in  Figure 
74  show  that  as  saturation  increases  the  product  steadily  increases.  In  Figure  75  the 
product  for  both  the  horizontal  and  vertical  planes  increases  to  about  40  percent  higher 
saturation  and  decreases  thereafter.  The  values  for  blow  count  used  in  both  NO  SHPB  and 
SHPB  specimens  are  the  same.  However,  values  for  vector  magnitude  are  greater  for 
SHPB  specimens  than  for  NO  SHPB  specimens.  Therefore,  the  curves  developed  in 
Figures  74  and  75  are  dependent  on  the  vector  magnitude  values. 

Overall,  the  NO  SHPB  specimens  showed  an  increase  in  the  product  as 
saturations  increased,  whereas  the  product  in  the  SHPB  specimens  reached  a  peak  value 
and  decreased  thereafter.  This  seems  to  correlate  with  compaction  energy,  cotr5)ressional 
wave  speed,  and  stress  transmission  ratio  as  a  function  of  saturation  which  generally  had 
maximum  values  at  intermediate  saturations  (see  Figures  59,  62,  and  63).  Both  the  vertical 
and  horizontal  plane  exhibited  a  similar  trend  for  both  NO  SHPB  and  SHPB  specimens. 
Overall,  the  product  of  vector  magnimde  and  average  total  blow  count  presented  in  this 
analysis  shows  small  variations  in  the  results  between  the  horizontal  and  vertical  plane. 
This  suggests  that  soil  microstructure  throughout  the  specimen  may  be  similar. 

Figtues  76  and  77  show  normalized  vector  magnitude  and  normalized  average 
total  blow  count  as  a  function  of  saturation  (normalized  to  results  at  S  =  0  percent)  for  NO 
SHPB  and  SHPB  specimens,  respectively.  Normalizing  vector  magnitude  and  average 
total  blow  count  (plotted  using  the  same  scales  on  both  axes)  allows  for  direct  comparison 
of  trends.  Both  rigures  show  an  inverse  relationship  between  vector  magnitude  and 
average  total  blow  count.  At  intermediate  saturations,  vector  magnitude  values  are  lower 
while  average  total  blow  count  values  are  higher.  The  lower  vector  magnitude  values 
indicate  a  more  random  orientation  within  the  soil  structure.  This  random  orientation  can 
be  attributed  to  particles  resisting  movement  relative  to  one  another  during  the  compaction 
process  when  moisture  is  present.  The  results  shown  in  Figures  76  and  77  suggest  that  a 
random  orientation  may  create  a  stiffer  soil,  and  thus  more  energy  is  required  to  compact 
the  moist  soil  to  a  constant  dry  density. 

b.  Layers  1  Through  4 

Figures  78  and  79  show  the  product  of  vector  magnitude  and  average  blow 
count  as  a  function  of  saturation  for  Layers  1  through  4  for  NO  SHPB  and  SHPB 
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Figure  74.  Product  of  Vector  Magnitude  and  Average  Total  Blow  Count  as  a  Function 
of  Saturation  for  Ottawa  20-30  Sand  (NO  SHPB). 
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Figure  75.  Product  of  Vector  Magnitude  and  Average  Total  Blow  Count  as  a  Function 
of  Saturation  for  Ottawa  20-30  Sand  (SHPB). 
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Figure  77.  Normalized  Vector  Magnitude  in  the  Vertical  Plane  and 
Normalized  Compactive  Energy  as  a  Function  of  Saturation  (SHPB). 
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Figure  78.  Product  of  Vector  Magnitude  and  Average  Blow  Count  as  a  Function  of 
Saturation  for  Ottawa  20-30  Sand  (NO  SHPB). 
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Figure  79.  Product  of  Vector  Magnitude  and  Average  Blow  Count  as  a  Funtion 
of  Saturation  for  Ottawa  20-30  Sand  (SHPB). 


specimens,  respectively.  In  each  figure,  the  prcxiuct  of  vector  magnitude  and  total  blow 
count  indicates  that  Layer  1  has  the  greatest  value  at  all  saturations.  In  Layers  2,  3,  and  4, 
the  product  of  vector  magnitude  and  total  compaction  energy  decrease  with  respect  to 
saturation.  However,  in  most  cases  values  for  vector  magnitude  are  greater  for  SHPB 
specimens  than  for  NO  SHPB  specimens  and  therefore,  the  curves  developed  in  Figures  78 
and  79  are  dependent  on  the  vector  magnitude  values. 

Overall,  the  NO  SHPB  specimens  showed  a  more  uniform  increase  in  the 
product  of  vector  magnitude  and  average  blow  count  as  saturations  increased,  whereas  the 
product  for  the  SHPB  specimens  varied  in  that  there  was  more  dispersion  within  each 
layer.  This  relates  to  compaction  energy,  compressional  wave  speed,  and  stress 
transmission  ratio  as  a  function  of  saturation  which  generally  had  maximum  values  at 
intermediate  saturations. 

3.  Product  of  Vector  Magnitude  and  Normalized  Compressional  Wave  Speed  and 

Stress  Transmission  Ratio 

Figures  80  and  81  show  the  product  of  vector  magnitude  and  normalized  wave 
speed,  and  the  product  of  vector  magnitude  and  transmission  ratio  as  a  function  of 
saturation,  respectively,  for  SHPB  specimen  sections  in  the  horizontal  and  vertical  planes. 
Results  in  the  horizontal  and  vertical  planes  for  Figures  80  and  8 1  show  a  similar  trend  in 
that  the  product  generally  decreases  with  increasing  saturation.  Figures  82  and  83  show 
the  product  of  normalized  vector  magnitude  and  normalized  wave  speed  and  transmission 
ratio  normalized  to  unity  at  zero  percent  saturation  as  a  function  of  saturation,  respectively, 
for  SHPB  specimen  sections  in  the  horizontal  and  vertical  planes. 

The  products  of  normalized  products  in  Figures  82  and  83  for  the  vertical  and 
horizontal  planes  show  similar  results  to  those  shown  in  Figures  80  and  81.  Figures  82 
and  83  show  a  similar  trend  in  that  the  product  generally  decreases  with  increasing 
saturation.  The  behavior  in  the  vertical  and  horizontal  planes  of  Figure  83  is  a  bit  more 
pronounced  than  that  of  Figure  82.  Overall,  the  product  of  vector  magnitude  and 
normalized  wave  speed  and  transmission  ratio  presented  in  this  analysis  show  small 
variations  in  the  results  between  the  horizontal  and  vertical  planes.  This  suggests  that  soil 
micTostructure  throughout  the  specimen  may  be  similar. 

Figures  84  and  85  show  normalized  vector  magnitude  for  NO  SHPB  and  SHPB 
specimens  as  a  function  of  saturation  (normalized  to  the  results  at  S  =  0  percent)  in  the 
vertical  and  horizontal  planes,  respectively.  For  SHPB  specimens,  vertical  and  horizontal 
planes  exhibit  the  same  general  trend  in  that  the  normalized  vector  magnitude  decreases  as 
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Figure  80.  Product  of  Vector  Magnitude  and  Normalized  Wave  Speed  as  a  Function  of 
Saturation  Ottav^a  20-30  Sand  (SHPB). 


Figure  81.  Product  of  Vector  Magnitude  and  Normalized  Transmission  Ratio  as  a 
Function  of  Saturation  Ottawa  20-30  Sand  (SHPB). 
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Figure  82.  Product  of  Normalized  Vector  Magnitude  and  Normalized  Wave  Speed  as 
a  Function  of  Saturation  for  Ottawa  20*30  Sand  (SHPB). 


Figure  83.  Product  of  Normalized  Vector  Magnitude  and  Normalized  Transmission 
Ratio  as  Function  of  Saturation  for  Ottawa  20-30  Sand  (SHPB). 
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Figure  84.  Normalized  Vector  Magnitude  in  the  Vertical  Plane  as  a  Function 
of  Saturation  for  Ottawa  20-30  Sand  Compacted  Moist  and  Tested  Moist. 
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Figure  85.  Normalized  Vector  Magnitude  in  the  Horizontal  Plane  as  a  Function 
of  Saturation  for  Ottawa  20-30  Sand  Compacted  Moist  and  Tested  Moist. 


saturation  increases.  However,  for  NO  SHPB  specimens,  the  normalized  vector 
magnitude  in  the  vertical  and  horizontal  planes  initially  decrease  to  about  20  percent 
saturation  and  increase  thereafter. 

Figures  86  and  87  show  normalized  vector  magnitude  (NO  SHPB)  and 
normalized  compressional  wave  speed  (moist-moist)  as  a  function  of  saturation, 
respectively,  (normalized  to  results  at  S  =  0  percent)  for  the  vertical  and  horizontal  planes. 
Both  figures  show  an  inverse  relationship  between  vector  magnitude  and  compressional 
wave  speed.  At  intermediate  saturations,  vector  magnitude  values  are  lower  while 
compressional  wave  speed  values  are  higher.  Figures  88  and  89  show  normalized  vector 
magnitude  (NO  SHPB)  and  normalized  transmission  ratio  (moist-moist)  as  a  function  of 
saturation  (normalized  to  the  results  at  S  =  0  percent)  for  the  vertical  and  horizontal  planes, 
respectively.  Both  figures  show  an  inverse  relationship  between  vector  magnitude  and 
transmission  ratio.  At  intermediate  saturations,  vector  magnitude  values  are  lower  while 
transmission  ratio  values  are  higher. 

Figures  90  and  91  show  normalized  vector  magnitude  (NO  SHPB)  and 
normalized  compressional  wave  speed  (moist-dry)  as  a  function  of  saturation  (normalized 
at  zero  percent  saturation)  fw  the  vertical  and  horizontal  planes,  respectively.  Both  figures 
show  an  inverse  relationship  between  vector  magnitude  and  compressional  wave  speed.  At 
intermediate  saturations,  vector  magnitude  values  are  lower  while  compressional  wave 
speed  values  are  higher.  Figures  92  and  93  show  normalized  vector  magnitude  (NO 
SHPB)  and  normalized  transmission  ratio  (moist-dry)  as  a  function  of  saturation 
(normalized  at  zero  percent  saturation)  for  the  vertical  and  horizontal  planes,  respectively. 
Both  figures  show  an  inverse  relationship  between  vector  magnitude  and  transmission 
ratio.  At  intermediate  saturations,  vector  magnitude  values  are  lower  while  transmission 
ratio  values  are  higher. 

The  results  shown  in  Figures  86  through  93  suggest  that  the  lower  vector 
magnitude  values  indicate  a  more  random  particle  orientation  within  the  soil  structure.  This 
random  orientation  can  be  attributed  to  particles  resisting  movement  relative  to  one  another 
during  the  compaction  process  when  moisture  is  present.  This  implies  that  a  random 
orientation  creates  a  stiffer  soil,  thus  producing  higher  values  of  wave  speed  and 
transmission  ratio  at  intermediate  saturations  for  Ottawa  20-30  sand  compacted  to  a 
constant  dry  density. 

Figure  94  shows  a  comparison  of  compressional  wave  speed  values  obtained 
firom  in  the  split-Hopkinson  pressure  bar  tests  and  theoretical  calculations  using  constrained 
modulus  values  obtained  from  the  quasi-static  compression  test  results  (Figures  38  through 
43).  Theoretical  values  of  compressional  wave  speed  were  calculated  using  Equation  (13) 
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Figure  86.  Normalized  Vector  Magnitude  (NO  SHPB)  in  the  Vertical  Plane 
and  Normalized  Wave  Speed  (moist-moist)  as  a  Function  of  Saturation. 
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Figure  87.  Normalized  Vector  Magnitude  (NO  SHPB)  in  the  Horizontal  Plane 
and  Normalized  Wave  Speed  (moist-moist)  as  a  Function  of  Saturation. 
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Figure  88.  Normalized  Vector  Magnitude  (NO  SHPB)  in  the  Vertical  Plane 
and  Normalized  Transmission  Ratio  (moist-moist)  as  a  Function  of 
Saturation. 


CXtawa  20-30  Sand 


Oliva  MoissmsMvai  aaznvwaoN 


Figure  89.  Normalized  Vector  Magnitude  (NO  SHPB)  in  the  Horizontal  Plane 
and  Normalized  Transmission  Ratio  (moist-moist)  as  a  Function  of  Saturation. 
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Figure  90.  Normalized  Vector  Magnitude  (NO  SHPB)  in  the  Vertical  Plane 
and  Normalized  Wave  Speed  (moist-dry)  as  a  Function  of  Saturation. 
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Figure  91.  Normalized  Vector  Magnitude  (NO  SHPB)  in  the  Horizontal  Plane 
and  Normalized  Wave  Speed  (moist-dry)  as  a  Function  of  Saturation. 
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Figure  92.  Normalized  Vector  Magnitude  (NO  SHPB)  in  the  Vertical  Plane  and 
Normalized  Transmission  Ratio  (moist-dry)  as  a  Function  of  Saturation. 
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Figure  94.  Comparison  of  Average  Experimental  and  Calculated  Wave 
Speed  as  a  Function  of  Saturation  for  Moist/Moist  Tests. 
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Figure  93.  Normalized  Vector  Magnitude  (NO  SHPB)  in  the  Horizontal  Plane 
and  Normalized  Transmission  Ratio  (moist-dry)  as  a  Function  of  Saturation. 


and  replacing  the  modulus  of  elasticity,  E,  with  the  constrained  modulus,  D.  At  all 
saturations,  the  theoretical  wave  speeds  are  greater  than  experimental  values.  The 
theoretical  wave  speed  slightly  increases  from  zero  to  about  20  percent  saturation  and 
gradually  decreases  thereafter.  At  saturations  greater  than  about  40  percent  both 
experimental  and  theoretical  wave  speed  curves  exhibit  the  same  general  trend.  Table  9 
shows  results  of  constrained  modulus  and  theoretical  compressional  wave  speed. 

TABLE  9.  RESULTS  OF  CONSTRAINED  MODULUS  AND 
CALCULATED  WAVE  SPEED  FOR  OTTAWA  20-30  SAND 
COMPACTED  MOIST  AND  TESTED  MOIST 


s 

(%) 

Total 

Unit  Weight 
(g/cc) 

Constrained 

Modulus 

(kPa) 

Normalized 

Constrained 

Modulus 

Calculated 

Wave  Speed 
(m/sec) 

Normalized 

Wave  Speed 

0.0 

1.715 

13595 

1.000 

614.3 

1.000 

18.7 

1.781 

14463 

1.064 

621.7 

1.012 

37.4 

1.847 

14740 

1.084 

616.3 

1.003 

56.2 

1.913 

13814 

1.016 

586.3 

0.954 

74.9 

1.979 

14041 

1.033 

581.1 

0.946 

93.6 

2.045 

13131 

0.966 

552.8 

0.900 

4 .  Comparison  of  Vector  Magnitude  and  Rayleigh  Values 

Vector  magnitude  and  Rayleigh  test  of  significance  values  can  be  related  as 
follows.  A  preferred  particle  orientation  exists  in  a  soil  when  the  Rayleigh  value,  p,  is  less 
than  0.05  (see  Equation  (10)).  The  corresponding  vector  magnitude  value  was  generally  in 
the  range  of  0.1  to  0.5  for  this  investigation;  however  the  highest  obtainable  value  is  1.0. 
A  correlation  between  the  confidence  angle  and  Rayleigh  value  can  also  be  observed  such 
that  a  smaller  confidence  angle  is  present  when  the  Rayleigh  value  is  below  0.05.  This 
implies  that  the  dispersion  within  the  azimuth  readings  is  small  at  lower  Rayleigh  values. 

a.  Horizontal  and  Vertical  Planes 

Similar  vector  magnitude  and  Rayleigh  values  were  obtained  in  the  horizontal 
and  vertical  planes.  All  Rayleigh  values  were  less  than  0.05  for  NO  SHPB  and  SHPB 
specimens.  The  only  exception  was  the  Rayleigh  value  at  saturations  of  18.7  percent  for 
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NO  SHPB  specimens  in  the  vertical  plane.  One  possible  reason  for  this  is  that  the  number 
of  particles  measured  were  less  than  those  measured  at  other  saturations,  or  that  no 
preferred  particle  orientations  exist  for  that  section.  The  number  of  particles  measured,  N, 
does  affect  the  Rayleigh  value.  In  comparing  the  vertical  and  horizontal  planes,  the 
sections  that  showed  a  preferred  particle  orientation  had  a  number  of  particle  measurements 
greater  than  or  equal  to  320.  This  coincides  with  the  findings  from  other  researchers  such 
as  Oda  (1972b)  and  Campbell  (1985)  who  said  that  a  minimum  of  200  to  3(X)  particles  need 
to  be  measured  in  an  orientation  analysis. 

The  lower  Rayleigh  values  and  higher  vector  magnitude  values  obtained  in  the 
horizontal  and  vertical  planes  may  be  attributed  to  the  location  of  Section  A-A  (see  Figures 
22  and  23).  Ottawa  20-30  sand  is  well  rounded  and  elliptically  shaped  in  two  dimensions, 
thus  particles  will  tend  to  rest  on  the  longest  axis  parallel  to  the  plane  of  deposition  (90  or 
270-degrees)  for  drv  conditions  only  (see  previous  discussion  in  SECTION  V).  If  a 
majority  of  the  particles  follow  this  trend  then  it  is  expected  that  lower  Rayleigh  values  and 
higher  vector  magnitude  values  will  be  obtained. 

Stresses  developed  during  compaction  using  the  standard  Proctor  hammer 
force  the  sand  to  rearrange  itself  into  a  denser  packing  after  initial  placement.  During 
SHPB  testing  the  density  increases  further  and  the  soil  particles  are  forced  into  a  higher 
degree  of  preferred  orientation  by  rotating  a  majority  of  the  particles  from  their  original 
orientation  back  along  the  plane  of  deposition  plus  or  minus  a  few  degrees  for  dry  and 
unsaturated  conditions.  This  may  explain  why  lower  Rayleigh  values  and  higher  vector 
magnitude  values  are  obtained  in  the  horizontal  and  vertical  planes  for  SHPB  specimens 
when  compared  with  the  NO  SHPB  results. 

Overall,  the  vertical  plane  results  had  the  highest  vector  magnitude  values  for 
both  NO  SHPB  and  SHPB  specimens.  The  vertical  plane  on  Section  A-A  (Figtires  22  and 
23)  is  at  the  centerline  of  the  specimen  and  at  a  constant  distance  from  the  boundary,  i.e., 
the  wall  of  the  specimen  container.  Therefore,  any  boundary  condition  effects  should  be 
uniform  over  the  entire  vertical  plane  since  the  specimen  is  cylindrical  in  shape.  The 
horizontal  plane  on  Section  A-A  is  not  at  a  constant  distance  from  the  boundary  and 
boundary  condition  effects  would  not  be  expected  to  be  uniform  over  the  entire  horizontal 
plane.  This  may  be  the  reason  for  lower  vector  magnitude  values  in  the  horizontal  plane. 
However,  the  trends  for  both  the  vertical  and  horizontal  directions  are  similar  (see  Figures 
80  through  83). 
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b.  Layers  1  Through  4 

The  vector  magnitude  and  Rayleigh  values  obtained  in  Layers  1  through  4 
were  typically  higher  than  those  for  the  horizontal  and  vertical  planes.  The  data  can  be  seen 
in  Tables  6  and  7  for  NO  SHPB  and  SHPB  specimens,  respectively.  In  each  layer  the 
Rayleigh  values  were  less  than  0.05  for  NO  SHPB  and  SHPB  specimens.  The  only 
exceptions  were  the  Rayleigh  values  at  a  saturation  of  18.7  percent  for  NO  SHPB  results 
for  Layer  3,  and  at  saturation  of  18.7  percent  and  37.4  percent  for  NO  SHPB  results  for 
Layer  4.  The  number  of  particles  measured  for  each  individual  layer  was  not  very  large, 
typically  on  the  order  of  120  to  220  particles.  This  does  not  coincide  with  the  findings  of 
Oda  (1972b)  and  Campbell  (1985),  who  stated  that  200  to  300  particles  need  to  be 
measured  for  particle  orientation  analysis.  A  possible  reason  for  higher  vector  magnitude 
values  and  lower  Rayleigh  values  for  Layers  1  through  4  can  be  attributed  to  the  location  of 
Section  A- A  (see  Figures  22  and  23).  Orientation  data  for  Layers  1  through  4  were 
obtained  from  the  vertical  plane,  developed  from  Section  A-A  (see  Figure  27),  and  a 
preferred  particle  orientation  was  evident  in  both  NO  SHPB  and  SHPB  specimens  along 
the  vertical  plane  (see  Hgure  72, 73, 78,  and  79). 

In  each  specimen,  the  first  layer  was  subjected  to  the  greatest  amount  of 
compaction  energy  (see  Figure  61).  Therefore,  it  was  expected  that  Layer  1  would  have 
low  Rayleigh  values  and  higher  vector  magnitude  values  than  Layers  2,  3,  and  4  due  to  the 
increased  number  of  blows  (i.e.,  higher  compaction  energy).  Even  though  Layer  1 
received  the  greatest  amount  of  compaction  energy  it  was  also  the  layer  closest  to  the  steel 
pedestal.  The  effect  of  the  reflected  tensile  wave  may  have  been  greatest  in  this  layer  and 
may  have  rearranged  the  particles,  thereby  reducing  the  vector  magnitude  value.  However, 
a  preferred  particle  orientation  does  exist  within  this  layer  for  both  the  NO  SHPB  and 
SHPB  specimens. 

Layers  2  and  3  were  subjected  to  an  intermediate  amount  of  compaction  energy 
(see  Figure  61).  Therefore,  it  was  expected  that  these  layers  would  have  similar  Rayleigh 
values  and  lower  vector  magnimde  values  than  Layer  1.  As  the  total  amount  of  compaction 
energy  decreases  (compared  to  Layer  1),  there  would  be  less  rearrangement  of  the  particles 
leading  to  lower  vector  magnitude  values.  However,  even  though  Layers  2  and  3  received 
an  intermediate  amount  of  compaction  energy  they  are  further  away  from  the  steel  pedestal, 
thus  the  effect  of  the  reflected  tensile  wave  may  be  lower  in  these  two  layers.  Therefore, 
the  particles  may  have  less  rearrangement  and  the  vector  magnitude  values  would  be  lower. 
A  preferred  particle  orientation  does,  however,  exist  within  both  layers  for  both  the  NO 
SHPB  and  SHPB  specimens. 
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Layer  4  was  subjected  to  the  least  amount  of  compaction  energy  (as  compared  to 
Layers  1,  2,  and  3,  see  Figure  61).  It  was  expected  that  Layer  4  would  have  higher 
Rayleigh  values  and  the  lowest  vector  magnitude  values.  Layer  4  is  the  furthest  away  from 
the  steel  pedestal  and  the  effect  of  the  reflected  tensile  wave  may  be  the  least  in  this  layer 
and  due  to  this,  particles  may  have  less  rearrangement,  therefore  vector  magnitude  values 
would  be  lower.  However,  a  preferred  particle  orientation  does  exist  within  this  layer  for 
both  the  NO  SHPB  and  SHPB  specimens. 
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SECTION  VI 


SUMMARY 

This  study  has  investigated  the  effects  of  saturation  on  soil  microstructure, 
compaction  energy,  compressional  wave  speed,  and  stress  transmission  ratio  in  unsaturated 
soils.  Soils  are  generally  in  an  unsaturated  state  in  the  field,  therefore  a  better 
understanding  of  their  behavior  during  dynamic  loading  is  useful. 

Specimens  of  Ottawa  20-30  sand  were  compacted  to  a  constant  dry  density  of  1.715 
g/cm^  in  a  2.54  cm  (1.0  inches)  wall  specimen  container  in  four  layers  using  a  standard 
Proctor  hammer.  Saturations  varied  from  zero  to  about  94  percent.  Some  specimens  were 
tested  under  quasi-static  loading  conditions  using  a  Material  Testing  System  (MTS)  device 
to  evaluate  the  soil's  constrained  modulus.  Dynamic  wave  propagation  tests  were 
performed  using  a  split-Hopkinson  pressure  bar  (SHPB)  which  applies  an  intense, 
transient  compressive  loading  to  the  specimen.  Compressional  wave  speed  and  stress 
transmission  ratio  can  be  evaluated  through  the  use  of  the  SHPB. 

The  amount  of  compaction  energy  required  to  compact  the  Ottawa  20-30  sand  to  a 
constant  dry  density  of  1.715  g/cm^  varied  with  saturation.  This  trend  has  been  observed 
by  several  other  researchers  for  Ottawa  20-30  sand  and  other  sands  (Ross  et  al.,  1986; 
Pierce,  1989;  Veyera,  1989).  There  appears  to  be  an  increase  in  soil  stiffness  at 
intermediate  saturations  which  is  most  likely  due  to  capillary  tension  forces  present  during 
compaction.  Capillary  tension  forces  may  resist  particle  rearrangement,  thus  requiring 
more  compaction  energy.  As  saturations  increase  beyond  about  60  percent,  the  additional 
moisture  may  "lubricate"  the  soil  particles  as  surface  tension  forces  are  reduced  which 
would  allow  particles  to  move  relative  to  one  another  with  less  resistance  decreasing  the 
amount  of  compactive  energy  required  to  obtain  a  constant  dry  density. 

More  compaction  energy  was  required  to  compact  the  first  layer  at  all  saturations  as 
compared  with  Layers  2,  3,  and  4.  In  all  four  layers,  less  energy  was  required  to  compact 
the  soil  to  a  constant  dry  density  at  zero  percent  saturation  (dry  soil).  More  energy  was 
required  to  compact  each  layer  to  a  constant  dry  density  at  intermediate  saturations  (20 
percent  to  60  percent).  The  effect  of  moisture  (capillary  tension  forces)  during  compaction 
for  each  individual  layer  is  also  evident  (see  Figure  61). 

Dynamic  split-Hopkinson  pressure  bar  (SHPB)  tests  were  performed  on  specimens 
compacted  moist  and  tested  moist  (moist-moist),  and  also  on  specimens  compacted  moist 
and  tested  dry  (moist-dry)  to  evaluate  the  influence  of  moisture  during  SHPB  testing. 
Results  showed  that  specimens  compacted  moist  and  tested  moist  had  lower  values  of 
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compressional  wave  speed  and  stress  transmission  ratio  at  all  saturations.  Thus,  the 
presence  of  moisture  during  SHPB  testing  does  influence  the  compressional  wave  speed 
and  stress  transmission  ratio. 

The  results  presented  herein  are  similar  to  those  obtained  by  previous  researchers 
who  performed  SHPB  tests  on  the  same  and  similar  soils  compacted  moist  and  tested  moist 
(Veyera,  1989;  Ross,  1989;  Felice  et.  al.,  1987;  and  Ross  et.  al.,  1986).  Results  of 
compaction  energy,  compressional  wave  speed,  stress  transmission  ratio,  and  constrained 
modulus  as  a  function  of  saturation  all  show  a  similar  trend.  At  intermediate  saturations 
(20  percent  to  75  i>ercent)  these  values  were  greater  than  those  at  zero  percent  saturation, 
and  may  be  due  to  an  increase  in  soil  stiffness  from  capillary  tension  forces.  However, 
Pierce  (1989)  noted  that  capillary  stresses  are  negligible  compared  to  stress  levels  during  in 
SHPB  testing  and  therefore  probably  do  not  influence  the  dynamic  response. 

Results  of  compressional  wave  speed  and  stress  transmission  ratio  as  a  function  of 
samration  in  this  investigation  differed  from  results  presented  by  Pierce  (1989)  who 
compacted  specimens  dd  and  saturated  and  desaturated  them  to  the  desired  level  by  means 
of  a  pressure  plate  apparams,  thus  the  influence  of  moisture  during  compaction  was  not 
present.  In  this  investigation  specimens  were  compacted  moist,  thus  the  influence  of 
moisture  was  present  during  compaction.  Differences  between  results  obtained  by  Pierce 
(1989),  Ross  (1986),  Veyera  (1989),  and  this  investigation  can  be  attributed  to  differences 
in  sample  preparation  techniques  particularly  the  presence  of  moisture  during  compaction. 

S[}ecimens  were  compacted  using  the  standard  Proctor  hammer  only  (NO  SHPB)  and 
also  compacted  using  the  standard  Proctor  hammer  and  tested  using  high  strain  rates  in  the 
split-Hopkinson  pressure  bar  (SHPB).  The  microstrucmre  from  both  NO  SHPB  and 
SHPB  specimens  was  preserved  by  the  use  of  a  two-part  epoxy  resin.  Sections  were  taken 
through  the  cylindrical  specimen  to  view  the  microstructure  along  different  planes. 
Photomicrographs  of  each  specimen  were  subdivided  into  three  different  planes,  i.e.,  the 
horizontal,  vertical,  and  semicircular  directions.  The  vertical  plane  was  further  subdivided 
into  four  2.54  cm  (1.0  in)  layers  such  that  the  microstructure  of  each  compacted  layer  could 
be  analyzed.  A  complete  microstructural  analysis  was  performed  including  an  axial  ratio 
and  particle  long-axis  orientation  study.  A  total  of  70  rose  plots  were  generated  by 
digitizing  the  two-dimensional  photomicrographs.  A  relationship  between  the  vector 
magnitude  (a  result  of  the  statistical  analysis  from  the  Rosy  program)  and  saturation  was 
analyzed.  The  product  of  vector  magnimde  and  compaction  energy,  vector  magnitude  and 
compressional  wave  speed,  and  vector  magnitude  and  transmission  ratio  was  also  analyzed 
as  a  function  of  saturation. 
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The  microstructural  analysis  produced  several  important  findings.  The  horizontal 
section,  vertical  section,  and  Layers  1  through  4  showed  that  a  preferred  orientation  exists 
for  both  NO  SHPB  and  SHPB  specimens.  Most  of  the  specimens  had  vector  mean  values 
in  the  range  of  260-  to  280-degrees.  The  semicircular  section  showed  no  preferred 
orientation  for  both  NO  SHPB  and  SHPB  specimens  which  is  most  likely  due  to  the  small 
number  of  particles  measured.  An  inverse  relationship  exists  between  compaction  energy, 
compressional  wave  speed,  and  stress  transmission  ratio,  and  vector  magnitude  as  a 
function  of  saturation  (i.e.,  at  intermediate  saturations  lower  values  of  vector  magnitude 
were  observed). 

Soil  microstructure  developed  during  compaction  is  influenced  by  the  moisture 
present  during  compaction.  This  can  be  seen  in  vector  magnitude  values  obtained  from  the 
particle  orientation  analysis  in  that  lower  values  are  obtained  at  intermediate  saturations  for 
all  parameters  investigated.  The  lower  vector  magnitude  values  can  be  related  to  soil 
stifftiess  and  the  data  implies  that  if  a  more  random  orientation  of  particles  exists  within  the 
grain  structure,  a  stiffer  soil  may  be  present.  Therefore,  it  appears  that  the  presence  of 
moisture  during  compaction  significantly  affects  the  resultant  soil  microstmcture  which  in 
turn  influences  compaction  energy,  compressional  wave  speed,  and  stress  transmission 
ratio. 

Results  from  quasi-static  compression  tests  and  constrained  modulus  data  were 
developed  from  stress-strain  curves  generated  by  the  MTS  device  using  four  compacted 
layers  and  one  compacted  layer  at  each  saturation  level.  Stress-strain  curves  at  zero  percent 
saturation  were  compared  to  stress-strain  curves  at  all  other  saturations.  There  were  some 
noticeable  differences  between  the  slopes  of  the  stress-strain  curves  with  variations  in 
saturation.  Stress-strain  curves  which  have  the  steepest  slope  indicate  a  stiffer  soil 
structure.  In  general,  as  saturations  are  increased  to  about  60  percent,  the  soil  became 
stiffer  as  compared  with  zero  percent  saturation.  As  saturations  increase  beyond  60 
percent,  the  soil  became  less  stiff.  This  compares  well  with  SHPB  test  results  for 
compressional  wave  speed  and  stress  transmission  ratio,  in  that  a  stiffer  material  is  present 
to  saturations  of  about  60  percent. 

Values  of  constrained  modulus  at  intermediate  saturations  (20  percent  to  60  percent) 
were  higher  than  at  all  other  saturations.  This  trend  was  also  observed  in  data  for 
compaction  energy  (Figure  59),  compressional  wave  speed  (Figure  62,  and  stress 
transmission  ratio  (Figure  63)  as  a  function  of  saturation.  The  behavior  may  be  due  to  the 
resistance  of  particle  movements  relative  to  one  another  at  intermediate  saturations.  By 
applying  a  quasi- static  loading,  the  soil  is  compressed  into  a  denser  packing  than  that 
produced  in  the  compaction  procedure.  At  lower  saturations,  capillary  forces  resist  particle 


152 


movements  relative  to  one  another.  As  saturation  increases,  capillary  forces  decrease, 
thereby  reducing  the  frictional  resistance  between  soil  particles  which  may  allow  particles  to 
move  relative  to  one  another  with  greater  ease. 
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SECTION  vn 


CONCLUSIONS 

Based  on  the  results  of  this  study  and  the  discussions  presented  in  previous  chapters, 
the  following  conclusions  were  drawn.  These  conclusions  pertain  to  the  laboratory 
investigation  of  Ottawa  20-30  sand  evaluated  under  the  experimental  conditions  previously 
stated  in  this  thesis. 

•  The  results  from  dynamic  SHPB  tests,  quasi-static  compression  tests  and 
microstructural  analyses  indicate  that  the  amount  of  moisture  present  during 
compaction  is  probably  the  most  significant  factor  influencing  both  the  dynamic 
and  static  behavior  of  Ottawa  20-30  sand. 

•  The  compactive  energy  required  to  produce  a  constant  dry  density  specimen  by 
dynamic  impact  compaction  is  strongly  dependent  on  the  moisture  present  during 
compaction.  The  greatest  amount  of  effort  was  necessary  when  compacting 
specimens  at  intermediate  saturations  (from  about  20  percent  to  60  percent).  This 
is  most  probably  due  to  the  influence  of  capillary  tension  forces  which  affect  the 
ease  with  which  grains  can  reorient  during  compaction. 

•  The  microstructural  analysis  results  suggest  that  an  inverse  relationship  exists 
between  vector  magnitude  and  compaction  energy,  compressional  wave  speed, 
stress  transmission  ratio  and  constrained  modulus  with  variations  in  the  degree  of 
saturation.  Lower  vector  magnitude  values  indicate  a  more  random  particle 
orientation  and  higher  specimen  stiffness  at  intermediate  saturations.  Higher 
vector  magnitude  values  indicate  a  more  preferred  particle  orientation  and  lower 
stiffness  for  dry  soil  (zero  saturation)  and  at  saturations  above  about  75  percent 

•  For  SHPB  specimens  compacted  moist  and  tested  moist,  the  transmission  ratio 
increases  from  zero  percent  to  about  40  percent  saturation,  and  then  decreases  with 
increasing  saturation.  The  wave  speed  increases  from  zero  percent  to  about  20 
percent  saturation,  and  remains  relatively  constant  thereafter.  This  implies  that  the 
presence  of  moistiu'e  during  testing  may  have  a  lubricating  effect,  allowing  the 
particles  to  move  somewhat  relative  to  each  other  during  loading.  However,  the 
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influence  of  noicrostructure  formed  during  compaction  is  a  significant  contributing 
factor. 

•  For  SHPB  specimens  compacted  moist  and  tested  dry,  increases  in  transmission 
ratio  and  wave  speed  are  generally  larger  than  for  specimens  compacted  moist  and 
tested  moist,  particularly  at  about  20  percent  saturation,  after  which  they  decrease 
with  increasing  saturation.  No  plateau  region  exists  at  intermediate  saturations  and 
it  appears  that  the  microstructural  characteristics  developed  during  compaction 
which  affect  the  dynanuc  properties,  remain  intact  even  after  the  moisture  in  the 
pores  has  been  removed. 

•  Quasi-static  uniaxial  compression  test  results  show  changes  in  specimen  stiffness 
with  changes  in  compaction  moisture  which  are  generally  consistent  with  trends 
observed  in  the  SHPB  dynamic  test  data. 
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SECTION  vni 


RECOMMENDATIONS 

After  a  thorough  review  of  the  results  obtained  in  this  investigation  and  those  from 
previous  work,  the  following  recommendations  for  further  research  are  made: 

•  The  influence  of  specimen  preparation  techniques  on  stress  transmission  ratio, 
wave  speed  and  vector  magnitude  results  should  be  examined.  Other  possible 
compaction  methods  include  pluviation,  vibration  and  static  densification. 

•  A  model  to  predict  compressional  wave  speed  and  stress  transmission  in 
unsaturated  soils  using  laboratory  results  for  wave  speed,  stress  transmission 
ratio,  and  vector  magnitude  should  be  developed.  The  model  should  be  verified 
and  improved  through  a  series  of  carefully  controlled  explosive  field  tests. 

•  Scaled  tests  should  be  performed  in  the  centrifuge  using  actual  explosive  charges 
in  compacted,  unsaturated  soils.  This  would  have  the  benefit  of  being  able  to 
simulate  a  variety  of  field  conditions  in  a  laboratory  environment  and  would  also 
to  provide  an  independent  verification  of  the  SHPB  results  reported  herein. 

•  Microstructural  analyses  should  be  conducted  to  study  the  effects  of  particle  shape 
and  size  compacted  to  varying  degrees  of  dry  density.  These  analyses  may  show 
that  soil  type  and  dry  density  affect  wave  speed,  stress  transmission  ratio,  and 
vector  magnitude  results. 

•  Analyses  of  two-dimensional  particle  preferred  orientations  developed  under 
different  confinement  conditions  (i.e.,  different  specimen  container  wall 
thicknesses)  should  be  performed. 

•  An  analysis  of  the  preferred  orientation  of  pores  within  the  soil  should  be 
considered.  The  orientation  of  pores  which  contains  both  water  and  air  may 
influence  stress  transmission  in  an  unsaturated  soil. 
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Three-dimensional  characterization  of  soil  microstructure  should  be  undertaken. 
This  will  require  developing  new  techniques  and  could  use  such  devices  as  a  CAT 
scanner  or  MRI.  Results  would  be  very  useful  in  model  development. 
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APPENDIX  A:  Rose  Diagrams 
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Figure  A,1 .  Rtsuiti  from  staiisocai  analysis  using  Rosy  for 
W«a^4  cm  (1 .00  incn)  and  S«0%  in  tn# 
honzomai  Diana  (NO  SHPe). 
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figure  A.2.  Raauits  from  stansocal  analysis  using  Rosy  for 
W-i2.54  cm  (1 .00  incn)  and  S-18.7%  in  tna 
honzomai  qtana  (NO  SHPB). 
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Figure  A.4. 


Resuits  from  statisocal  analysis  using  Rosy  for 
W«2.54  cm  (1 .00  inert)  and  S«56.2%  m  me 
honzontai  piane  (NO  SHPB). 
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Ftguf«  A^.  RMUlts  ffwn  smtiMtieM  analysts  using  Rosy  for 
W«i54  an  ( 1 .00  inch)  and  S»74.9%  in  tha 
honzontal  plana  (NO  SHP8). 
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Ftgure  A. 6. 


Raaults  from  statisocai  analysis  using  Rosy  for 
W-2.54  cm  (1 .00  inch)  ano  S-0%  in  ma 
vertical  plana  (NO  SHP6). 
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Rgurt  A.7.  RMulta  from  statistical  analysis  using  Rosy  for 
W-2.54  cm  (1.00  inch)  and  S«  18.7%  in  ma 
^«rocal  plana  (NO  SHP8). 
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Ftgure  A.8.  Rasuits  from  statistical  analysts  using  Rosy  for 
W-2.54  cm  (1 .00  incn)  and  S-37.4%  .n  tna 
vamcai  plana  (NO  SHPB). 
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F?guf#  A.9,  RmUIb  from  otataticii  anoiysm  using  Rosy  for 
W-2.54  cm  (1 .00  incn)  and  S-56.2%  in  m# 
vMcal  ptana  (NO  SHPB). 
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VERTXX*I,  (NOSHPS) 


N.  704 

CtM  intaiva  -  10  riigmi 
Mawnum  Ptreartaga-  fliS 


VMT  M«an  -  30.1 
CafS.Angia-  12.44 

R  MaoMud*  -  0.228 


Maan  PweMga  -  3.36  StandMd  Ommon  •  2.04 


Ftgure  A.  10.  Results  from  staosocai  analysis  using  Rosy  for 
W-2.54  cm  ( 1 .00  incn)  and  S-74.9%  m  tne 
verticai  piane  (NO  SHPB). 


W«2.54om  S^3%  SCMI-ORCUUAR  (NOSHP9) 

StaMn 

N.244 

dsMirtsTva*  lOdsgfM 

MaamumPflmmga* 

Mmh  Pwosntaos  •  &aa  Standwd  Oiwlnii  •  2.0* 

V«aar  Mmr  •  79.4 

Coftf.  >yioi»  •  227.13 
n  Majniiirta  «  0.029 
RayMgh  -  0.8SB7 

Figur«  A.  1 1 .  RMUts  from  sooisticai  anirtyta  using  Rosy  for 

W«2^cmO-00  inch)  and  S-0%in  tho 
somhoroiar  ptena  (NO  SHP8). 


V¥^54cni  S»t&7%  SEMhCIRCXaAR  (NO  SHPB) 

SIM. 

N.290 

ClSM  InMrva  m  10  dagiaas 

UsacmimPsreanoioao  9l2 

M««i  Pwoamags  -  5.96  SMndwd  Orn/mrrn  »  2.16 

v«onrM«ai.272J 

Cort.  Angia  •  63J29 

R  MagnWirts*  0.078 
r\m0m^  m  0.2146 

Ftgura  A.: 2.  Rasufts  from  statisbcal  analyse  using  Rosy  for 

W-2.54  cm  ( 1 .00  men)  and  S«  \  8.7%  in  tfia 
samt-orcuiar  plana  (NO  SHPB). 


W-^54om  S-37.4%  SEM^CWCUAfl  (NOS^«^ 

StMMdOi 

N-310 

OsM  irMivil  -  to  degr<M 

Msxmum  Pwwage  «  S.7 

M«an  PwoMMo*  -  •^Tsrrf  Tfirtinn  l  M 

VMrUM-2S.3 

Cof4.  Angie  •  37.98 

R  Meqneuoe  -  0.032 
Reyliign  -  0v42S6 

Ftguf*  A.1 3.  RMUtB  from  3tiiatical  arufysM  us«ig  Rosy  for 

cm  (1 .00  inch)  and  S-37.4%  in  tfia 
samMareUar  ptena  (NO  SHP8). 


SMtfce  ! 

N-28a 

Oeee  Interva  -  10  digw 

Menmum  Peroemage  -  flLl 

Mm  Pereantage  -  3.38  Stondve  OavMtan  -  2A4 

V«cHrM«an-291J 

Con.  Angia  -  33.89 

R  Megnaurti  •  OMI 
Rayla^- 0.1088 

Figure  A.  U.  Resuita  from  statiaticaf  anafyats  using  Rosy  for 
W-2.54  cm  ( 1 .00  incn)  and  S«56.2%  in  tna 
semi-orcuJaf  ofana  (NO  SHPB). 
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W«2.54  on  S.74.9%  SEMM:inCULAR  (NO  SHPS) 

N«237 

am»  imervii •  i0  aiyew 

Maaomum  Poewge  ■  9A 

CaiS.AnQie-dS.43 

R  Megniiirie  -  0.001 
RayM^.  0^4183  | 

Rgur«  A.  1 5.  RMUts  from  stuslieai  anafysia  us«>g  Rosy  for 
'HmZM  cm  <1 .00  incn)  ar)d  Si74.9%  in  m« 
somKarcUar  ptena  (NO  SHP6). 


W«i54cm  S-0%  LAYER  #1  (NOS«»01 

StMMOee  1 

N-ee 

Claee  Irtenrn  •  10  Oegroee 

Moianiiiin  Pwawtage- 21.7 

Mm  Pwoon^  -  5.50  StancM  0«n«lon  -  4.tt 

V«eiarMM-290j0 

CofS.  Angle -2S.S2 
RMi«n«udn«0J63 
RayMi^- 0.0001 

Figure  A.1 6.  Resuits  from  staos&cat  anaiysis  usrg  Rosy  for 
W«2.54  cm  (1 .00  incn)  ana  S^%  in  layer  #1 
(NOSHPB). 


W«2.54cm  S^ia.7%  UKYEflgi  (N0SHP9) 
N«  147 


Ct4M  InwvHa  10 
M«anumP« 


StMIlOM 

V«eiarMMn«2792 


13.6 

56  SdMMOMttan*3.96 


Cort.  18.22 

fl  Mignaudi  ■  0.346 
Raytti^- 0.0000 


Rgur*  A.17.  Rtcults  from  statistical  analysis  usmg  Rosy  for 
W»Z54cm(1.00  Inch)  aodS«l8,7%in  laysrfi 
(NOSHPB). 


w-2.54  cm  S-37.4%  LAYER  #1  (NOSHPB) 

StaMOBi 

N«  160 

Omm  tntmrvwi  -  10  dwgtmm 

ManrnimPvaMeoa-  11.2 

Main  Pwoamga  •  5.56  Standvd  DaraHon  -  2.77 

Vsear  Maan  •  78J 

Cow.  An9la-21  J4 

R  MagrMJda  *  0.286 
Rayiagh- 0.0000 

Figure  A.  18.  Results  from  statisocai  analysis  using  Rosy  for 

W-2.54  cm  (1 .00  inch)  and  S«37.4%  in  layer  #1 
(NOSHPQ). 


W«2.Man  SmS6S£%  L>YEB  <1  (NOSHPq) 

N«  164 

Otm  ittmMi  ■  i0  digWM 

HmamumPmnmnaom*  \ZB 

Mmi  Pticfwga  •  4.36  SiMMd  CMfllon  ■  3.61 


_ Sntt6dc« 

Vwnr  M«an  •  32.8 
Co(t.>Wi9l»«  22.68 
RM^n«id»»0287 
-  0  JOOO 


Rgur*  A.1 9.  R664jl!a  from  stitisocai  anaiyvt  using  Rosy  for 
Wi.2.54  cm  (1 .00  incft)  and  S-56.2%  in  laysr  #1 
(NOSHP8). 


W«2.34an  S-;4.9%  LAYER  41  (NO  SHPB) 

s.^  1 

N-  177 

CUM  inl«rv«l  «  10  dnqmm 

ManmumPorcafttaoe*  i0.7 

Hmn  Pmtnmom  «  3.36  SUMrt  OmSbii  •  2A^ 

VMr  Mean  -  d6v4 
Co(S.>fiol»*31.lS 

R  Meqneirn  »  0.190 
Raying*  0.001 7 

Figure  A.20.  Results  from  stausticaf  analysis  using  Rosy  for 
W-2.54  cm  (1 .00  incn)  and  S»74.9%  in  layer  si 
(NO  SHPB). 
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W»2.5icm  Si4)%  UAYEfliC  (NOS^^q^ _ 

N-  197 

Om*  iMval  •  10  rliBWM 
M«wnufnP«o»m«g®«  ll^ 

MMnf>«QMaQB«S.d6  SMMOMtan.ajt 


VaoiarMMn«i270J 
Cart.  Angti  «  29.M 
RMooMud*.  0.190 
noirtagti- 0.0008 


Rguf*  A.21 .  RMulti  frofn  statiatical  analysis  usvio  Rosy  for 
W-2.54  cm  (1 .00  incn)  and  S^%  in  layar  #2 
(NOSHPB). 


W.2.54  cm  S- 18.7%  LAYER  #2  (NO  SHP9) 

StMDCS 

N«  151 

Ctaulntofvrt-  lObegieaa 

Menmum  Rsrasrt^e  •  9,8 

Mm  Pmmaga  •  5.58  SlaniM  OMMn  -  2.17 

VMrtMM-  73.3 

Cort.  Angle  •39.64 

R  Meqnftioe  ^  0.164 
Reylert"- 0.0173 

Figure  A.22.  Results  /rom  stansocal  analysis  using  Rosy  for 
W«2.54  cm  (1 .00  incn)  and  S«i8.7%  in  layer  #2 
(NO  SHP8). 
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N*  179 


CtaMMaivit*  10  rtigw 
M«muinP<rcTiag#«i  }ZZ 


V«eHr  Mm  •  27&2 
Coft.  AngM  *  27  JO 
R  Maon<M-OJl2 


Mmh  PweOTtaga  «  9l98 


RMign*  0.0003 


Rgur*  A23.  RMults  from  stuisiicai  analysts  using  Rosy  for 
W-Z54  cm  <  1 .00  inch)  and  S-37.4%  in  iayar  #2 
(NOSHPB). 


w-2.54  cm  S-0&2%%  LAYER  #2  (NOSHPB) 

StsMoa  I 

N-  181 

C1«M  irMrvM  -  10  degreas 

Mtaanum  Pareamege  ■  13J 

Maan  Parcamag*  «  S.se  Standard  Pax  ■raw  «  3.51 

VaoorMaan-aaj 

Card.  Angta-  19.97 

R  Uagniuda  -  0  J88 
Rapaigh- 0.0000 

Figure  A.24.  Results  from  statisDcal  analysis  using  Rosy  for 
W-2.54  cm  (1 .00  incn)  and  S-56.2%  m  layer  #2 
(NO  SHPB). 


Hgur*  K25,  RMtits  from  sdtistal  anilyM  using  Rosy  for 
'HmZJU  cm  n  -00  incn)  and  S-74.9%  in  iayar  #2 
(NOSHPBV 


Figure  A^6,  Results  from  staiisocai  analysis  using  Rosy  for 
W-a.54  cm  { 1 .00  incn)  and  S-0%  .n  layer  #3 
(NOSHPQ). 


W.2.54an  3-18.7%  U^YER  «3  {NOSHRS) 

St«Miea 

N-90 

CtSM  imsrvel  -  10  d»0feM 

Maonum  Pmeanagm  •  8^ 

Maan  Paroamaga  -  S.S8  SMndwd  Oavattan  «  2.18 

Vacaor  Mam  -  290^ 

Oort.  >U«oii  •  60.07 
RMaQfHuai-0.137 
Rayiaigrt- 0.1833 

Rgur#  A^7.  RMuita  from  saunticif  analyse  uain^  Rosy  for 
W-2^  cm  (1 .00  inert)  and  S-l  8.7%  in  layer  #3 
(NOSHPe). 


W-2.54  cm  S-37.4%  LAYER  63  (NO  SHP9) 

StaMies 

N-  I8e 

ClaM  intarval  -  iQ  daqfaaa 

Maumum Paroafiaoa  -  n.8 

Maan  Pareantaga  -  3.28  Slandart  OainaSofi  •  2.S1 

Vae«arUaan-29S>i 

CorT.  AngM.31.98 
RMagntuda  .0.181 

Rayiai^  .  0.0022 

Figure  A.28.  Results  from  staosocai  analysts  using  Rosy  for 
W-2.54  cm  {1.00  Inert)  and  S«07.4%  m  layer  #3 
{NOSHP9). 
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W.i2.54cm  S-6S.2%  LAYER  S3  (^OSHPB) 

smmm  I 

N.202 

ClflM  ifsarvii  •  lO  dsqiSM 

MsnnumPwosftaoe*  9L9 

Vaonr  Mmh  •  40.7 

Cort.  Aneie«40.i7 

R  Mssnsude  •  0.137 
nsymga  «  0.0220 

Rgur*  A^.  RMutts  from  iWateM  analysis  usmg  Rosy  for 
'HmZS^  an  (1 .00  incft)  and  S«56^  in  taysr  #3 
(NOSHPB). 


W.2.94cm  3*74.0%  LAYERS3  (NOSHPB) 

SMtte.  1 

N.218 

Cl4M  intervtf  «  10  dcgraas 

Maxmjm  Parcsmage  ■  12A 

Mmh  Pernmage  •  5.SS  asndwd  Osvwton  -  Z6 

VscMT  Mean  .273.0 
Cort.Aneia.  19.04 
RMaQntusa-0.27Q 
Rayts^n.  0.0000 

Figure  A.30.  Resuits  from  statisacal  analysis  using  Rosy  for 
W-2.54  cm  (1 .00  incn)  and  S-74.9%  m  layer  #3 
(NOSHPB). 
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W«2.54om  Si^  LAVER  #4  {NOSW>8) 

Sumhoi 

N-  I2t 

CtsM  ifSSfvil «  to  dagpMB 

MaxmjmPweanaoe*  ^ZJ5 

M«fln  Pweamags  •  5.56  SttndMd  Omohi  •  3JS 

V«»Mm-275J 

CorS.  Angis*  16.28 
RMegniude- 0.414 

0.0000 

rtQum  A^1 .  RMtits  from  staosticil  anoiysta  ustf)g  Rosy  for 
W-2^  cm  ( 1 .00  incft)  and  S«a%  in  iayw  #4 
(NOSHPQ). 


W-2.54  an  S-ia.7%  LAYER  M  (NO  SHPB) 

4—« 

N-  15 

Oats  imrvtf  -  10  digriM 

Maumum  Pqreanttge  •  20.0 

Maan  Pafoancega  •  9.09  SUndafO  Oamalon  -  4  JO 

Vaoor  Mean  -  338 J 

Cow.  AfiQla  «  84/43 
RMagnauda-0J38 
Raytaigli  •  0.42B1 

Figure  A.32.  Results  from  staosQcal  analysis  using  Rosy  for 
W-2.54  cm  {1 ,00  incn)  and  S-i  8.7%  in  layer  S4 
(NO  SHPB). 


!  W«2.M  an  S>37.4%  LAYER  **  (HQ  SHP9) 


Nm  101 

CteM)rMfv««  10 
M«KmifnP«eaf«ao** 


>Xlt 


V«enrMMn«290J 
CorOnoM  «  53.17 
n  M«0Mua»«0.146 

•  0.1178 


Rgur*  A.33.  Rttuits  froffi  siattttical  anafysit  using  Rosy  for 
y^mZ54  cm  ( 1 .00  incn)  and  3-37.4%  in  tayw*  #4 
(NOSHP0). 


W-iMom  S.fle.2%  LAYEHM  (NOSHPB) 

Smmka 

N«  181 

ClM  iiMfvrt  «  10  dagr—i 

MexvnumPwoantage*  i2.7 

Mm  Pwomge  .  5.5S  aMvd  Orntmion  .  3^i 

V«aarMM.271>i 
Cart.>ngt».2au}7 
n  MeoMub*  -  0.206 

0.0003 

Figure  A.34.  Results  from  statisocaf  analysis  using  Rosy  for 
W-2.54  cm  (1 .00  incn)  and  S-56.2%  in  layer  S4 
(NOSHP8). 
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V¥^5»  cm  S«74.9%  LAYEff  M  (ISQ  SHPB) 
N-  179 


Clam  imwvn  •  i0  aipwM 

Ce(«.  Angle  •27.8S 

ManmumPtroMUge*  11.2 

R  M«entud»«021l 

M— n  nwDwwgi  ■  5.39  StMert  OewMen  «  2.72 

nayiggh  -  0.0003 

Rgurt  AJ35.  RMuits  from  stuisbcal  analysis  using  Rosy  for 
W.2^  cm  (1 .00  incn)  and  S-74.9%  in  layar  #4 
(NOSHPB). 


W^SA  cm  S-0%  HOAtZOKTAI.  (SHP9) 

SWMiea  ! 

N-420 

Clam  iramrn  •  10  degme 

MenmumPwoMnge*  iA.3 

t4mn  Pmeatn^  m  S.S0  Standvd  0«wMDn  •  3.97 

VocHr  Mm  •  279>i 

Cove.  Angle  •d.OS 

R  Uegneirte  •  0.407 
Re^iwqn  -  0.0000 

Fgure  A.36.  Results  from  statisiical  analysis  using  Rosy  for 
W«2.S4  cm  (1 .00  incn)  and  S«0%  m  tne 
honzomai  plane  (SHPS). 


W«i54cni  S-»a.7X  HORIZONTAL  (SHf»0> 

SMoa  I 

N«  742 

Vaoor  Maan  -  2740 

CtaM  liaarva  •  10  tSwqmm 

CortAnoa*  iZOS 

MmrvjmPweamaoa*  ^2A 

R  Mapnajot  ■  0.Z34 

MMnPanMmaoflSO  SMmI  OavMon  -  ZSS 

Rayiaqft*  0.0000 

Ftguf»A.37.  RMtits  from  stiiisiicai  analysts  usng  Rosy  for 
W«2^  cm  (1.00  inert)  and  S-18.7%  in  trta 
honzontai  plana  (SHP8). 


I  W-254  cm  S^7.4%  HORIZONTAL  (SHP9) 

SWMOe* 

I  N-55a 

I  ClaM  irtMvn  «  10  riaqw 

1  MaamumPwoaraaoa*  10.4 

I  M«fln  Pwowaaqa  .  5.08  StandMd  OMtton  •  Z4a 

Vmr  Mm  «  2718 
Ccfi.Angm*  i3.79 

R  Idapnaim  «  0.241 
nayMi»«  0.0000 

Fgure  A.38.  ResuJts  from  staostical  analysts  using  Rosy  for 
W-2.54  cm  (1 .00  inert)  and  S»37.4%  in  ana 
honzontal  plana  (SHPB). 
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Rgur*  A^.  RmUIs  from  statiscicil  anoiytm  using  Rosy  for 
W«2^  cm  (1 .00  incn)  and  S-56^  in  tna 
honzontu  plana  (SHPB). 


W.2.54  cm  3-74.9%  HORIZONTAt  (SHPB) 

Sttatddoa  \ 

N-510 

CUm  (rttiM  -  10  dsgrM 

Msamum  Parowsags-  yQA 

Mm  Pmmag*  -  S.sa  Standard  OcvMtan- i.9S 

Vaew  Mm  -  275.8 
Ctay.An9ta- 20.29 
RMagnduda-0.174 

nayiinn- 0.0000 

Fgurs  A.40.  RaauJts  from  statistical  analysis  using  Rosy  for 
W-2.S4  cm  f  1 .00  incn)  and  S-74.9%  m  (tie 
honzontai  plana  (SHPB). 


W.2.54  cm  S^3%  VERIICAL  (SHP8) 

StBMDCS 

N«S44 

V«cnrM«an«27da 

OaM  tflimi  •  10  dagmas 

Cord.  Angia  •  9v43 

Usnmum  PvaMags  •  13^ 

RMagr«uda«0246 

Mm  Psforaoi  •  19tt  Standard  Omdgn.  122 

Raytwtfi*  0.0000 

F^urt  A.41 .  Rtstits  from  jtaiisticai  analysts  using  Rosy  for 
W-2^  an  (1 .00  inch)  and  3-0%  in  tf)a 
vart)cai  qlana  (SHPS). 


W«2.54om  S-18.7%  VERTICAL  (SHPW 

SMMtQB 

N.73S 

1  CUM  Intsrva  •  10  dagrMS 

1  MMmjmPwoamgsa  i02 

Mm  PMosmaof  S.0O  SlanMrd  Omdon  -  lOS 

V«cMrMm.d82 

Cord.  Angia*  10.70 
RMaontiidf02e8 
RavMgn  *0.0000 

Ftgure  A.42.  Results  from  sia&socai  analysis  using  Rosy  for 
W-2.54  cm  (1 .00  inch)  and  S-1 8.7%  »n  me 
verticai  piane  (SHP0). 


t 
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W«2.Man  S^7.4%  VERTICAL  (SHPq> 

N-4fl9 

Ctesimwvit-  lOdigw 


N-4S9 

Vaeur  Mm  «  273^ 

CtesinMvit-  lOdigw 

CoiS.  Angie*  11.25 

Mrnnum  Pwcsrmse  •  14J 

R|y^«ude-0.33S 

)49m  -  5.30  SMidvd  OmSbh  •  3J2 

RayW^aOUXXW 

Ftgur*  A.43.  RmuIIb  from  sadstical  anaiytit  usvig  Rosy  for 
W-2^  cm  (1 .00  inert)  and  S-37.4%  in  ms 
vsracii  plans  (SHP6). 


W-2.5*cm  5-i6&2%  VERTICAL  (SHPS) 

StaMOoB  1 

N.632 

CteM  UMivM  •  to  dtgreie 

McnmumPoraeniaoe*  10.4 

Mean  Pereamaoe  *  5.SS  SttndeidO«n«ion-2.15 

Vaenr  Mm  •  272.1 

Cofd.  Angie*  13.S7 

R  Megnsirte  ■  0238 
Reytai^  *  0.0000 

Ft^ure  A.44.  Results  from  siatisacsJ  analysis  using  Rosy  for 
W-2.S4  cm  (1 .00  inen)  ano  S-56.2%  in  me 
vertical  plans  (SHPB). 


W.2.5icRi  Sm7A.9%  VERTICAL  {SHPm 

SttOMIOS  ! 

N-a07 

OSM  immri  •  1 0  dooTM 

MaxmimPsroMaga-  iU 

MMnP«t»niaoB«S.M  SMMd  OMton  •  2J0 

Vsenr  Maan  •  2812 
CofS.Angl»»  12.07 

R  Magnsuda  ■  0209 
«  0.0000 

Figur*  A.45.  RMcjits  from  satisfical  anafysit  using  Rosy  for 
"HmZjSA  cm  0  *00  inch)  and  S»74.9%  in  thd 
varbeai  piarw  (SHPS). 


W«2.5«cm  S-0%  SEAMORCUAA  (SHP8) 

SUMUea  ! 

N«  190 

CUaa  iraanral «  i0  diBfaaa 

MaxmumParoaniaQa*  10.7 

Uaan  Pwoanuioa  •  9.56  SMndvd  OaviMon  •  277 

VacHrM«an«0l20 
CofCAngla-  159.58 

R  UaQMiida  «  0.030 
Raytai^  •  0.7813 

Fguro  A.46.  Rasults  from  statisocal  analysts  using  Rosy  for 
W•^54  cm  (1 .00  tncn)  and  S^%  in  m# 
samt-circuiar  piano  (SHP6). 
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W«2.5*CIW  S»ta.7%  SEMK:WCU>a  (NOSHPq) 
N.320 


CtaM  inwva  «  10  degma 

Celt.  Angti  •  40.78 

Maonum  P«raanaoo  •  ar 

niitagMudo«0.107  | 

M««iP«Mniage.9.9a  SMndvd  1.78 

Raytaiofi- 0.0293  | 

f^uf  AA7.  RMtits  from  scttitticil  mmy%m  using  Rosy  for 
W«^54  cm  (1 .00  inch)  and  S-t  8.7%  in  (ho 
som-orcuiar  piano  (SHP6). 


W^^9»aT^  3^37.4%  SEMM^IRCULAB  (SHP« 
N«218 

CtaM  Inlarvfli «  10  dtqwM 
Maaomum  PwcmtmQm*  11.9 

M— n  PtwramsQt  *  5.56  SMndsid  0«Hadon  «  ^42 


V«r  Umm  -  278.9 
CoiS.An«l»*  139.59 
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Figure  A.48.  Results  from  statistical  analysis  using  Rosy  for 
W-2,54  cm  (1 .00  inch)  and  S-37.4%  m  the 
semt-orcuiar  plane  (SHPB). 
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Ffqxjtm  A.49.  RMiits  from  stilisiicai  anatyw  using  Rosy  for 
W«Z54  cm  (1 .00  Inch)  and  S-56^  in  m« 
MfTi-orcUar  pim  (SHPB). 
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ngure  A.50.  Rasults  from  staiumcaJ  anaiysis  usmg  Rosy  for 
W-a.54  cm  ( 1 .00  incn)  and  S-74.9%  in  tne 
sami-arcuiar  plans  <SHPB). 
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F;gur«  1 .  RmUIs  from  aiatitticai  analysis  using  Rosy  for 
W-2^cm(1.00  inch)  and  S-0%  in  layar  #1 
(SHPB). 
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Rgure  A.52.  Results  from  statisiicai  analysis  using  Rosy  for 
W-2.54  cm  ( 1 .00  incn)  and  S- 1 8.7%  in  layer  #1 
(SHPB). 
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Figuro  A.53.  RMuits  from  ot«»ticat  anofyiia  usmg  Rooy  for 
W-2.54  cm  <1 ,00  inch)  and  S^7.4%  in  layar  #1 
(SHP8). 
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Figure  A.54.  Resuits  from  statuocai  anaJysis  using  Rosy  for 
W-2.54  cm  (1 .00  incn)  and  5-56.2%  *n  layer  #1 
(SHPB). 
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Ftgur*  A^.  Rmuh  from  statisbcai  analytta  using  Rosy  for 
W-Z54  cm  (1 .00  incn)  and  S«74.9%  in  laysr  #1 
(SHP8). 
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Figure  A.56.  Results  from  staosocal  analysts  using  Rosy  ^or 
W-2.54  cm  (1 .00  incn)  and  S^%  in  layer  #2 
(SHPB). 
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Rgur«  ^S7,  Rmuis  from  stitisocii  aniiysis  usin^  Rosy  for 
W-2-54  cm  (1 .00  incft)  and  S«1 8,7%  in  tayar  #2 
(SHPB). 
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Figure  A.58.  Reauits  from  statisocai  anaJyaia  using  Rosy  for 
W.2.54  cm  (1 .00  men)  and  S«a7.4%  m  layer  #2 
(SHP8), 
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F»gijf#  A^,  RmUIs  from  atidsiicai  anaiysM  using  Hoay  for 
W«Z54  cm  n.OO  indi)  and  S^S^in  layar  #2 
(SHPB). 
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Ftgure  A.60.  Rasuits  from  stata6cai  anaiysis  using  Rosy  for 
W-2.54  cm  { 1 .00  incn)  and  S-74.9%  m  layer  92 
(SHPB). 


W«2.5i  cm  LAYEH  a  (SKP81 
N«  174 


CtSM  lM¥«l «  10  digieM 

CofS.  Angle-  182 

Manmum  FsreMaee  «  14.9 

R  Megntiirti  •  0.397 

MMPmmaga.a^  aMvd  OMtfon  •  XOS 

naylltfi- 0.0000 

Ftqur#  A-61 .  RMirits  from  stadtocai  analysts  using  Rosy  for 
W-2^cm  (1.00  incn)  and  5-0%  in  iaysr«3 
(SHPB). 
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Figure  A.62.  Results  from  stattsocal  analysis  using  Rosy  for 
W-2.54  on  ( 1 .00  incn)  and  S- 1 8.7%  tn  layer  #3 
(SHPB). 
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Rgur*  A.63.  RMtiti  from  stuiaociri  anaiysa  usmg  Rosy  for 
W.2^  cm  (1 .00  ifKn)  and  8-37.4%  in  layw  »3 
(SHP®. 
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Figure  A. 64.  Reeuits  from  statistical  analysis  using  Rosy  for 

W»2.S4  cm  (1 .00  «ncn)  and  S-56.2%  *n  layer  #3 
(SHPB). 
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Ftgur9A.65.  RMtitt  from  ataosbcari  anaiysit  us«>g  Rosy  for 
'HmZJU  cm  (1 .00  inch)  and  S«74.9%  in  iaysr  #3 
(SHPB). 


Rgure  A. 66.  Results  from  staosocai  analysis  using  Rosy  for 
W-2.54  cm  (1 .00  incn)  and  S-0%  in  layer  #4 
(SHPB). 


Rguf»  A.67 .  RMuits  from  staiitocai  aniiystt  u«ng  Rosy  for 
cm  (1 .00  inch)  and  Sat  8.7%  in  iayar  «4 

(SHP8). 


Figure  A. 68.  Reauits  from  staiiaocai  anaiysia  us»ng  Rosy  for 
W-2.54  cm  (1 .00  incn)  and  $-37.4%  (n  layer  9A 
(SHPB). 
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Hgur*  A.69.  Riauiti  from  stuitticii  an«iyBt6  uunq  Rosy  for 
W-2^  cm  (1.00  tncO)  and  S-56^  in  layw*  #4 
(SHPB). 
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Figure  A.70.  Reeuits  from  statistical  analysts  using  Rosy  for 
W-2.54  cm  ( 1 .00  ifx:n)  and  S-74.g%  m  layer  #4 
(SHP8). 
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